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THE SPECTRUM OF NOVA SAGITTARII 1954 
(HARO-HERRARO) 
M. W. Feast 
(Communicated by the Radcliffe Observer) 


(Received 1955 July 25) 
Summary 


A general description is given of the spectrum of Nova Sagittarii 1954 
(Haro~Herraro) from 1954 July 7 to August 20. Strong C1 lines are present 
in the principal absorption spectrum which has a velocity of 729 km/sec 
The hydrogen lines show several absorption components whose velocities 
are listed in Table Il. ‘The structures shown by the broad emission lines 
are discussed A distance of 3 kps and an absolute magnitude at maximum 
of about —7™-8 are deduced from the equivalent widths of the interstellar 
Catt lines 
Introduction and observations.—News was received at the Radcliffe 
Observatory on 1954 July 7 of the discovery on 1954 July 4 of a 7th magnitude nova 
(17" 50™-5 — 36° 15’ (1875)) in Sagittarius. Several spectrograms were obtained 
that night and the nova was observed from time to time till 1954 August 20. The 
nova had then become sufficiently faint that further spectroscopic work would 
have caused too much interference with other programmes. ‘The spectrum of 
the nova underwent considerable changes during the period covered by the obser- 
vations and the present paper is intended to give a summary of these changes. 
The observations were all made with the 2-prism spectrograph at the Casse- 
grain focus of the 74-inch Radcliffe reflector. Between July 7 and August 5 the 
{/3°7 camera (49 A/mm at Hy) was used and after the latter date the f/2 camera 
(86 A/mm at Hy). ‘The spectrograms obtained are listed in ‘Table I. 
Generally Kodak 103aB plates were used and the focus was adjusted to the 
optimum value for the blue and violet spectral regions. ‘Though in such cases 
the red region is not in the best possible focus, the lose in resolution is small and 


the plates are usable from Ha to H{ (Ha ts actually outside the normal working 


range of 103aB plates but its great strength in the nova ensures its appearance on 


the spectrograms). 

Each spectrogram was measured for wave-lengths of absorption lines as well as 
widths and positions of structural features in the broad emission lines. It will be 
found simplest to discuss the absorption and emission spectra separately though 
in fact of course the two are closely related. !t will also be convenient to use the 
nomenclature suggested by McLaughlin (1, 2) and to compare the present results 
with the hypothetical “typical” nova as described by McLaughlin. ‘Though 


33 





462 M. W. Feast Vol. 115 


rather rough magnitude estimates were generally made at the times of taking the 
spectrograms, it is not intended to discuss the light curve or its relation to the 
spectrum since accurate photoelectric measures of the nova are likely to become 
available in the near future (Dr A. B. Muller of the Leiden Southern Station, 
Johannesburg, informs me that he has obtained a good series of measures). In 
the present paper the time of development of various features is given so that when 
a light curve is published it should be relatively simple to relate it to the spectrum. 


Taste I 
List of Radcliffe Spectra of Nova Sagittarii 1954 


Date (U.T.) Emulsion Camera Date (U.T.) Emulsion Camera 
1954 
July 7 17° 30 «©103aF ig" July 19 20" 51 103a0 fl49-7 
7 «(37 Zs fs 21 36 103aF - 
17 4! ‘a : 22 o9 pe 
17 48 Ls se 22 20 2 103aB 
20 40) = 103a0) ; se 6S 
21 30 103aF , . § 55 ” 
21 55 a - 8 20 40 - 
23 «55 ” ‘ 10 20 626 
9 21 oc 103aB 13 42 
16 ©O 43 ie - 16 20 30 
18 22 24 - - 20 20 40 


The spectrogram on July 16 was kindly obtained by Mr H. van Woerden 
The spectrograms on July 18 and August 1o were kindly obtained by Dr A. D. Thackeray 
* spectrograph adjusted for best focus in red region 


The principal absorption spectrum.—When the nova was first observed on 
1954 July 7 the spectrum showed rather weak emission but a strong set of absorp- 
tion lines. ‘This absorption spectrum persisted (with some modifications) for a 
month or more whilst the emission spectrum gradually increased in intensity. 

On 1954 July 7 the main absorption spectrum matches moderately well a super- 
giant F type star except that all the lines are rather diffuse. ‘This set of absorption 
lines is pretty clearly the principal spectrum and the lines in it were identified 
without difficulty. ‘The compilation of Cecchini and Gratton (3) was very help- 
ful in this respect. ‘The stronger lines listed by Cecchini and Gratton in their 
Table IX are found in the principal absorption spectrum of the nova. The 
main contributors to the spectrum are H, ‘Ti, Fen, Cru, Sem, Mgu, Sri, Sit, 
Caut,C rand O1, Whilst many of the lines are blends, about 14 (other than H lines) 
were selected as being free from any serious blending and were used for radial 
velocity determinations, though not all the lines were available for use on each 
spectrogram. Over the period from 1954 July 7 to August 5 the velocity was 
constant within the accuracy of the measures, ‘The mean velocity was — 729 km/sec. 
It should be noticed that the hydrogen lines of the principal spectrum were 
not used in deducing this velocity, for reasons given below. 

An outstanding feature of the principal spectrum is the presence in consider- 
able strength of the lines of C1. The lines AA 4762, 4772-76, 5039-42, 5052 and 
5380 are all strong. ‘The absence of A 4932 is probably due to the overlying Fe 1 
4924 emission, ‘The lines appear to be of the same order of intensity as in the 
R Corona Borealis variables of types F-G which are held to have an anomalously 
high abundance of carbon in their atmospheres. ‘The presence of C 1 absorption 
in nova spectra was first established in Nova Herculis (DQ) 1934 (4) and 
McLaughlin (§) has shown the lines to have been present in several earlier novae. 
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In view of this it has been suggested (for instance by Bidelman (6)) that novae are 
in some way related to carbon stars. However, it is not certain that strong C1 
absorption lines are always present in nova spectra and they were apparently 
missing from Nova RR Pictoris 1925. 

Unfortunately none of the spectrograms is strong enough at 3883 A for the 
presence of CN to be detected (as it was in Nova Herculis). The band at 4216 is 
not detectable on the Radcliffe spectrograms but this is a less sensitive test for CN. 
The general similarity of the principal spectrum to that of the hot R Cr B variables 
suggested that weak C, bands might possibly be present in the nova spectrum. 
Whilst it is just possible that the bands at 5165 A and 4737 A are weakly present, 
no very definite conclusions can be drawn from the available spectrograms. 

Fairly strong absorption lines of O1 were observed in Nova Herculis (4) and 
several other novae (§) but notin Nova Pictoris. In the present nova the presence 
of Ot absorption is not quite so clear as with the C1 absorption. Probably 
the O1 blend 6158 accounts for a line which would otherwise remain unidentified. 
016455 is probably blended with Fe tl 6456. ‘The Tit line 4368 interferes with 
the detection of O1 4368 but the line appears to be stronger in the nova than in 
n Sgr (cF2) so that a contribution from QO 1 is indicated. 

The principal absorption spectrum is observed from July 7 to August 5 apart 
from H which persists longer. A few other lines may also persist after August 5, 
but, due to the fading of the continuum and the lower dispersion used after 
this date, they are not very clearly seen. Whilst the velocity remains sensibly 
constant throughout, various changes are noted in the principal absorption spec- 
trum. In particular the lines all gradually sharpen and there are also various 
changes in relative intensity. ‘The most noticeable change is the rapid fading of 
Mgt 4481 and the increase in prominence of various lines of ‘Tit. Both these 
phenomena have been well observed in other novae and are probably connected 
with the dilution of the radiation in the shell as it moves away from the star. It 
does not appear that it would be very profitable to study the changes in detail on the 
present spectrograms due to the difficulties that would be encountered with the 
large number of blends present. 

The hydrogen absorption components.-—It now remains to discuss the various 
absorption components of the hydrogen lines. ‘The various components that have 
been distinguished are denoted by roman numerals for convenience in discussion, 

The hydrogen lines associated with the principal spectrum (1) are observed 
throughout the entire period July 7 to August 20. However, unlike the rest of the 
principal spectrum, the velocity of the H lines is not constant. At the start of the 
series on July 7 the velocity agrees fairly well with the rest of the principal spectrum 
(~729 km/sec) but the lines gradually increase their displacement towards the 
shorter wave-lengths, reaching —88qkm/sec on July 22. ‘The displacement then 
gradually decreases and is again close to —729 km/sec at the end of the series 
(August 20). In view of the constancy of the rest of the principal spectrum it 
appears unlikely that this is due to a real velocity change. ‘The shift is more 
probably caused by an unresolved component (I1) which has a greater negative 
velocity shift than I, and which gradually rises and then falls in intensity, being 
weak or absent on both the first and last spectrograms of the series. 

Also present on July 7 is II] with a velocity of ~ 1700km/sec. On July7 I is the 
strongest absorption but III gradually increases in strength and is stronger than I 
on July 18. III becomes rather broad and between July 7 and July 19 shows a 


/ 
progressive shift towards shorter wave-lengths, reaching —1895km/sec. This is 
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apparently the effect of blending of components of varying relative intensities 
since on July 22 the broad line is seen separated into two components, at 
1775 km/sec and —2053km,sec. ‘The former is considered as a continuation 
of II whilst the latter is a new component IV. By August 5 both III and IV 
are weak and, with the continual fading of the continuum and the change to a 
lower dispersion camera, they are not seen again after this date. On August 5 a 
new component V is seen for the first time with a displacement of — 1164 km/sec. 
It is possible that V was actually present at an earlier date but blended with I. 
By August 8 this component has increased in intensity but is not subsequently 
seen, due presumably to both true fading and to weakening of the continuum. 

Whilst it appears that I can be assigned to the principal spectrum and II, IIT 
and IV are in all probability part of the diffuse enhanced spectrum, which is 
known to be frequently multiple, the place of V in McLaughlin’s scheme is not 
quite so clear, V is probably another component of the diffuse enhanced spectrum 
though its assignment to an Orion absorption phase is not impossible. 

As noted earlier the principal spectrum (1) is quite rich in lines and a consider- 
able number of elements have been identified in it. On the other hand, apart from 
the strongest lines of Fe 11 which are present in absorption III no lines other than 
those of H have been identified with certainty in any but the principal spectrum. 
It should, however, be borne in mind that such lines could easily be masked by 
lines of the principal spectrum. ‘The absorption lines of N 111 (4097, 4103) are 
frequently seen in the later stages of nova development. Whilst displaced N 111 
has sometimes been suspected in the present nova, the continuum near 4100 A is 
weak at the time when these lines would be expected to be strong. 

‘The velocity measurements in the nova are summarized in Table II. 

The emission spectrum.—In this section a summary is given of the lines identi- 
fied in the emission spectrum and their changes in relative intensity and structure. 
As with the absorption spectrum the emission lines are sufficiently close to the 
“typical’’ nova for it only to be necessary to point out the principal features of 
the spectrum. 

Below are listed the main lines identified at one stage or another of the nova 
development : 

H : a, B, vy, 0, «, | 

Het: 5876, 4471 

Cit: 4267 

NI: §942~§2, 5667-86, 5530-44, 5045, 5005-4994, 4774-4803, 4447, con- 

tribution to 4640. 

Nill: 4097, 4103 

[N tt] : 5755 

Or: 6158 

Oit: 4415, contribution to 4640 

[O1] : 6364, 6300, 5577 

[¢ ) rt} :§007, 4959, 4363 

Nal: 5890, 5896 

Catt: 3986 (blended with He) 3933 

Fett: 6456, 6248, §317, 5276, §235, 5197, §169, 5018, 4924, 4583, 4233, 4174 

The following lines are of doubtful identification : 

6004: A fairly strong line through the series of spectrograms. It has also been 
found in some other novae. Its identification is discussed by Larsson—Leander 
(7) who found the line to be present in Nova DK Lacertae. The identification 
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with NI 5999, 6008 is considered by him to be unsatisfactory in view of the absence 
of the other N1 lines. This view is confirmed by the observations of the present 
nova. 

5814: A diffuse band. ‘The variations in the intensity of this line suggest that 
itis part of the Orion spectrum (see below). It is possibly C 1v 5802—12 which is 
found in Wolf—Rayet stars and has been suggested (8) as a possible identification 
for a line in the Orion spectrum of Nova Herculis. The identification should, 
however, be considered rather doubtful. 

A broad emission from 4030 to 4063 ; Contributions from O tt and N tare possible. 


Tassie Il 


Velocity Measurements in Nova Sagittarii 1954 


Absorption | 729 km/sec 
I! (<I 
ill 1700 km/se« 
IV 2053 km/sec 
V 1164 kin /sec 


Interstellar Catt velocity 6 km /se« 


A broad band stretching from 4565 to 4497 gives the impression, especially on 
the last spectrogram of the series, as consisting of two overlapping emissions with 
centres at approximately 4515 and 4539. ‘These may be 4510-4514 Nitiand 4541 
Heit. (4686 Hel is probably present.) Fe 4508-4549 may also contribute. 

In the main essentials at least the development of the emission spectrum 
follows that shown in the 


‘ 


‘synthetic spectra”’ of a “ typical’’ nova as constructed 
by McLaughlin (2). The first spectrum obtained is probably a little earlier than 
McLaughlin’s “ spectrum”’ for 1™-o below maximum and the last is between his 


sé 


™-6 below maximum. 


spectra’’ for 3-8 and 4 
The actual drop in magnitude during the period covered by the observations 
appears, however, to have been greater than is indicated by these figures 
The chief changes found may be listed as follows : 
(1) ‘The gradual fading of Fe 11 and Cau and the increase of [01] [Om] [N11] 
Het, and C1. 
The moderately rapid rise to maximum about August 5 of diffuse lines of 


N 
— 


Nu. This is the Orion spectrum. 
(3) The rapid development of the wide 4640 Orion emission about this time and 
its persistence throughout the rest of the seri 
(4) ‘The presence of wide H and Fe lines of the diffuse enhanced spectrum 
which overlie the principal emission spectrum lines from early in the seri® 
but which fade out with the fading of the diffuse enhanced absorption 
McLaughlin (g) recognizes various rather definite stages in nova development 
characterized by the change in relative intensities of various lines. An attempt 
has been made to estimate these stages in the present case. ‘They are: 
(1) The [O1] flash July 18 
(2) The [Nu] flash July 26 or slightly earlier 
(3) The Helium flash July 26 
(4) Cul 4267 Feil 4233 =1 about August 13 
(5) First trace of 5007 [O 111] August 13 
(5) [O11] 5007/Fe 1 4924 = 1-0 between August 13 and 16 
Later stages are not covered by the present series of spectrograms 
On the first spectrograms of the nova the emission spectrum is weak, even the 
hydrogen lines being relatively inconspicuous. Hy and 4 have the appearance 
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of broad hazy bands extending from the principal absorption line (1) (at 

729 km/sec) to about + 1050km/sec with broad maxima at +630km/sec and 

200 km/sec and a minimum at + 100km/sec. ‘The relative intensities of these 
features seem to change along the Balmer series. ‘The maximum at — 200 km/sec 
strengthens from H4 and y to Hf and «; in fact, Hf and « appear to show this 
maximum alone as a moderately sharp line displaced —2zo0okm/sec. It seems 
unlikely that Feit 4351 can markedly affect the appearance of Hy since other 
Feu lines are very weak at this stage. The interpretation of these lines 
is not quite clear as the well-known profile consisting of the H lines of the 
principal spectrum on which are superimposed the wider and weaker lines 
of the diffuse enhanced emission spectrum does not develop properly till between 
July 9 and 1s. 

By July 15 the emission lines have become strong and thereafter the inter- 
pretation is relatively simple. ‘The diffuse enhanced emission never becomes 
very strong and the main features, apart from the Orion phase shown mainly by 
Nut, are the wide but fairly sharply defined lines of the principal spectrum. 
The H lines of the principal spectrum, as seen about July 15, show two broad 
maxima at about +415 km/sec and — 420 km/sec, the shortward peak being the 
stronger. A minimum is seen between the two peaks. 

The main development in these bands from July 15 onwards is the gradual 
fading of the shortward peak and its further displacement towards the short wave- 
length edge of the line (to about —620km/sec) together with the strengthening 
and sharpening of the longward peak, till from August 8 onwards the hydrogen 
lines have the appearance of broad uniform bands with a strong line superposed 
at +415 km sec. 

Since the majority of the lines observed are in the visual region of the spectrum 
where the dispersion is rather small, especially with the f/2 camera, it is not possible 
to discuss the structures of these lines in any detail. And in any case a detailed 
discussion could only be attempted after accurate spectrophotometry. However, 
it may be said that all the lines (apart from the diffuse lines of the Orion spectrum) 
tend to show structures similar to the hydrogen lines and to develop peaks approxi- 
mately 800 km/sec apart. Nevertheless all the lines do not have the same structure, 
as is well illustrated by the lines 4267 Cu and 4471 Het. ‘These lines increase 


rapidly in intensity towards the end of the series and are strong on the last spectro- 
gram obtained (August 20). Since these lines are in the blue region their struc- 
ture is more easily studied than is the case with the majority of the lines. These 
lines appear as broad uniform bands with fairly sharp peaks near each edge (at 


+ 423 and —570km/sec). The peaks with positive velocity displacement agree 
in position with the prominent peaks shown by the H lines but in 4267 and 4471 
this peak is considerably weaker than the peak of negative velocity displacement, 
the latter being present, if at all, only very weakly in the H lines. 

The interstellar Ca 1! H and K lines.— Strong sharp interstellar Cat H is seen 
throughout the series of spectrograms. Cat K is only seen near the beginning of 
the series, as later the fading of both the nova Catt K broad emission and the 
continuum leave no background in this region of the spectrum. ‘The persistence 
of broad He emission ensures a background on which Cat H absorption may 
be seen. A mean velocity of —6km/sec has been derived for these lines from the 
various spectrograms. The constancy of this velocity as well as the general 
appearance of the lines leaves little doubt that not only are they interstellar in 
origin but also that they are not affected by blending with nova features. 
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The position of the nova is sufficiently close to the direction of the galactic 
centre for the galactic rotation effect to be negligible. ‘Thus a distance for the 
nova cannot be derived from the radial velocity of the interstellar lines. An 
attempt has therefore been made to estimate the distance of the nova from their 
intensities. Microphotometer tracings of the H and K region were obtained 
from two of the most suitable f 3-7 spectrograms and equivalent widths were 
derived. Spot sensitometer marks were used to calibrate the photographic plates. 
The mean of the two plates is: 

Equivalent width Cait H o-7 A 
Equivalent width Cau K 11 A 

The ratio of the equivalent widths is therefore K/H=1-6. ‘This ratio is in 
good agreement with the observations of Sanford and Wilson (10) who find that 
with increasing equivalent widths the ratio approaches 1°56, being practically 
equal to this value for H>o-4 A. 

Using the relationship between the intensity of the K line and distance, re- 
cently derived by Beals and Oke (11), a nova distance of about 3000 ps is deduced, 
If a general space absorption of o”-8/kps is assumed the distance modulus is 
M —m-=14"™-8, so that if the nova reached 7™.0 at maximum (it may perhaps have 
become even brighter) the absolute magnitude at maximum was — 7:8. In view of 
the considerable uncertainties involved in this type of calculation the value is in very 
satisfactory agreement with the absolute magnitude of other novae. McLaughlin 
(12) deduced a mean absolute magnitude of —7™-0 + 0-2 for novae but since he 
used, amongst other data, the absolute magnitudes of novae in both the Magellanic 
Clouds and the Andromeda nebula which are now known to be in error by about 


1™ in view of the new distance scale (13), his mean required changing to about 
_ i. @ 


7 °5 
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AN OCCULTATION OF ANTARES OBSERVED AT ‘THE 
UNION OBSERVATORY 


David S. Evans 
(Communicated by H.M. Astronomer at the Cape) 


(Received 1955 January 28)* 


Summary 


The reappearance of Antares on 1953 March 5 was observed by Union 
Observatory and Leiden Southern Station personnel using the Rockefeller 
Twin Photographic Telescope of the Leiden Station, Union Observatory 
Johannesburg. A photoelectric trace of both components was obtained 
The angular diameter of the principal star is approximately 0o”'o39. An 
auditory record was also obtained on magnetic tape in the form of a con- 
tinuous note, frequency-modulated according to the output of the photocell 


Personnel.—'The work described in this paper was done at the request of the 
author, and the analysis of the observation and preparation of this report were 
undertaken at the request of the Union Astronomer. ‘The work was organized 
by Dr W. S. Finsen; Dr Th. Walraven was responsible for the photometer and 
associated amplifiers; Mr J. Hers for the frequency-modulation recorder ; 
Mr P. C. Seligman for the Hathaway recorder; Mr J. Churms acted as time-keeper. 


Equipment.—The telescope used was the Rockefeller ‘Twin Photographic 
telescope of 16 inches aperture, and focal length 2:25 metres. ‘The guiding 
telescope was of 6 inches aperture with a micrometer capable of being offset by 
about one degree. ‘The focal plane diaphragm was of 4 mm diameter, and behind 
it was placed a simple lens of 1 inch focal length. ‘The photocell was an R.C.A. 
931 A operated at go volts per stage with go volts on the collector stage. A filter, 
Schott BG 5, which, apart from infra-red transmission, has a maximum trans- 
mission of about 14 per cent at about 4300 A was employed, for the following 
reasons: first, to reduce the secondary spectrum and to permit the use of a small 
diaphragm (the peak transmission of the filter corresponds closely with the 
wave-length for which the object-glass of the telescope was corrected); second, 
in preliminary experiments it was found that with no filter and full aperture 
the light from Antares produced overshooting of the recording galvanometer. 
In these circumstances the initial deflection was not maintained, but fell rapidly 
over a period of several seconds, finally settling down to a value about half as 
great as the initial one. Application of a d.c. voltage to the amplifier produced 
no such effect. When the filter was introduced, the effect was much reduced, 
but there was still a change in deflection with time after every strong change in 
light, in the reverse direction from that initially seen (i.e. an undershoot instead 
of an overshoot), somewhat slower, and of a much smaller magnitude. 

The explanation of this behaviour has not been found: it was considered by 
the Union Observatory team to be due to excessive light from Antares, and, 
for the given photocell, and in the given circumstances, there may be validity 


* Received in original form 1954 August 8. 
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in this explanation. On the other hand, exactly parallel tests with much higher 
fluxes have been made on the occasion of other observations of Antares occultations 
without showing any behaviour of this kind. All that has ever been observed 
when the cells have been deliberately over-tried on very bright stars has been the 
normal progressive type of fatigue, in which there has been a loss of response 
of a few hundredths of a magnitude per minute, and only slow recovery. The 
incident flux for the full light of Antares can readily be computed from Whitford’s 
formula (1) as 1-7 x 10°-* lumens for photoelectric light. ‘The photocell current, 
with the filter in, was o'1 microamperes. Even allowing for a reduction of as 
much as 5 magnitudes by the filter, if the values for the unobscured light of 
Antares are brought into comparison with data for the performance of 931 A 
photocells (2), it can be seen that even extremely slight, very slow, fatigue ts 
only just possible, while saturation phenomena should only have occurred with 
fluxes tens of thousands of times greater. 

What concerns us now is the validity of the observation. It is considered 
that, although the necessity for correcting the remaining instrumental drift 
introduces an element of personal judgment, an approximate value of the diameter 
of Antares can be derived from it. 
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To return to the equipment: the output from the photocell could be fed to 
three recording devices. ‘The first of these was a Brown recorder, the feed being 
brought through a 100 megohm series resistor to the 10 megohm input of the 
amplifier of this instrument. ‘The second was a Hathaway recorder. The 
output of the photocell was fed to it through a d.c. amplifier, and the output of 
this at low impedance fed to three galvanometers of the Hathaway connected in 
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series and providing a load of about 150 ohms. ‘Two of these galvanometers 
deflected in one direction, the third in the opposite direction. They had 
sensitivities of the order of 10 microamp/mm, and with damping resistors their 
response was flat to 150 c/s (Fig. 1). The galvanometer block was removed from 
the recorder case and was mounted separately to avoid motor vibration. Low- 
speed galvanometers were used to record one-second pips and alternating current 
at 50 c/s from the quartz crystal clock. ‘The frequency of the latter is accurate 
to better than one part in 10’ and at the time of the occultation the clock correction 
was +0°007 seconds. ‘The recording speed in the Hathaway recorder was 
5 inches per second: the internal time-marker was run from the 50 c/s mains 
merely in order to provide lines for ease of measurement later. 
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Brown recorder readings 


Fic 2.—Calibration tests 


The third recording device was the frequency-modulation recorder devised 
by Mr Hers. ‘This was connected in parallel with the three galvanometers of the 
Hathaway recorder. The principle of this instrument consists in frequency 
modulation of a low-frequency carrier, the modulation frequency being a function 
of the photocell output. The record is made on magnetic tape and it has been 
found that when the tape is played back through a suitable demodulator to the 
Hathaway recorder, the result does not differ materially from that obtained 
directly. It is hoped to publish details of the method elsewhere. Its principal 
advantages appear to be: tape economy (for the tape can be left running for 
thirty minutes and can be re-used); wide frequency response; reduced liability 
to accident as compared with many types of recorder; and provision of an audible 
monitor, convenient during the event itself, as well as during preliminary 
adjustments and for demonstration purposes. 
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The equipment was calibrated by comparing the responses of the galvano- 
meters with those of the Brown recorder (assumed linear) to the same source. 
The source was Antares itself immediately after the occultation, reduced by 
various filters. ‘The filters were (in addition to the Schott BG 5): empty aperture, 
Schott NG 3, Schott NG 4, and Schott NG 5. ‘The simultaneous records on 
the galvanometers and on the Brown recorder for these settings, repeated several 
times, enable the overall linearity of the equipment to be tested. Using the 
deflections attained after instrumental drift has ceased, one finds the linearity 
of the three records given by the three galvanometers to be very good (Fig. 2), 
and no correction 18 necessary. 

The observation._-On nights previous to the occultation the micrometer 
was carefully set on the guiding star CPD — 26° 5642, with Antares centred in the 
diaphragm. ‘This adjustment was checked shortly before disappearance and 
again immediately after reappearance, when Antares was found to be accurately 
centred in the diaphragm. ‘The photometric seeing at the time of the observation 
was much better than average and decidedly better than on previous nights. 
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Time (seconds) 
I ; Reduced trace of reappearance of Antare 


The leneths of the vertical lines sh the maximum displacement duc to 
seeing at those levels 


This is borne out by comparison of the traces obtained for the occultation with 
those obtained at rehearsals, the latter showing far greater fluctuations, As 
a check, the diaphragm was placed at such a distance from the bright limb of the 
Moon that the intensity of the sky light approximated to that of Antares itself. 
The photocurrent fluctuations were then much less than in the case of Antares, 
showing that atmospheric conditions are the principal cause of the fluctuations 
recorded. 
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The companion reappeared at 2" 20 11%93 U.T. and the centre of the 
principal star 8-78 later. This last-named figure has already been used in 
combination with other observations to derive a separation of 3°-16 and a P.A. of 
274 38’ for the pair (3). The trace of the principal star as recorded on the three 
galvanometers is shown in Plate 6. ‘The reappearance took slightly over # 
second, so that there is no distortion produced by lag in the galvanometers. 
On galvanometer II, which was used for the reduction, the deflection for Antares B 
was 2°6 mm, and for A, 26 mm. ‘Iwo or three trial reductions from measures 
of the traces of galvanometers Nos. II and III have been made. They indicate 
that the chief source of uncertainty in drawing the final curve lies in the allowance 
to be made for drift immediately after the reappearance of the principal component. 
The extent of this uncertainty in the estimated diameter is three or four units in 
the third place of decimals. ‘The measures shown in Fig. 3, which are reduced 
to percentage deflections, represent what is considered to be the best choice for 
the value of the drift. ‘The rate of advance of the lunar limb normal to itself at 
the point of emergence was taken to be 0"-360 per second. ‘The continuous 
curve drawn in is that for the advance of a straight edge across a uniform disk of 
diameter 0°:039, which is taken as the value of the angular diameter indicated by 
this observation. Diffraction effects are negligible for disks of this order of 
diameter, ‘The numerical differences introduced by taking this effect into account 
never exceed 1 per cent, and are confined to very short sections at the toe and 
shoulder of the curve (4). 


The observation team wish to acknowledge their thanks to the Director of 
the Bernard Price Institute for Geophysical Research, Johannesburg, for the 
loan of the Hathaway recorder, and to Messrs ‘Thos. Barlow and Sons (S.A.) Ltd, 
for the loan of accumulators. 


Royal Observatory, 
Cape of Good Hope 


1955 January 21 
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SURFACE PHOTOMETRY OF THE MAGELLANIC CLOUDS 
A. R. Hogg 
(Communicated by the Commonwealth Astronomer) 
(Received 1955 May 3°) 


Summary 


Two-colour photoelectric observations have been made at a number of 
points in the Magellanic Clouds and results reduced to the (B, V) system 
of Johnson and Morgan. Isophotes have been drawn for B magnitudes 
and sectional diagrams for the variation of colour B-V across the brighter 
regions of the Clouds. 

The observational results give, for certain limited regions of the Clouds 
the following figures : 

5.M.¢ L,.M.C, 
Area sq. degrees (measured) 10°35 27°5 
B magnitude 2°6 14 
B-V colour 0°33 
Central brightness mag. sec™* 21°4 


The integrated brightnesses given are obviously lower limits to the actual 


higures 
The colour of the Small Cloud is found to be systematically bluer in the 
brighter central regions than nearer the edge. ‘The Large Cloud shows no 


regular variation of colour over its surface 


Introduction.—Integrated magnitudes of the Magellanic Clouds have been 
given by G. van Herck (1), by H. Shapley (2), by J. 5. Paraskevopoulos (3), by 
k. Hubble (4), and by G. de Vaucouleurs in an unpublished revision of the 
Shapley-Ames catalogue (5). 

The magnitudes on photographic scales are as follows : 
5.M.C 


van Herck 4 2m.g 
> 


Shapley +2~3 
Paraskevopoulo (a) 
(b) 


Hubble < bors 
de Vaucouleurs ob 


and show a scatter reflecting the difficulties surrounding the measurement of such 
objects. 

The results given by van Herck were obtained by extra-focal photometry and 
are probably the most satisfactory of those so far published. ‘The figure given by 
Shapley was derived from star counts, whilst the results attributed to Paras- 


kevopoulos were obtained by comparison of in-focus images of the Clouds with 
| | 


* Received in oriinal form 19 January 19 
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extra-focal images of bright stars. With a lens of 1in. focal length Paras- 
kevopoulos’ result was 1™-8 whilst with a lens of } in. focal length the result was 


2-3, ‘The figures mentioned by Hubble are described by him as “ estimates 
based upon published descriptions ”’ 


The present investigation aimed at defining isophotes in the Cloud regions by 
scanning with a photoelectric cell to determine a magnitude at least for the central 
regions of the Clouds. Partly because of the small acceptance cone of the photo- 
meter and the large angular dimensions of the object, this method becomes very 
tedious if applied to the detailed examination of the peripheral regions of the 
Clouds; and this has not been aimed at in the present series of results. 

Observational._-The observations were made at the Cassegrain focus of a 
20-inch reflector. Stray light from the sky was excluded by means of an iris 
diaphragm mounted on an extension of the draw tube. A diaphragm in the 
focal plane defined an area of sky approximately 3-4 minutes of arc in diameter 
and admitted light via a Fabry lens to the cathode of an EMI photomultiplier 
(Type 5060). ‘The photocurrent was amplified and fed to a Brown recording 
potentiometer, Filters employed were b= 3mm BG 18 and y=2mm GG 5. 
Apparently dark regions of the sky adjoining the Clouds were used as control 
areas and measures were made of the magnitude difference between the control 


areas and regions of the Clouds. ‘The control regions chosen were as follows : 


for the 5.M.C 4(1950)=02" oom 6(1950) 
L.M.¢ ee 
20 


40 


“Iw J 
www uv 


In the case of the L.M.C, the mean of the three areas was used as a control. 

Ihese control areas were also compared with the star F of the Harvard F1 
region (6), the magnitude of which had been found by comparison with a number 
of stars (particularly 40 Leonis, which is very nearly the same colour as F) for 
which magnitudes have been given by Johnson and Morgan (7) (UBV system) 
Magnitude—colour relations were so obtained to convert the “ natural”’ system of 
filters and cell to the UBV scale. ‘This gave for F a result on the B scale of 7™-g1 
and a colour of +0™-33 on the B—V scale. ‘The magnitude of the control areas 
by comparison with F amounted to 


S.M.C. control area 12™-58 per sq. minute 
L..M.C. control area (mean of 3 positions) 12™-67 per sq. minute 


from measures on five and eight nights respectively ‘These magnitudes are on 
the B seal quoted above 

Magnitude measures and isophotes.-Measures with the blue filter (6) were 
made at 113 points in the 5.M.C. on 8 nights (total 397 observations) and at 
84 points in the L.M.C. on 7 nights (total 302 observations). In setting the 
telescope for these measures care was taken to exclude any obvious bright stars 
from the field The limit of such stars was probably 12™—13! ‘The measures 
were partly traces at constant « with varying 4 and partly traces at constant 4 
with varying «. Results were arranged in order of 6 and by interpolations 
appropriate curves were drawn to define isophotes (Fig. 1). The central 
regions of the Small Cloud give average total readings (sky + Cloud) o™-8 brighter 
per unit area than the control area, i.e, they are 2-09 times the brightness of the 


‘ 
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control area which, as mentioned above, is 12™-58/min®. The brightness due 
to the Cloud alone (assuming that the control area truly represents the 
background) is 1-09 times that of the control area, i.e. 12™.58 —o™-og 12™-5/min* 

21™-4/sec* (5.M.C.). The corresponding figure for the Large Cloud is 
12™-67-0"-19 = 12-5 min® = 21™-4/sec? (L.M.C.) Uhese figures are on the 


B scale 


O of2c* O10 


LMC 


Fic. 1 lsophotes in the Magellanic Clouds (B system) The contours refer to the magmtud 
difference between the control regvions and the total light + d from object hy 


This result shows the measured parts of the central regions of the two Clouds 
to have about the same surface brightness, although the bright central portion 
of the Small Cloud covers much less area than that of the I irye The isophote : 
give some indication of the westward extension of the Small Cloud often referred 
to as the “ wing”. 

Integrated magnitudes.—Integrated magnitudes were obtained by measure- 
ments with a planimeter of the areas included by the isophotes in Fig. 1 Partly 
because of the extended area of the object some difficulty was experienced 
in locating the o™-o isophote. Accordingly the following arbitrary limits of 
measurement were adopted : 


S.M.C. (1950) « =o00"25™ —orh2ch 0 18 sq. deg 
LL.M.C. (1950) « = 04"40™ — o6"00™ 
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In the case of the Small Cloud the limits seem to include most of the obviously 
brighter portions of the object, but with the Large Cloud the adopted limits 
excluded some obviously bright but irregular areas on the northern edge. ‘The 
magnitudes given by this process, viz. 

S.M.C, = + 2™-6 
L.M.C, = +14 


are therefore to be regarded as lower limits to the total brightness of the Clouds 
especially in the case of the Large. 
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Pia, 2.--Colours in the Magellanic Clouds (B-V system). The figures alongside the points 
indicate the number of observations and the vertical lines through the points denote standard 
deviations for those cases where more than two observations are availabl 


The regions measured are notably less than the recognized dimensions of the 
Clouds, e.g. H. Shapley (8), gives angular diameters of 8° and 12° for the Small 
and Large Clouds respectively and the latter figure has been increased to 20° by 
G. de Vaucouleurs (9). Furthermore the measured integrated brightness will 
also be affected if the adopted control area does not represent the true sky 
background, i.e. the sum of (a) the amount of light from terrestrial sources, (4b) the 
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amount of light coming from stars between the observer and the Cloud and (c) 
the amount of unobscured light coming from behind the Cloud. This is 
obviously a difficult quantity to measure with certainty and is not necessarily 
represented by readings taken around the edges of the Cloud. Some readings 
were obtained which indicated the presence of regions which were some one or 
two hundredths of a magnitude fainter than the control areas although in the 
case of the Large Cloud the mean of several peripheral regions showed zero 
magnitude difference against the control on the V system whilst the B result 
showed this mean to be brighter by o™-o1 than the control. If the adopted control 
areas were actually o™-o2 brighter than the true background then the integrated 
magnitudes given above would need to be decreased by o™1. 

Whilst there was no reason to suspect that the present photoelectric results 
were greatly in error (the integrated magnitudes being in reasonable agreement 
with van Herck’s figures) it seemed desirable to undertake certain check observa 
tions. 

Therefore measures were made, after the isophotes had been drawn, by 
observations on sections through the brighter parts of each Cloud. In the case 
of the Small Cloud new measures (O) at 13 points were made with the same 
apparatus and compared with results (C) interpolated from the diagram. ‘This 
gave a mean O —C = —o0"-00 + 0°02 (s.d.). 

\ new photometer (EMI cell Type 5659, filter combination 1 mm Chance OBI 
glass+2mm Schott GGis5 glass) with a different photometer head having an 
aperture of 3° was used to obtain similar check. readings on both Clouds. 

For the Small Cloud O—C for 18 points was +o-o1 + 0-02 and for the Large 
Cloud for 10 points O—C = +0°04+0°03. The grand mean of the O —C for 41 
check points is therefore a little in excess of one hundredth of a magnitude. In 
view of the small-scale irregularities in the Clouds this can be regarded as a 
satisfactory result. 

Colours.—The colour of the Small Cloud was determined mainly from two 
N-S and E-W traces centred near the brightest area using 70 observations on 
five nights. ‘The figures were brought to the (B—V) scale of Johnson and Morgan 
by comparison of 13 stars. The results clearly show that the Small Cloud 
exhibits a definite gradient of colour, being bluer towards the centre than at the 
edges. ‘The same effect is also shown in Table I, where the mean colour is 
determined for regions of varying brightness. The standard deviations and 
general drift of the means show that this is not a chance result. 

Colour variations across a galaxy (the M31 nebula) have been detected by 
Stebbins and Whitford (10), and also by Fricke (11), but in this case the colour 
increases towards the centre and also exhibits a marked intensification in a dark 
urea (see Fricke, Fig. 17) Ihe colour of the Small Cloud at the centre is notably 
bluer than the results for other irregular nebulae quoted by Holmberg (12) 
Qualitatively a discrepancy of this sort could arise from an unsuitable choice of 
background region In the present case it is found that, in the unlikely « it of 
the true background being fainter than the corrected background by as much 
o"-o2 in the yellow, the colour should be o™-os bluer than indicated ; a correction 
which does not improve the discrepancy. 

The distribution of colour in the Small Cloud is difficult to account for by 
absorption due to dust intermingled with the stars in any regular fashion Th 
distribution suggests that the hotter stars are concentrated towards the centre 


34 
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of the Cloud. The overall colour of the Small Cloud does not support the idea 
of a Type II population, which according to Baade (13), would be much redder. 
This fails to confirm a preliminary result obtained by Gascoigne and Kron (14) ; 
a result which is now being withdrawn by the authors (private communication). 


Tasie I 


Colours B-V of the Magellanic Clouds compared with surface brightness in magnitudes 
brighter than control areas 
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brightness B-t 
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The colour results for the Large Cloud are arranged according to surface 
brightness in Table 1, which summarizes 155 observations on 7 nights at 45 points. 
Traces across the Cloud at « = 05" 40" and 6= — 70° are shown in Fig. 2. This 
figure also displays the number of observations used for each point and when this 
exceeds two, the standard deviation is indicated by lines. Although the trace in 
declination suggests a slight reddening at — 70°:5 the standard deviations are such 
as to make it unlikely that the result is real. ‘There is also no systematic change of 
colour with the surface brightness (‘Table 1) and it seems appropriate to represent 
the results by a grand mean giving a colour of +0™-33 for the Large Cloud; even 
though one area at «(1950) = 5"40™ and (1950) =69°-5 was found, where B-V 

+o™12, The average colour transferred to a photoelectrically defined Inter- 
national Colour Scale (7) is +0"-19, which is rather bluer than the result given 
by Holmberg for the resolvable irregular nebulae ( + 0™-3). 

Conclusion.— The foregoing results clearly indicate a distinction between both 
Magellanic Clouds and other irregular nebulae for which colours have been found. 
Whilst it is not impossible that despite the relatively high galactic latitudes of 
the Clouds ( — 32° for the Large and — 45° for the Small) some galactic absorption 
occurs, it would not qualitatively affect the conclusions that the Small Cloud is 
bluer than a Type II population and that the Large Cloud is also bluer than the 
usual resolvable irregular nebula. ‘The colour of the central portion of the Large 
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Cloud is somewhat, and that of the central portion of the Small Cloud is 
considerably bluer than the colour of other irregular nebulae. Further, the 
systematic variation of colour across the Small Cloud contrasts with the practical 
uniformity of colour over the Large. 


The assistance of Mrs A. de Vaucouleurs in measuring and reducing the 
records is gratefully acknowledged. 


Commonwealth Observatory 
Canberra, A.C.T 
1955 April 27 
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THE PERIOD-LUMINOSITY RELATION IN POPULATION II 
V. C. Reddish 


(Communicated by Professor W. M. H. Greaves) 


(Received 1955 April 5) 


Summary 
Correlations between the periods, luminosities and radii of Population II 
variables with periods between one and a hundred days are investigated 
It is shown that the recently suggested division of the period—luminosity 


relation into three discrete lines (1) probably cannot be justified on the basis 
of existing data 

A further analysis of the published observations on twelve variables (1, 2) 
leads to the conclusion that two of the variables differ considerably in com- 
position, though not in mass, from the remainder, and have therefore 
probably evolved along different lines. 


The results of accurate observations of the colours and magnitudes of thirteen 
Population I] variables have recently been published by H. C. Arp (1). ‘The 
variables have periods in the range one to a hundred days. In an analysis of the 
data, Arp suggests that the period—luminosity relation consists of three discrete 
lines, each differing from its neighbour by a factor of two in period, and 
representing different modes of pulsation. 

In this paper the analysis of the observations is extended, and considered in 
conjunction with pulsation and stellar structure theory. It is shown that there is 


probably no justification, on the basis of existing data, for supposing that these 
variables pulsate in different modes. It appears that two of the variables (those 
with the shortest periods) differ considerably in composition, though not in mass, 


fromthe remainder. ‘The remaining variables specify a unique period—luminosity 
relation. 

The period-luminosity relation* is reproduced in Fig. 1 on an M,,, scale, 
together with Arp’s lines. The data is given in Table I, The mean M,, are 
those given by Arp (1). Mean spectral types are from observations by Joy (2). 
Luminosity class and effective temperature from the Table in Hynek’s 
Istrophysics (3); bolometric corrections by Kuiper (4). 

Arp noted that in four of the variables successive cycles differed but alternate 
ones reproduced. He suggested that the true period was therefore the time 
for two cycles. The effect of doubling the previously used periods is shown 
on Fig. 1 by dashed arrows. 

The period of pulsation of a variable is related to its mass and radius by the 
relation 

P(M) R*)"? = 0-041 (1) 
where the period is in days and the mass and radius in solar units, the constant 
being that derived by Epstein (5). 


* Only the twelve stars for which spectral types are known are plotted, since only for these 
can reasonably accurate bolometric corrections be made (see next footnote) 
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The period luminosity diagram according to Arp The circléd dots are those for which 


Arp doubled the pe riod, as shown by arrows and crosses 
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Sources of the data are given in the text 
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From the point of view of the structure related to pulsation, a radius—period 
diagram is therefore more fundamental than a period-luminosity relation. Such 
a diagram is given on Fig. 2 from the data in Table I.* The radius increases 








0 








——————— 


0 0.5 ‘10 ter 1S 


Fic. 2.—Radius-period diagram. The circled dots are those for which Arp used double the period, 


continuously with period. It therefore follows from equation (1) that there is no 
discontinuity of mass along the sequence. 
Fig. 2 gives the relation 


log R = 0°74 log P+ 0°75. 
Equation (1) may be written 
log R = 0°67 log P + 0°33 log M + 0°925. 
i:quations (2) and (3) lead to the relation 


log M=o 21 log P oO’e4. 


he radii have been obtained using surface temperatures derived from the spectral types 
If the colours are used to give surface temperatures, the resulting relation similar to equation (4) 


predicts masses a mall as o-o1rM for the brightest variables Considerations of age and 
evolution show that this is ridiculous. The colours result in too high surface temperatures, the 


error increasing with luminosity. The use of spectral types gives much more agreeable results 


such a difficulty 
is a result of insufficient knowledge of colour—temperature 


It is not suggested that the measured colours are in error; the possibility of 
occurring was pointed out by Arp; it 
relations 
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For log P=o, P=1 day, this gives M=o-3, and for log P=1-5, P=30 days, 
M-=o0'6. These values are only about, or less than, half those to be expected 
according to recent estimates of the masses of the brighter Population II stars (6). 
This discrepancy (if in fact it is) could result from the spectral types being too 
early by 0-1 to 0-2 spectral classes or (and more probably) from over-estimating 
the luminosity class. Since the masses of Population II stars are considerably 
less at a given magnitude than those of Population I, their densities at a given 
luminosity and spectral type are less; thus those stars given luminosity class II 
in Table I may perhaps be Ib.* Such a change would lower the calculated 
effective temperatures and increase the bolometric corrections; both these 
changes would increase the calculated radii and hence the masses. 

Theories of stellar structure give relations involving mass, luminosity, radius 
and composition. A radius—luminosity diagram is therefore of interest, and is 
shown on Fig. 3. All but two of the stars fall close to a straight line. ‘The 
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Radwus-—luminosty diagram The circled dots ave those for which Arp used double the 
pe riod 


radius and luminosity vary continuously along the remaining sequence of stars. 
Since the mass also varies continuously there is no discontinuity in composition 
along this sequence. Wowever, since the masses and radii of the two stars lying 


off this sequence are continuous with those on it, they must differ from the remainder 


in luminosity and hence in composition. Yor a given radius and period, they are 
brighter by about three-quarters of a magnitude than the sequence formed by the other 
ten stars. ‘They deviate from the line by six and eight times the standard deviation 
of the remaining points: such large displacements are unlikely to be due to 
chance. 


* Popper (7) classified some of the brighter stars in Vii and Mii as I 
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It is now instructive to re-examine the period—-luminosity diagram, given 
in Fig. 4 (the points are the same as on Fig. 1). The two stars which do not 
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Fic. 4 Peritod-luminosty diagram The circled dots are those for which 
Arp used double the pe riod 


form part of the continuous sequence on Fig. 3 are the two of shortest period 
Neglecting these two points a straight line is drawn through the remainder. 
The standard deviation of the ten points about the line is o™-22. Much of the 
error probably results from the calculation of mean luminosities from the observed 
magnitudes, uncertainty of the zero point of the magnitudes in each cluster, 
and the conversion to bolometric magnitudes. ‘The dispersion is about the value 
that should be expected (1). The distribution of the errors is in good accord 
with the law of distribution of random errors, as can be seen from ‘Table II 


this supports the view that the ten points define a single linear relation, 
M1 0°07 — 2°67 log P. (5) 


Tasie Il 
Deviations of the ten points from the straight line on Fig. 4, compared to those to be expected 
from a normal error distribution 
x = deviation Observed number Expected number 
in magnitudes with deviation 2x with deviation > x 
o'4 i 
3 
4 
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The two stars of shortest period are again about three-quarters of a magnitude brighter 
than the sequence defined by the remaining stars. 

It is evident from Fig. 3 that these two stars differ from the remainder not, 
as Arp assumed, in period as a result of a different mode of pulsation, but in 
luminosity as a result of a difference in composition. ‘The distribution of the 
remaining ten points gives no justification for a division into different modes of 
pulsation. ‘The differences of successive cycles, and similarity of alternate 
ones, which were noted in four of the variables, may result from interaction between 
the fundamental and other frequencies rather than from high excitation of the 
second harmonic. 

The opacity in such high luminosity, low density stars will be due to electron- 
scattering (unless the abundance of heavy elements in their envelopes is very high). 
The luminosity then depends on the composition in proportion to ut/(1 +X) 
where » is the molecular weight in the envelope and X the hydrogen abundance. 
For a factor of two in luminosity this implies a difference of about 15 per cent in 
molecular weight and a difference of between 10 per cent and 30 per cent in 
hydrogen abundance. ‘Thus the two stars of shortest period probably have about 
20 per cent less hydrogen in their envelopes than the remaining ten stars. 

Such a difference in the composition of the envelopes may be supposed to have 
existed throughout the lifetime of the stars, or to have developed comparatively 
recently. ‘The first alternative suggests that the luminosity L has always been 


* and since age is proportional, for a given stage 


higher by a factor of at least two, 
in the evolution, to M/L where M is the mass, the age would be only about half 
that of the remaining stars. However, if a group of stars forms with considerably 
different composition and at a time different from the majority, it is extremely 
unlikely that they would evolve at such a rate that at some particular later time 
those which had reached a certain stage in their evolution would have the same 
masses as the remaining stars which had also reached that stage It is much more 
probable that the two stars have only recently strayed from the path of evolution 
followed by the rest, and that the change in the composition of their envelopes 
has occurred fairly quickly Ihat such a splitting of the evolutionary paths may 
occur has been predicted by the writer (8), but of course the existence of two 
stray stars cannot be considered as strong evidence of this; they merely point 
to the need for more observations 

On his final period-luminosity diagram, Arp included six variables in w 
Centauri for which the values of m,, were known. Since the mean spectra 
ire not known, it is not possible to include them tn Fig. 4. It appears likely 
however, that four of them will fall close to the straight line, while the other 
two may be close to the pair which are three-quarters of a magnitude brighter 
than the line 
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* This f e may not be correct if there are differences in heavy element abundances, since 
the earls fe of the irs opacity may be due partly to photoelectric absorption; but arguments 


imular to that which follows will still hold 
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POLARIZATION MEASUREMENTS ON THREE INTENSE 
RADIO SOURCES 


R. Hanbury Brown, H. P. Palmer and A. R. Thompson 
(Communicated by A, C. B. Lovell) 
(Received 1955 May 2) 


Summary 


The radiation from three intense extra-terrestrial radio sources, in 
Cygnus, Cassiopeia and Taurus, has been tested for polarization at a 
wave-length of 1-gm. ‘The results show that, in the flux from the sources 
in Cygnus and Cassiopeia, there is no plane polarized component greater 
than 1 per cent and no circularly polarized component greater than 4 per 
cent. In the case of the source in Taurus (the Crab Nebula) the corresponding 
lumits are 24 per cent and 4 per cent respectively 


1. Introduction.—At the present time the Sun is the only extra-terrestrial 
radio source which is known to emit polarized radiation. It is well established 


that, under certain conditions, a significant fraction of the solar radio flux is cir- 
cularly polarized. ‘The radiation from the discrete sources of galactic and extra- 
galactic radio emission (radio stars) is generally taken to be unpolarized, although 
there is little evidence to support this assumption. ‘The only published experi- 
ments appear to be those of Ryle and Smith (1) who examined the intense sources 
in Cygnus (1gN4A)* and Cassiopeia (23N5A)* with an interference polarimeter, 
They concluded that, at a frequency of 80 Mc/s, the flux from these two sources is 
unpolarized and they set an upper limit of about 5 per cent to the intensity of any 
polarized component. 

The present paper reports some measurements which were intended to check, 
and if possible improve, the existing limits set by Ryle and Smith, and also to 
extend the observations to the intense source in ‘Taurus (os; N2A)*. ‘his source 
has been identified with the Crab Nebula, and the recent report (3) that the light 
from this nebula is polarized suggested that the radio flux should also be tested for 
polarization. 

2. The Equipment.—An outline of the equipment is shown in Fig. 1. The 
two acrials X and Y were arranged to form an interferometer with a baseline of 
about 200 metres. ‘This distance was chosen so that, at a wavelength of 1-gm, 
the interferometer had a resolving power of about half a degree and was therefore 
capable of separating the discrete sources from the general background radiation, 
The two aerials were circular paraboloids with diameters of 30ft and 25 ft. A re- 
flector with double curvature, such as a paraboloid, is not ideal for polarization 
tests ; nevertheless the present aerials were selected because of their unusual con- 
struction, which allowed each paraboloid to be rotated about its axis of symmetry 
without altering its direction. ‘This particular property made it possible to reduce 
the systematic errors in the measurements 


* "The number wiven awainst then yurce refer to thw (_atalogue ofl Reliab Known 


Discrete Sources of Cosmic Radio Waves ’’ prepared by a commuttee of 1LA.U. (see reference (2)) 
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The 30ft paraboloid (X) was equipped with two identical primary feeds at right 
angles, each feed consisting of a dipole and reflector. The 25 ft paraboloid (Y) was 
equipped in the same way for the detection of plane polarization (Fig. 1(a)); but, 
for the detection of circular polarization (Fig. 1(b)), a single dipole and reflector 
were substituted. In the arrangement shown in Fig. 1(a) the two dipoles in each 
aerial were connected by short cables of equal length to a coaxial changeover switch 
(5,,5,) which in turn was connected to a preamplifier feeding a long cable to a 
central receiver. ‘The two switches were remotely controlled from the receiver in 


Fic. 1 Outline of apparatus for measuring polarization, 


(a) for plane polarization, (b) for circular polarization 


such a way that either of the dipoles in each aerial could be selected at will. In the 
second arrangement, shown in Fig. 1(b), only one coaxial switch was employed. 

The central receiver was tuned to 158-2 Mc/s and had a pre-detector band- 
width of 400 kc/s and an output time-constant of 5 seconds. It operated on the 
“ rotating-lobe ’’ principle and has been described elsewhere (4) ; in effect it multi- 
plied the two signals together and displayed their product on a pen recorder. 

3. Method of observation. (a) The measurement of plane polarization.-F or the 
measurement of plane polarization the apparatus was connected as shown in 
lig. 1(a). ‘The two aerials were fixed in azimuth (due south), and were tilted in 
elevation so that the axes of the paraboloids were directed to the declination of the 
source. ‘The paraboloids were then rotated about their axes so that the primary 
feeds were in the relative positions shown in Fig. 2(a). With the aerials in this 
position the amplitude of the interference pattern produced by the source was 
observed during its transit through the aerial beams, the observations being mad 
alternately with the dipoles AD and BC 
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In order to minimize any errors in the relative amplitude of the patterns obser- 
ved, the receiver was switched between the two aerial combinations (AD,BC) 
every five minutes during the period of observation (Fig. 3). This switching 
procedure eliminated the effects of any slow changes in the receiver gain, such as 
might be expected from day to day. It was found convenient to increase the 
frequency of the interference pattern, by means of a rotating phase-shifter (4), so 
that at least ten complete cycles occurred in each switching interval, thus enabling 


the amplitude of the pattern to be measured with satisfactory accuracy. 
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Fic. 2 Relative positions of the primary feeds in the two aerials as viewed from the source 


The first observations were made on the most intense of the sources, namely 
that in Cassiopeia. After a sufficient number of satisfactory transits had been 
observed with the aerials in the relative positions shown in Fig. 2(a), the two para- 
boloids were rotated through 45° (Fig. 2(b)) and the measurements were repeated. 
Finally a third set of observations were made with the paraboloids rotated through 
go” (Fig. 2(c)). 

The sources in Cygnus and Taurus were observed in the same way, except 
that the measurements were confined to two positions of the aerials 45° apart 


(Figs. 2(a) and 2(b)). 
A f\ \\\ f\ 
| i | | | 


> 


1 simplified diagram of an interference pattern, 


(b) The measurement of circular polarization.-Yor the measurement of cir- 
cular polarization the apparatus was connected as in Fig. 1(b) and the paraboloids 
were adjusted so that the primary feeds occupied the relative positions shown in 
Fig. 2(d). ‘The amplitudes of the interference patterns produced by each source 
were then observed using alternately the parallel dipoles BC and the orthogonal 
dipoles AC; as in previous tests the receiver was switched between these two 
combinations every five minutes. 

4. Plane polarization-analysis and results. (a) Method of analysis.—If the 
radiation incident on the aerial system of the apparatus shown in Fig. 1(a) is ran- 
domly polarized or circularly polarized, then the amplitude of the interference 
patterns (/, 5, /»,,.) observed with the orthogonal pairs of dipoles (AD or BC) will 
be equal and independent of the rotation of the paraboloids about their axes. But, 
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if a significant component of plane polarization is present, then there will be a 
difference in the amplitudes of the patterns given by 


Iyp~! 


Be 

a =m, COs (20) (1) 
p 

lay , ly, 


where m, (assumed to be <1) is the fraction of the incident flux which is plane 
polarized, and @ is the angle between the plane of polarization and that of the 
dipoles A and D. 

The records of each source were analysed to give values of J,» and Jy, for 
each angular position of the aerial system. ‘The upper limit to the value of the 
polarization (m,) was then estimated by means of equation (1). It was assumed 
in making this estimate that the plane of polarization might have any value, but 
that it remained constant throughout the period covered by the observations. 
This assumption is discussed later in Section 6. 

It was assumed in equation (1) that the gain of each aerial is independent of 
which of the two primary feeds is in use, Since the beam formed by each para- 
boloid is elliptical, this assumption is only valid for directions close to the axis. 
Thus, in principle, observations of the amplitudes /,), /,,~ should be restricted 
to the instant when the source transits the centres of the aerial beams. In the 
present experiment this difficulty was overcome by recording the mean values of 
1,» Jo in alternate five minute intervals for about half an hour before and after 
the transit of the source. ‘The amplitudes /,5, /,¢ corresponding to the time of 
transit were then obtained from smooth curves drawn through these observed 
values. 

It was also necessary to ensure that there was no systematic difference in the 
gain of the two aerials with the two orthogonal sets of primary feeds. A difference 
in gain is to be expected if there is any lack of symmetry in the paraboloids or if the 
primary feeds are not identical. ‘The necessary check was provided by the obser- 
vations made on the source in Cassiopeia. It was found that for all three posi- 
tions of the aerials (Figs. 2(a), 2(b) and 2(c)) there was no detectable difference 
between the amplitudes of the patterns observed with the dipoles AD and BC. 
It followed from this result that any difference between the gains of the aerials 
with the two combinations of primary feeds could be neglected 


Taste I 
Results of plane polarization measurements 

Percentage of flux 

Source which is plane 
polarized (m, x 100) 

Cassiopeia (23N5A) I per cent 

Cygnus (19N4A) I per cent 

‘Taurus (og N2A) 2¢ per cent 


(6) Results.—The results of the tests for plane polarization are shown in 
Table 1. ‘They are based on six transits of the sources in Cassiopeia and ‘Taurus, 
and five transits of the source in Cygnus. All the observations were made within 
a period of four weeks in December 1954 and January 1955. ‘The observations 
on Cassiopeia were divided between the three aerial positions shown in Figs. 2(a), 


(b) and (c), while those on Cygnus and ‘Taurus were taken with the two aerial posi- 
tions shown in Figs. 2(a) and (6). 
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No evidence of plane polarization was found for any of the three sources. In 
all cases the difference in the amplitudes of the two patterns (J, p —/,,) was not 
significantly greater than the uncertainty introduced by random noise. ‘Thus 
the limits given in Table I were determined entirely by the signal to noise ratio, 
and to improve them substantially would require the use of larger aerials. 

5. Circular polarization-analysis and results. (a) Method of analysis.—lf the 
radiation incident on the aerial system (Fig.1(d)) is randomly polarized, then an 
interference pattern will be observed with the parallel dipoles BC but not with the 
orthogonal dipoles AC. However, if there is a component with circular polari- 
zation, then a pattern will be observed with both combination of feeds. Thus 
if I.) ge are the amplitudes of these two patterns then 

As SS mM, (2) 
Be 


where m, is the fraction (assumed to be <1) of the total incident flux which is 
circularly polarized. 

The records of the three sources were analysed in the same way as those 
obtained in the previous tests. ‘The resulting values of /,,,and /,., corresponding 
to the time when the source was in the centre of the aerial beams, were then used 
to estimate the upper limit to any circular polarization by means of equation (2) 

It must, of course, be noted that a plane-polarized component can also give 
rise to an interference pattern with orthogonal dipoles. However, this possi- 
bility could be neglected in the present analysis, since the presence of plane polari- 
zation had already been ruled out by the previous tests. 


Taare Il 
Results of circular polarization measurements 


Percentage of flux 
Source which is circularly 
polarized (m 100) 
Cassiopeia (23N5A) 4 per cent 
Cygnus (19N4A) 4 per cent 


Taurus (o5 N2A) 4 per cent 


(b) Results.—No significant evidence of circular polarization was found for 
any of the three sources. The results of the tests are shown in Table I]. They 
are based on two transits of Cassiopeia, three transits of Cygnus and four transits 
of Taurus. All the observations were made between 1955 January 12 and 17 

The upper limits given in Table I] were set, not by the signal-to-noise ratio, 
but by imperfections in the aerial system. In practice it is difficult to ensure that 
the response of each aerial is zero to waves polarized at right angles to the plane of 
the feed ; in consequence a large component of unpolarized flux produces a small, 
but observable, interference pattern even when the dipoles in the two aerials are 
orthogonal. The presence of this spurious pattern sets a limit to the minimum 
fraction (m,) of circular polarization which can be measured. An estimate of the 
errors in the aerial system, together with a comparison of the records obtained on 
the three sources, showed that the minimum value of m, which could be detected 
reliably was about 4 per cent, and hence this limit has been given in Table II 

6. Discussion. The results presented above how that any component of 
plane polarization from the sources in Cygnus and Cassiopeia must be less than 
1 per cent of their flux at 1-gm. They also show that any component of circular 
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polarization must be small (<4 per cent) and therefore confirm the conclusion 
reached by Ryle and Smith (1). 

No evidence of polarization, plane or circular, was found in the radiation from 
the source in ‘Taurus (the Crab Nebula). ‘Thus the suggestion that the radio flux 
from this object, like the light, might be polarized is not confirmed at 1-gm. If, 
as has been suggested by several authors, the radio emission from the Crab Nebula 
is due to fast electrons in a magnetic field, then it seems rather unlikely that a sub- 
stantial percentage of polarization would be observed at this wavelength. Thus, 
in the presence of the necessary magnetic field (~10°* — 10 * gauss), the plane of 
a polarized wave would be rotated many times in the nebula; however, it must be 
noted that this rotation is proportional to the square of the wavelength, and it is 
therefore important to test for polarization at considerably shorter wavelengths. 

In conclusion one peculiar difficulty of the present type of measurement is of 
interest. ‘Throughout this paper it has been assumed that the position angle (@) 
of the plane of polarization remains constant from day to day. However it has 
been shown elsewhere (§), (6) that, in the presence of the Earth’s magnetic field, 
the polarization of an extra-terrestrial radio wave is modified by the ionosphere. 
A very rough estimate (§) shows that, at a frequency of 100 Mc/s, a plane-polarized 
wave may be rotated through as many as 3 complete revolutions. Since this 
rotation depends upon the total electron content of the ionosphere, the apparent 
plane of the wave may vary significantly from day to day and even from hour to 
hour. Not only does this phenomenon complicate the problem of measuring 
the plane of a wave, but it may also reduce the chance of detecting polarization 
with the technique described here. A detailed analysis of the present results has 
therefore been made on the assumption that the plane of a polarized component 
remained constant throughout a single observation but that it was completely 
random from day to day. ‘The results show that the chance that the limits given 
in Table I were exceeded is increased by a negligible amount, and that therefore 
the conclusions reached above are not significantly altered by the effects of the 
ionosphere 
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CENTRE-LIMB VARIATIONS IN THE EQUIVALENT WIDTHS OF 
SOME INFRA-RED FRAUNHOFER LINES 


B. E. 7. Pagel 


(Received 1955 June 14) 


Summa ry 


Using tracings recorded on Mount Wilson with the McMath-—Hulbert 
Observatory duplicate infra-red spectrometer, the equivalent widths of 
nine Fe lines in the neighbourhood of A 15 500 with excitation potentials 
close to 5°60 volts have been measured at angles of emergence from the 
solar surface corresponding to p= 1°0, 0°8, o°6, o's, 0°4, 0°3 and oz, The 
results were compared with those which would be predicted on the basis 
of three photospheric models recently put forward by Aller and Pierce 
de Jager, and Michard. The line profiles were calculated in detail for each 
model with the aid of modern techniques of high-speed computation, using 
the EDSAC at the University Mathematical Laboratory, Cambridge The 
centre-limb variations in equivalent width calculated from the models diffes 
mainly as a result of differences in the first and second derivatives of the 
source function close to the surface. ‘The behaviour of the stronger lines 
observed is best reproduced from the temperature distribution of de Jager 
while that of the weaker lines (except close to the limb) is more consistent 
with those of Aller and Pierce, and Michard. Finally, the observations wer« 
ised for an mductive determination of the temperature distribution; this 
listribution was then corrected to give agreement with the observational 
data. ‘The final model thus obtained agrees closely in the deeper layers 
(7T,20°06) with Michard’s model In the outer layers, the temperature 
gradient is considerably steeper, as in de Jager’s model, but the additional 
feature of a temperature inversion appears very close to the surface 
(tT, 0°015). It is suggested that a model of this type may be consistent 
with observations of the intensity of the continuum at the extreme limb, as 
well as with the observed behaviour of Fraunhofer lines. 


1. Introduction.—The determination of the detailed structure of the solar 
photosphere is of great interest in connection with two separate problems—the 
examination of the physical character of the outer regions of the Sun, and the 
quantitative chemical analysis of the solar atmosphere from measurements of 
Fraunhofer line intensities. Much work has been carried out with the object 
of determining the distributions of temperature, total gas pressure, and electron 
pressure with optical depth in the photosphere, and the models that have been 
put forward may in general be divided into three classes according to the manner 
in which the temperature distribution is derived. One class consists of theoretical 


models in which the temperature distribution is governed by the theory of 


radiative equilibrium (Chandrasekhar, 1950), with or without modifications 
introduced to take into account the occurrence of convective instability at optical 
depths 7, >1-0 and the “ blanketing effect’’ of radiation scattered or reradiated 


35 
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in Fraunhofer lines near the surface. A second class consists of the limb- 
darkening models, in which the source function (and hence the temperature) 
is derived empirically from the observed limb-darkening curves at various 
wave-lengths, combined with the relative energy—wave-length distribution at 
the centre of the disk, by means of trial solutions of the integral equation of 
limb-darkening (Chalonge and Kourganoff, 1946; Pierce and Aller, 1951; 
Peyturaux, 1952). In the third class, finally, we have the empirical models in 
which the temperature distribution is derived from the behaviour of line profiles 
or equivalent widths across the disk (Voigt, 1950; de Jager, 1952; Houtgast, 1952). 
All models assume the value of the solar constant as a boundary condition, and 
the distributions of pressure and electron pressure are usually derived from the 
temperature distributions ; for this purpose, local thermodynamic and hydrostatic 
equilibrium are assumed, together with a chemical composition consistent with 
the observed spectrum (Stromgren, 1940). 

\ detailed survey of the present state of studies of the photosphere is given 
by Minnaert (1953). ‘The derivation of a theoretical model is complicated by 
three factors: the difficulty in calculating the Chandrasekhar mean of the con- 
tinuous absorption coefficient (Chandrasekhar and Miinch, 1946; Aller and 
Pierce, 1952), the remaining uncertainty in the precise effect of the onset of 
convective instability, and the inadequate convergence of the iterations employed 
in calculating the blanketing effect. Hence, at present, it is necessary to find 
the true temperature distribution empirically. Observations of continuous 
limb-darkening at different wave-lengths, interpreted through expansions of the 
source function in exponential integrals as suggested by Kourganoff (1949), 
apparently lead to very satisfactory results for the layers with monochromatic 
optical depths 7,>0°2, but the source function nearer the surface is uncertain 
owing to the difficulty in obtaining limb-darkening measurements closer to the 
limb than = 0-2, where y is the cosine of the angle of emergence. ‘The study 
of line intensities offers supplementary information of great value, since the 
residual intensities at the different points of a line profile are determined by 
the variation of the source function over the entire range of depths visible 
from outside: in particular, the sharp increase towards the surface in the selective 
absorption coefficients of neutral metal lines leads to a high weighting of the outer 
layers (7, < 0°2) in the calculation of the central intensities and equivalent widths 
of these lines. By observing their variation across the disk, it is possible to 
eliminate some of the uncertainties in the amount of the selective absorption 
coefficient and in its precise broadening mechanism that would complicate the 
interpretation of a line profile measured only at the centre of the disk. If, as now 
seems likely, the mechanisms of non-coherent scattering and pure absorption 
predominate in the formation of Fraunhofer lines (Houtgast, 1942; Spitzer, 1944; 
Unsdld, 1948; Houtgast, 1952), the source function in any layer is directly 
related to the temperature through Planck’s radiation formula, though corrections 
must be applied in the cases where scattering predominates. The temperature 
distribution with continuous optical depth 7, may then be deduced if the relation 
between 7, and the selective optical depth 7, is known for each portion of a line 
profile; the most convenient method of obtaining such a relation is to calculate 
the continuous and selective absorption coefficients as functions of 7, in 
a previously given model that is known to be approximately correct. If the ratio 
» of the two absorption coefficients is critically dependent on the assumed 
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temperature, this procedure must be repeated several times before the final 
model can be obtained; if not, » may be assumed to be the same in the final! 
model as in the initial one. 

The results of recent investigations of the infra-red solar spectrum with high 
resolution, using photoconductive detectors, have been described by Goldberg 
(1950). ‘The variation across the disk in the profiles of the Sit lines near A 16 000 
has been investigated at the McMath—Hulbert Observatory by R. M. Page 
(1951), and the variations in the equivalent widths of the CO lines at A 23 000 and 
at A46000 have been studied by G. Newkirk using tracings obtained with the 
Snow ‘Telescope (Goldberg and Miiller, 1953). A full list of solar lines in the 
wave-length range 1°4-2°5 », together with their identifications and equivalent 
widths, has been published by Mohler, Pierce, M@Math, and Goldberg (1953). 
The tracings obtained with the McMath—Hulbert Observatory duplicate infra-red 
spectrometer (Mohler, Pierce, McMath, and Goldberg, 1950) now have a resolving 
power of the order of 60000 (R. R. McMath, private communication), and are 
very suitable for equivalent width measurements, although accurate studies of 
line profiles must still be confined to the wings of stronger lines (E. W. > 250 mA) 
In the present paper, the observed centre-limb variations in some infra-red 
equivalent widths are analysed in terms of three empirical photospheric models 
recently put forward, and a new model, consistent with the observational material, 
is derived. 

2. Observational results.—The plan in this investigation has been to collect a 
series of data on a number of lines formed under similar conditions of wave-length, 
excitation potential, and ionization potential in order that the trends may not 
be obscured by a number of independently varying parameters. ‘The wave-length 
range A15 530-A 15 670 1s ideal for this purpose, as it contains several Fet lines 
with excitation potentials in the neighbourhood of 5-60eV and with equivalent 
widths ranging from 12mA to 350mA, so that their behaviour across the disk 
may be studied as a function of line intensity only; small differences in the 
degree of collision and Stark damping have a negligible effect on equivalent 
width in this range of line intensity. ‘Telluric lines are almost completely absent 
in this wave-length range. 

lhrough the courtesy of Professor R. R. McMath and of Professor I, S. Bowen, 
it was possible for the writer to obtain a long series of tracings of this wave-lengt! 
region (Fig. 1) at seven different points of the solar disk, corresponding to 
/ 1°0, 0°8, 0°6, 0°5, 0° 4, 0°3, 0°2, taken between 1952 December and 1953 February 
with the McMath-—Hulbert Observatory duplicate infra-red spectrometer; the 
tracings were recorded by Mr W. Livingston. ‘The spectra were taken undet 
excellent atmospheric conditions with a slit width of o25A (or ot1em'; 
resolving power 60000) and the effective resolving power, though difficult to 
measure precisely, certainly does not fall far short of this. ‘The calibration curve 
of the recording system is substantially linear apart from a depression of the zero 
(Mohler, Pierce, McMath, and Goldberg, 1950); this depression was com 


pensated during the tracings by permitting an auxiliary beam of light of a suitable 
constant intensity to reach the photoconductive cell. ‘The deflection of the 
records is thus directly proportional to the intensity of incident radiation; an 


inspection of the instrumental profile recorded with light from a mercury discharge 
tube indicates that a negligible amount of radiation is scattered outside the range 
of the normal profile (which will leave equivalent widths unaffected). ‘The 
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wave-lengths and identifications of the lines were taken from the published list 
(Mohler, Pierce, McMath, and Goldberg, 1953) and the dispersion on the tracings 
is very uniform ; its value is 0°620 + 0-002 cm/A (0161 A/mm), and the equivalent 
widths were found by measuring each area ten times with a planimeter. ‘The 
rate of scanning employed was 4-0 A/min. 

The results of the measurements are given in ‘Table I, where the equivalent 
widths in mA, W, are recorded together with the r.m.s. errors, the number of 
tracings on which each result is based (in parentheses), and the apparent central 
residual intensities r,’ (uncorrected for instrumental broadening). The first 
four columns of Table I are taken from the table given by Mohler, Pierce, McMath, 
and Goldberg (1953). The equivalent widths obtained by the writer at the 
centre of the disk are appreciably larger (~20 per cent) than those published 
by the Michigan group of workers. The percentage discrepancy increases 
systematically with line strength, and is probably due partly to changes in 
resolving power and partly to the manner in which the line wings have been taken 
into account in the present series of measurements. ‘The height of the continuum 
on our tracings is about 20 cm at the centre of the solar disk, so that the discrepancy 
is too large to be accounted for by an error in the position of the zero; furthermore, 
such an error would result in the same percentage discrepancy for all lines, It 
is interesting to note that the discrepancy tends to be much smaller in the case 
of the more blended ‘‘B”’ and ‘‘C”’ lines; the reason for this is not clear, but 
it seems to point to differences in the treatment of the unblended line wings. 
Such differences will not affect the degree of variation across the disk in a line of 
given strength. 


i] 


| } 


v Wala 
/\( | “W™ VV [| 


} | } 


Y\ \\ i ‘i an ee , isan 


, : : : aed diesel 
555160 OTS! WISTT! 19542 15 





Fic. 1.—Tracings of Fraunhofer spectrum at 4 15540, taken with the MeMath-Hulbert 
Observatory duplu ate infra-red spectrograph 


An inspection of Fig. 1 shows that a certain allowance must be made fer 
blending in the profiles of nearly all the lines, though the blending is not severe 
in most cases. When only one half of a line is blended, the Michigan workers 
suggest the reconstruction of the blended half by symmetry. However, such 
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a procedure would introduce serious errors if applied to faint lines such as these, 
as it 16 not easy to draw the central ‘‘axis”’ of a line with sufficient accuracy ; 
furthermore, any asymmetry in the effective instrumental profile would introduce 
a systematic error. Since the degree of blending is always so slight that the 
residual intensity of either line wing would exceed o-go in the absence of the other 
line, it 18 sufficiently accurate to neglect “ saturation”’ effects and write 
T= Te, 

where r is the resultant true residual intensity due to two separate line wings 
which would have residual intensities r,, 7, in the absence of blending; this 
formula is strictly valid in the case of blending between a solar and a telluric line 


(Mohler, Pierce, McMath and Goldberg, 1953). Its adequacy in this context 


is confirmed by the calculations of Section 4, since the residual intensity r = ogo 
(corresponding to W=30 mA) turns out to be on the linear portion of the curve 
of growth for yu 20-2. 


BE 3 


mp peed 
F1G, 2.--Observational results ; log W is plotted against yp, where W is equivalent width in mA 
The sizes of the black circles represent the r.m.s, errors. 

The observational results are shown logarithmically in Fig. 2, the r.m.s. errors 
being represented by the sizes of the representative circles. It is clear that all 
the lines increase in intensity out to p= 0-3, while all but the strongest continue 
to increase out to »=0°2. ‘The mutual consistency in the behaviour of ditferent 
lines with similar intensities is very satisfactory; Fig. 2 may be compared with 
a similar diagram given for lines with a wide range of excitation potentials in the 
visible region by Houtgast (1952). 

3. Calculation of equivalent widths in a given model.-Owing to the high 
excitation potentials of the Fe1 lines studied here, it is appropriate to assume that 
the re-emission is governed by local thermodynamic equilibrium or “ pure 
absorption’ (Unséld, 1948; Pagel, 1954). The equation of transfer may then 
be expressed in the form 

dI{r,, AA, 1) 
4 


} i; "(7 ,, MA) HT Ar,, \A, uw) — B,(r,)}, (1) 
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where /(7,, \A,) is the intensity of radiation at a given point of a line profile 
distant AAA from the line centre in a layer with monochromatic optical depth r,, 
and B,(7,) is the Planck emission function for wave-length A corresponding to the 
temperature 7(7,). ‘lo determine the line profile in a given model, it is necessary 
to calculate the functions B,(r,), 7(7,, AA) for several different values of AA, and 
then to integrate equation (1) subject to the boundary condition 


dl, 


-O as 
dr 


A 
For a given yu, such an integration leads to the surface intensities /,(0, AA, «) for 
the line and /,(0, ) for the adjacent continuum ; the equivalent width is deduced 
by integrating the degree of absorption 
R(AA, 1) r( MA, 2) LAH) (0, AA, 1) (3) 
1 ,(0, j2) 7 
over the line profile. 

Over a limited range of optical depth, it is possible to represent the functions B, 
(AA) by simple algebraic functions of 7,, and this method of * schematic models” 
has been applied by Houtgast (1952) to the analysis of line profile variations 
in terms of a number of different photospheric models In this investigation, 
however, the observational material consists of equivalent widths, the values of 
which depend on the behaviour of B, 7 over a wide range of optical depths. It is 
therefore preferable to use numerical procedures exclusively, employing high- 
speed automatic techniques to deal with the large amount of computing involved. 
The integrations of the transfer equation were carried out using the EDSAC at the 
University Mathematical Laboratory, Cambridge, tables of 7,, B,, and (AA) 
being prepared in advance by hand computation 

The data published by the authors of the various models were first reduced 
to standard form by computing tables of 7, (for A15 500), ©, log P, and log P, 
corresponding to the published temperature distributions and to a standard 
chemical composition with a hydrogen-metal ratio 4=10* and a negligible 
proportion of helium Reasons for preferring this composition are given by 
Michard (1953); its precision is not crucial here, as the electron pressure is very 
nearly cancelled out in the calculation of » for neutral metals, while the lines 
observed are not intense enough for their equivalent widths to be affected 
appreciably by changes of the order of 50 per cent in the collision and Stark 
damping constants (cf. Section 4). The source function B was then calculated 
at each tabulated interval from Planck’s formula, which leads to 


1°3327 x 10! 


[ - so Stieiin 2 
10° wOle 


Bi ssool) erg cm 


sec cm sterad. (4) 


The calculation of 7 may be divided into the following steps 
a) The continuous mass-absorption coefficient» is found from the equation 
j { 


! 
K,= {Pa )+a,(H)}, (5) 


my 


where my is the mass of the hydrogen atom, a,(H_ ) is the total cross-section of 
H- per neutral H atom at unit electron pressure, and a,(H) is the atomic cross- 
section for absorption by neutral H (Chandrasekhar and Miinch, 1946). The 
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values of a,(H~) were obtained by interpolating for A15 500 and for the relevant 
values of © in the tables of Chandrasekhar and Breen (1946); a,(H) was calculated 
from the formula 

aH) =6-211 x 10 “Or! 19 196 (6) 
deduced from the standard formula (Unséld, 1938). The correction for 
stimulated re-emission was not included in these calculations, as it applies equally 
to the continuous and selective absorptions. ‘The ionization of hydrogen does 
not affect the opacity appreciably in the relevant photospheric layers. 

(b) The fraction of Fe atoms in the neutral state was calculated from the Saha 
equation, assuming an ionization potential of 7-90 volts (Moore, 1952). The 
relevant partition functions in the solar atmosphere have been tabulated by 
Claas (1951). 

(c) The relative number of atoms in the lower state of a line with excitation 
potential 5°60 volts was calculated from the Boltzmann equation. Since the 
statistical weight of the lower state varies from one line to another, and since 
the / values for the relevant transitions are in any case unknown, the ratio n,/g, 
of the occupation number to the statistical weight was calculated as a function 
of r, rather than n, itself. The effect of possible deviations from the Boltzmann 
and Saha equations is discussed later (Section 6) 

(d) The Dopple r width AAp was calculated from 


AAp * Fir +#): (7) 


where K is the gas constant, M the atomic weight of Fe, and £, the r.m.s. component 
of small-scale isotropic turbulent velocity along the line of sight. We assume 
€,= 1°7 « 10° cm sec 
(e) In order to obtain the profile of the selective absorption coefficient with 
due regard to damping, it is necessary to obtain the broadening parameter 
r ww T 
(8) 


on 47 Avy ss 4c AAs” 


where I is the total damping constant, due tonatural, collision, and Stark damping, 
and Avy is the analogue in frequency units of AAp. Formulae for Py and 
V'scarx, 48 deduced from the Lorentz—-Weisskopf theory of line broadening, are 
given by Aller (1953), and these are used here with the modification that the 
van der Waals interaction constant ( is replaced by the difference C, — C’, between 
the C's of the upper and lower states (Unséld and Weidemann, 1955). It has 
also been assumed that the a’s for all the lines may be calculated with sufficient 
accuracy from the formulae deduced for the transition &D, — fD,° at Ars 542°15. 
The calculations of Toy and Iggy, details of which are given elsewhere (Pagel, 
1954), lead to the formulae 
log 4.o, = log P + 0-7 log © — log AAy — 6°80, (9) 
lO Agra, = log P, + § log © — log AApy — 3°30. (10) 
Natural damping contributes at the most a few per cent to the total value of a and 
it has been neglected. ‘The Stark interaction constant C, — C, has, in the absence 
of more specific data, been assumed equal in order of magnitude to the value 
deduced by Aller (1953) from the experimental data published by Panter and 
Foster (1937) for the line A 5162-28 (which has the largest Stark effect observed in 


the visible region) 
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(f) The atomic cross-section for selective absorption at an interval 
AA = uAAy 


from the line centre has been calculated as a function of a, u by D. L. Harris 
(1948). From Harris’ tables one may deduce the ratio « x, where « is the relevant 
atomic cross-section for absorption and «, is the “ fictitious absorption coefficient ”’, 
given by 

wt zp? I 


ty= —— f — =1'2010° 8 —— cm* (11) 
mic Avy Mp 


atA1s 5ooif AApis measuredin A. ‘The value of « «, was found at each tabulated 
optical depth for the following values of AA: 0°00, 0°05, 0°10, 0°15, 0°20, 0°25, 
and o-30 A. In the case of two models calculated without damping (Section 4), 
the value AA=0°30 was omitted. The variation of AAy with depth is so gradual 
that «, may be treated as constant. 

(g) If the f values for the relevant atomic transitions were known, it would 
now be possible to calculate 4(AA) as a function of optical depth for each AA. 
However, neither theoretical nor experimental f values are available for Fe 1 in this 
wavelength region and we have therefore calculated relative absorption coefficients 

10°) 


(12) 


VE Mo’ 
where N is the number of Fe atoms per gram of photospheric material. ‘The 
integrations were carried out for a series of assumptions as to the ratio 9 /n, ascending 
by powers of 2 from »=(1/16)n to y= 4n. In this way, curves of growth were 
developed for each » in the various models. 

To carry out the integration, the equation of transfer was rewritten in its 
finite-difference form 


h= AJ = ——! (1-—B)Ar, (13) 
iL 


where the suffixes of equation (1) have been dropped, and use was made of the 
Runge—Kutta subroutine in the EDSAC library (Wilkes, Wheeler, and Gill, 
1951)*. In the Runge-Kutta process, an approximate increment k in the 
dependent variable / corresponding to the whole interval Ar in the independent 
variable is calculated from equation (13) four times, using values of /, B, » 
calculated once for the initial point, twice for the mid-point, and once for the 
end-point of the interval. The final value of A/ is a weighted mean of the four 
k's (Scarborough, 1950) 

A] = §(k, + 2k, + 2k, + hy), (14) 


and the accuracy of the process corresponds to that of Simpson’s rule, though 
the version in the EDSAC library is slightly improved. In accordance with the 
boundary condition (2), the starting value /(r,) is equal to B(r,) if 7, is sufficiently 
large; 7,=4 gives sufficient accuracy in the continuum, and smaller values may 
be used if large values of » occur, as is shown below 

To carry out the computations, the tabulated values of 7, with the 
corresponding values of B and n(AA) were punched on a data tape; the values of 
B, n corresponding to a series of pre-assigned integration intervals in 7 were then 


* The writer is greatly indebted to Professor D. K. Hartree for pointing out the advantages 
of using this subroutine as against attempting to evaluate the integral of the transfer equatiwr 
directly 
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found in the EDSAC by Lagrangian interpolation using Neville’s process 
(W. E. Milne, 1949). ‘The results of the interpolations were used as the initial 
data in the integrations, » being obtained from n by multiplication by a suitable 
powerof2. ‘The values of «.—1-0, 0°8, 0°6, 0-5, 0°4, 0°3, 0-2, 0: 1 —-were successively 
supplied to the machine on a separate data tape, and in this way a “‘ limb-darkening 
relation”’ for one AA was built up in about 20 minutes of machine time. The line 
profiles were obtained by repeating the procedure with the n’s corresponding 
to a new AA, and curves of growth were developed by repeating the procedure 
but doubling the ratio y/n. ‘The equivalent widths were obtained from the line 
profiles by means of Weddle’s rule. 

The accuracy of computations of this kind is ultimately limited by the accuracy 
with which the constants B(r), n(r, AA) are calculated. ‘These constants are 
known at best to an accuracy ott per cent (bearing in mind the compensation of 
errors in the damping when the profiles are integrated to obtain the equivalent 
widths); we therefore aim at a similar accuracy in the computed equivalent 
widths. ‘This requires an accuracy of a few parts per thousand in the computed 
intensities /(o, AA,y). Errors in programming, machine errors, and round-oft 
errors were eliminated by standard methods. ‘Truncation errors, inherent in 
the numerical procedure, may arise from the following causes : 

1. The starting value / = B was taken at a finite optical depth 7, instead of at 

~ a8 implied by the boundary condition (2). 

2. The values of B, n for the half-intervals were found, and the integrations 
carried out, by third-order numerical procedures. 

The effect of these approximations may be found by considering the solution 
of the transfer equation for the case p= 1-0, 7=0. We have 


K(o) = Kr," + |  B(r)e *dr. (15) 


Ir,)=e"( Bir) *dr. (16) 


For sufficiently large values of 7, we may write 
B(r)=a+hbr; (17) 


e.g., in de Jager’s Model VII (de Jager, 1952), equation (17) holds for 7, >34°5 
with a= 4°91 x 10! erg cm? sec cm sterad., b= 0°17 x 10" erg/cm?*/sec/cm/sterad. 
Hence, from equations (16), (17), 
I(7,)— B(r,)=6b 
and the error in /(o) due to the approximation /(7,) = B(r,) 1s, from equation (15), 
jlo) he 7 
Che relative error is thus 
51(0) o17e” aw. 
I(o) 41 


ifr,=4°0. In the general case, this condition is replaced by 


" dr 
| (1+) — ogre (18) 


The starting point of the integration was chosen by the programme itself in 
accordance with equation (18). 





No. 5, 1955 equivalent widths of some infra-red Fraunhofer lines 503 


‘The second approximation inherent in the numerical procedure is harmless 
if the intervals of integration are chosen to be sufficiently small. In addition 
any error in / that is present at optical depth 7, is attenuated according to equation 


(15) by a factor exp | (1 +%)dr «} by the time that the surface is reached, 
“0 


Hence the highest accuracy is needed in the layers closest to the surface. Ideally, 
the lengths of the intervals should be reduced to compensate increases in 1 
and 7, but to do this would make the programme too complicated to be of practical 
utility. ‘Thus the intervals near the surface were chosen to be very small 
(\7r=0-01), while larger intervals were permissible at greater depths (Ar = 0-5) 
The whole range 4:0 7-0-0 was divided into 25 intervals at convenient points 
and the accuracy checked by repeating some of the integrations with different 
intervals of the same order of magnitude. ‘The difference in /, was found to be 
at the most o-1 per cent for «= 0-2. 

4. Results of the computations for three recent models and comparison with the 
observational data.-The procedure described in Section 3 was applied to 
photospheric models corresponding to the temperature distributions put forward 
by Aller and Pierce (1952), de Jager (1952), and Michard (1953). 

Aller and Pierce use the continuous limb-darkening data at various wave- 
lengths to deduce the distribution of temperature with the Chandrasekhar mean 
optical depth r and give tables of log P, log P,, and « (the Chandrasekhar mean 
absorption coefficient) corresponding to a series of different chemical compositions 
lheir model was reduced to standard form by interpolating in their tables to 
obtain the values of log P, log P,, « corresponding to a composition with log A = 4:0 
and no helium; the values of 7, (for A 15 500) were found from the relation 
[ KAY) (1g) 
0 x(t) 

, being found from the tables of 
Chandrasekhar and Breen (1946). In the outermost layers (r<o'1), the 
temperatures were lowered in accordance with the results of blanketing-effect 


T,(7) 


the monochromatic continuous opacity « 


calculations by Labs (1951)*. Some uncertainty exists in the exact relation 
between 7, and 7, owing to the discrepancy between the values of « originally 
calculated and the values required to obtain the correct value of the solar effective 
temperature; Aller and Pierce remedied this discrepancy by multiplying the 
values of « by a constant factor of 1-35. ‘This uncertainty affects the calculated 
intensity of the continuum (cf. Fig. 2), but should have no serious influence on 
line strengths. ‘The constants for the Aller—Pierce model are given in ‘Table I1; 
the m values were calculated without allowing for damping (a4 = 0) 

de Jager (1952) uses the variations in the profiles of the hydrogen lines across 
the disk to deduce the temperature distribution as a function of the optical depth 
7, corresponding to A 5000 and gives the values of log P, log P, for an atmosphere 
with 20 per cent helium. The corresponding values of +, were deduced from 
a table given by Neven and de Jager (1951) and the values of log P, log P, for our 
standard composition were computed by the iteration procedure described 
by Aller and Pierce (1952). The nm values were calculated twice-—with and 
without damping—and the two sets of equivalent width calculations give an 
indication of the importance of this effect. ‘The results are given in Tables III, IV. 


® The writer is indebted to Dr L. H. Aller for suggesting this procedure 
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Tasie Il 


Constants in Aller-Pierce model 
n (without damping) 
0” log P log P, B, (7, 4A 
0°00 0°05 O'10 0°20 
176 752 x 10" 
o8o 3 Bi 112 16°40 13°11 6°85 
‘072 4°06 “148 16°32 13°09 6°83 
054 4°20 233 14°78 11°91 6°26 


020 44°39 493 9°52 9°86 
979 4°55 629 9°52 "70 
956 4°62 766 $°26 

923 4°71 ’ ‘978 6°80 5°50 


577 4B 307 5°06 4°13 
541 4 'BB 3°544 4°13 3°40 
525 4°91 . 3°732 3°56 2°92 
797 4°94 ‘O38 2°90 2°35 


731 4°97 "145 2°56 
754 5°Ol "97 426 2°02 
718 503 840 1°42 
644 5°04 5°265 0°96 


652 §°05 5°722 0°67 
605 5°06 3°08 6°493 0°37 


Tasie Iii 
(Constants in de Jager model (1): n without damping 
n (without damping) 
log P log P, B,(r,) AA 
0°00 O'1O O15) 0°25 
700 x 10" 
028 16°40 6°86 
254 15°40 6°54 
398 13°37 5°89 


2°455 9° 5°49 
664 12°39 5°56 
354 10°20 4°4¢ 
254 9°43 4°19 


570° 6°55 2°94 
‘O89 5°O7 2°38 
923 4°31 
097 3°43 


255 2°35 
554 2°21 
Sur 
955 


154 
632 
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In his table (Model VII), de Jager does not give the value of the surface temperature 
T, (corresponding to r,=0°00). Our computations have therefore been carried 
out with two assumptions as to the value of 7, : 7, = 4270" (corresponding to 
Ty=O°O1, T,=0°003) and 7,= 4000" (close to the value deduced by de Jager 
from his observations, but not finally adopted by him owing to the higher 
temperature indicated by observations of the continuum at the extreme limb). 


Taste IV 
Constants in de Jager model (2): n with damping 
A= 10"n/Neg oo 
AA AA AA AA MA 


0°00 0'0s 0°20 


16°55 13°42 
14" I1'g! 
12°7 10°47 


0g7 
255 
416 
554 32 16 


S31 5386 442 0'300 
Oss 696 0°264 
"154 604 242 
632 *306 207 164 


Michard (1953) deduces the temperature distribution as a function of 7, from 
the observational data on continuous limb-darkening at various wave-lengths, 
but corrects it in the outer layers to give the observed excitation temperatures 
of Fe1 and Tit. His value of 7, is 4270°. ‘The values of 7, corresponding to 
each 7, were again deduced from the table of Neven and de Jager (1951) and the 
values of log P, log P, were taken over directly from Michard’s table for the 
adopted chemical composition. ‘The constants corresponding to Michard’s 
model are shown in Table V, damping being taken into account 

Figs. 3, 4 show the absolute and relative limb-darkening curves in the 
continuum calculated from the three models. ‘The observational data given for 
comparison were deduced by combining the value /,(o, 1) for the intensity at the 
centre of the disk given by Minnaert (1953) with the limb-darkening coefficients 
measured by Pierce and Aller (1951). ‘The agreement is satisfactory, bearing in 
mind the uncertainties in the solar constant and in the precise values of + 
(particularly in the model of Aller and Pierce). 


4 
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Taste V 
Constants in Michard model 
(with damping) 
log P log P. B,(r,) AA 
O15 0°25 
700 

2-83 

44 


c 


179 
156 
1607 


174 


204 
199 


o'156 


ae 


<—— pw 


Fic. 3 Intensity of the continuum at A115 500 calculated from t three models Dotted lan 
Aller—Pierce. Full line; de Jager. Broken line Michard. Black circles : observed values 


< po 


Fic. 4.--Relative limb-darkening curves caiculated from the three models, taking 14 (0, 0°5) 
as umty. Dotted line Aller-Pierce. Full line: de Jager. Broken line Vichard. Black 


crcles observe d values. 
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lable VI gives the line profiles and equivalent widths calculated trom de Jager’s 
Model VII (with damping; 7,=4000°). ‘The results obtained from the othe: 
models are omitted for the sake of brevity. 

In order to compare the computed centre-limb variations in equivalent width 
with the observational results of Table I, the computed curves of growth wer: 
interpolated for values of » n chosen to make the computed and observed equivalent 
widths coincide, in general, at »=0o°5. ‘This interpolation was carried out in 
practice by taking weighted means, for each j«, of the equivalent widths calculated 
for two successive values of 7» n, the weighting factors being chosen to give the 


coincidence at »=0°s. The results are shown graphically in Fig. 5 


Taste VI 


Line Profiles in de Jager model (T4000 ) 
entry gives the value of (1 —7r) corresponding to the given 


(a) 7] (1/10)" 


0544 oor. | 0059 
0459 Osis os50 
0273 | 0°0307 0325 
0120 | 0°0132 0139 
1043 0045 0044 
0018 OO1s | © OOo1s 


10 | oo!1oc 


35 


o’b 


o0g20 | ror! 0°1126 1199 125 ‘ 1540 1664 


0754 0563 00962 1027 Os 1335 | 1400 
2°0479 "0526 0°0557 | 0620 | ] 0824 0939 
0216 0235 0°o255 | 0273 | p o" 04S4 0407 
0083 0088 | 0°0093 0096 | ovorr2 ol 
0036 | 0037 0°003%5 | 0036 ' »*( 0045 007 
oo19g 0022 o’0018 | 0020 | ‘ ; 00 0016 0016 





69 22°07 y 20° 25°62 ta 4! 
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Taste VI (cont.) 


{c) 4 (1/4)n 


1624 
1405 , 
okg7 ° 1072 | 1133 
0420 ) | 0°0497 0526 
0157 o1go o1g9 0213 


1904 1997 | 2106 0'2225 
1666 | 1752 | 1852 1993 
1210 1302 
oss6 | 0602 
pas 76 | oot 
0073 0073 0070 9009 


0044 


0°0036 0039 0034 0033 


g2 


“V 


"9-00 


2571 2955 31558 | o' 3082 
2555 | 2077 | 
1803 | 0'1883 | O°2194 2218 
| © 0947 1160 1243 

0365 | | O'O4I5 0453 
0137 O14! 
0072 0073 


2926 2569 


1§$7 


o782 < 8 2? Ogos 
782 


0318 03§2 
a , 
O1349 O14) O140 


0072 0072 0073 





6 00 
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Taste VI (cont.) 
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An inspection of Fig. 5 shows the following principal features ;:— 

(a) All models predict the observed variation of the fainter lines (W<40 mA) 
within the observational error, except close to the limb («4 <0°3), where the 
Aller—Pierce temperature distribution does not predict a strong enough increase. 
Since the nm values (apart from damping) are nearly the same in all three models, 
we conclude that the differences are due to the different gradients of the source 
function B,(7,) in the region r,<0°3. The reduction in intensity due to a given 
amount of selective absorption at a given yj is determined by the gradient of 
B,(7,) in the relevant optical depth range (cf. equation (21)); the gradient with uv 
is therefore related more directly with the curvature of B,(r,) than with the 
gradient itself. ‘The influence of damping is shown by a comparison of the two 
curves corresponding to de Jager’s model; near the centre of the disc, the effect 
is approximately equal to the influence of the transition from the Aller—Pierce 
model (without damping) to the Michard model (with damping), so that we may 
regard the difference between the two curves in this range as being predominantly 
due to the fact that damping was taken into account in the latter; near the limb, 
the difference between the line intensities calculated from these two models is 
considerably greater than the damping effect, so that here the lower temperatures 
in the neighbourhood of the surface envisaged by Michard and de Jager seem to be 
confirmed by the observations. 

(6) For a line of intermediate strength (W=4omA; Fig. §(ii)), the curve 
based on Michard’s model reproduces the trend of the observations with 
remarkable precision out toy. =0°3. The variation deduced from the Aller—Pierce 
model is not sufficiently sharp (though it would be improved if damping were 
taken into account), while the increase deduced from de Jager’s model is somewhat 
too rapid. In accordance with the effect of the curvature in the source function 
previously mentioned, the high gradient with » deduced from de Jager’s model 
may be associated with the “‘ knee "’ in his temperature curve at r, = 0-22 (cf. Fig. 8). 
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Fic. 5.--Observed and computed centre-limb variations in equivalent widths, Open circles 
observed values for (i) A15565°38; (4) mean of A15§593°71, A1ss66°B8; (ii) Atssg7°71; 
(ww) Atss4z'15; (v) mean of Aits531°Bo, A1gs34°31; (ot) Ars662-05 Dotted line 
Aller—Pierce (without damping) for (i) qn-o-oBon; (ii) neo2sin; (Hi) Heo77s", 
(iv) y= 2°02n. Broken line: Michard (with damping) for (1) n«o-o64n; (ii) neo'17bn ; 
(11) n= o°540n ; (iv) ne t-2on; (v) geaics7n. Light full line : de Jager (i) without dampin 
for y=0o'075n ; (ti) without damping, for 4~0'194n ; (iti) without damping, for 40° 585n 
(iv) with damping, for 4 1°35n (T4= 4270") ; (v) with damping, for 9 1-Bon (Ty 4270") ; 
(vi) with damping, for n= 12n (T,~4270°). Heavy full line: de Jager (with damping) for 
(i) nw ov075n ; (4) ne ov1g4gn; (ii) ne Oo'sBsn; (iv) g=1-j30n (T,—4000") ; (0) ne 1-67" 
(T,= 4000°) ; (vi) »=9°6n (Ty 4000"). 


36* 
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The smooth temperature curves obtained by Michard and Aller and Pierce 
seem to be in much better agreement with the observations in this equivalent-width 
range. 

(c) For somewhat stronger lines (W=80mA; cf. Fig. 5(iii)), de Jager’s 
model reproduces the observations very well (bearing in mind the considerable 
observational errors for A15537°71), while the gradients deduced from the other 
models are considerably too low. 

(d) For still stronger lines(W=100 mA; cf. Fig. 5 (iv), (v)), none of the models 
considered seems to give as much strengthening towards the limb as is observed, 
although the representation by de Jager’s model (especially with 7, = 4000") 
is considerably superior to the others. ‘The agreement between de Jager’s model 
and the observations for A 15542°15 would be improved by the adoption of a still 
lower surface temperature, e.g. 7, = 3800" as suggested by Béhm-—Vitense (1954), 
but further computations show that one result of this would be an increase in the 
predicted intensities of the lines A15531°80, A15534°31 from p=0-3 to o-2 
(cf. Fig. §(v)); such an increase would be inconsistent with the observed behaviour 
of these two lines, so that it is no longer possible to improve the overall agreement 
merely by changing 7,. We prefer to use the observational material for the 
inductive derivation of a new model (Section 5). 

(e) The two strongest lines observed (200 mA <W <350 mA) both show almost 
precisely the same variation with yu (cf. Fig. 2; A15662-05 is shown in Fig. 5 (vi)), 
which is tolerably well reproduced by de Jager’s model. ‘The sudden drop from 
ji = 0°3 tO 4 = 0-2 is a Somewhat awkward feature to explain, but there can be little 
doubt of its reality since all four lines with equivalent widths in excess of 120 mA 

how the same feature (Fig. 2). Its significance will be discussed below. 

We conclude that the present series of equivalent width observations may be 
adequately represented, for the most part, by a temperature distribution corres- 
ponding to that of Michard for 7, >0-2 and to that of de Jager’s Model VII for 

,0°2, Certain adjustments will be required very close to the surface (7, <0°05) 
to account for the behaviour close to the limb of the strongest lines observed. 

\ study of the line profiles given in Table VI is of some interest in connection 
with the suggestion by Allen (1949) that the observed increase from centre to 
limb in the widths of faint lines may be due to some form of anisotropic turbulence. 
Doubt has been cast on this suggestion by Houtgast (1952), who finds that Allen’s 
observations of line broadening are consistent with isotropic turbulence if the 
detailed mechanism of line formation is taken into account. In Table VI, 
values are given of W/(1—rgo9), where W is the equivalent width in mA and 
Yo.o9 the central residual intensity. For lines with » =n in Table VI, 7,4, changes 
very little from p= 1-0 to yp =0°3, so that these lines are comparable with those 
considered by Allen; according to Table VI, the increase in line width is about 
10 per cent, compared with 30 per cent as found by Allen in the neighbourhood 
of 44300. ‘To compare the behaviour of the line wings in the two wave-length 
regions, we use two approximate relations given by Unsdéld (1948): 


1,(0, ») = Br, =) 
1(0,1)  B,{7,=1) 
din B, 
Riu) (Fe efesen (21) 
Ra) (=) Br(t,=1)’ 
dr, Jry~) 





(20) 
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where the function Ar(7,) represents the small additional optical depth present 
in a faint Fraunhofer line as compared with the neighbouring continuum. If we 
ignore changes in the second factor on the right of equation (21), it is clear that 
R(«), R(1) for the faint line wings in a given wave-length region is related to the 
degree of continuous limb-darkening at that wave-length. An inspection of the 
limb-darkening curves given by Minnaert (1953) for various wave-lengths shows 
that the gradient din B,/dr, increases rapidly towards short wave-lengths, as is 
to be expected from the increasing sensitivity to temperature of the monochro- 
matic Planck emission function ; this fact, incidentally, offers a simple explanation 
of the general faintness of infra-red Fraunhofer lines. The relative change in the 
intensity of the continuum from j= 1-0 to «= 0°3 is three times as great in Allen's 
wave-length region as at A15 500; the line broadening shown in Table VI may 
therefore be taken as a confirmation of Houtgast’s conclusion that Allen's data 
may be satisfactorily accounted for without the assumption of anisotropic 
turbulence. 

5. Inductive derivation of a new model.—¥For lines formed in local 
thermodynamic equilibrium it is possible to derive the distribution of temperature 
with r, from limb-darkening curves observed for any given portion of the line 
profile, provided that the relation between 7, and the relevant selective optical 
depth 7, is known. This relation is 


114) = jou + n(t,)}dt,, (22) 


so that 7, can be found from 7, if (7,) is known. ‘To derive our first model 
(‘‘ Model I’’), we assume the form of the line profile to be given by the values of 
W (1 —fo49) shown in Table V1; for a given yu, this ratio is nearly independent 
of the assumed model and of 4. ‘The observed equivalent widths were 
converted into values of ry.,5, and a limb-darkening curve was deduced for each 
line at the point of the profile corresponding to AA =o-10 A, using the values of 
I ,(0, x) derived in Section 4. The results are illustrated in Fig. 6. 
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Fic. 6.—Semi-empirical limb-darkening curves corresponding to the point 4A=o°to 
of each line profile 
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To derive the source functions B,(7,) corresponding to the limb-darkening 
curves, it is convenient to express B,(7,) either as a polynomial in 7, or as a series 
of exponential integral functions; each expression leads to a corresponding 
series for I(o,) and the coefficients can be found by least squares. In the 
continuous spectrum, the second type of expansion is greatly preferable (Pierce 
and Aller, 1951; Peyturaux, 1952), three terms being generally sufficient to 
represent the observational data well within their r.m.s. error. For the present 
problem, however, we are concerned with the source function very close to the 
surface; since the temperature distribution in this region is determined by the 
blanketing effect rather than by the normal conditions of radiative equilibrium 
in a nearly grey body, the physical significance of the expansion of the source 
function in exponential integrals disappears and what is required is purely an 
interpolation formula valid in the restricted ranges 


o2 €7,;< 1°, 


The data were expanded in series corresponding to both the polynomial and the 
exponential integral type, the coefficients being obtained by least squares. ‘The 
various values of ~ were assigned equal weight and in each case the expansion 
was adopted which gave the smallest sum of squares of residuals. The resulting 
source functions are given in Table VII, together with the type of expansion 
adopted for each line. The mean sum of squares of residuals is 600 (the unit 
being the third decimal place in 10°"/(0,«)), compared with a mean value of 
1000 from the r.m.s. errors in equivalent width. In Table VII, E.1. refers to the 
three-term exponential integral expansion. 


Tasie VII 


Source function distribution with selective optical depth 7, with 10" erg/cm*/sec 
as umt 


Continuum 15565°38 15593°71 155606°88 15537°71 15542°15 15531°8 
E.I E.I E.I E.I E.I. quadratic cubic cubs 
n=O Hoon neo Bsnyn-o2zinyn=o'shsn neign ne167n Hg bn 
(2°058) (1°379) (1°712) (1°645) (1°472) (1°643) (1°749) (1660) 
(2°B21) (2°665) (2°484) (2°433) (2°113) (1°B44) = (1855) = (1 04) 
20 3°230 3'093 2°906 2861 2°475 037 966 1368 
jo 3°534 3401 3°21! 3168 2°744 ‘223 084 iss 


‘40 8=—- 43°75 639 446 405 2957 456 204 
5° § 3°949 829 636 595 134 ‘736 2329 
60 4°10! 3985 792 751 284 ‘Bo1 457 
70 4°228 "IIS 3°923 881 414 040 588 


So 8— 44335 4°226 4°034 3°99! 3°529 3180 721 

9° 4427 4°32! 4°13! 4085 3631 3°313 555 546 
4°5°7 4°4°3 4°215 4°167 3°72 3612 990 2083 
4°34) (4°330) (3°924) (3°864) (3°325) (2°361) 


oo 
25 (4°667) (4°568) | 


To convert the values of 7, in Table VII into the corresponding values of r,, 
the assumption was made that the relevant values of n(7,, 0-10) and of » were 
the same as calculated for these Fraunhofer lines from de Jager’s model. +, was 
tabulated as a function of 7, for each line by using Simpson's rule to integrate 
equation (22) in the EDSAC. The resulting source function (Model I) is 


illustrated in Figs. 7, 8. 
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Fig. 7.--Source function deduced in the range o'00< 14° 0°15 Dotted line: Model I 
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Fic. 8. Source function deduced in the range 0-1 << 74< 10. Heavy line: Michard and 
suggested model, Light line : de Faver 
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To test the accuracy of Model I, the equivalent widths were calculated as 
a function of 4 by the method of Section 4. “The comparison is shown in Fig. g. 
The agreement, especially for the fainter lines, is rather poor, which is not 
surprising in view of the approximations used in deriving the model (particularly 








FIG. 9.--Observed and computed centre-limb variation in equvalent widths, Open circles : 
observed values for (i) mean of A15593°71, A 1556688; (4) Aigs37°71 ; (i) A tgg4garts ; 
(te) mean of A 1553180, A15534°31 ; (©) A1s662°05. Dotted line : Model I for (i) »~0-178n ; 
() n~0°-4Q90Nn ; (Hi) ye 1-009n ; (tv) ge 1°968n ; (0) p= 708m. Heavy full line: Model If 
for (4) no 1Bgn ; (i) qeosagn; (i) qe iaggn ; (fe) neiss6n; (v) n= 7-76m. Laight 
full line : Model III for (i) qeortgyn ; (i) neos7sn; (Hi) ge 1-g8on ; (ie) qei-738n ; 
(v) neg tan 
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the assumption concerning the values of »/n). For r+, 0-1, Model I coincides 
very closely with the temperature distribution of Michard, and the results of 
Section 4 indicate that the temperature distribution found in this depth range is 
probably correct. In the outer layers (7, <0°1), the temperature distribution 
was modified in an attempt to remove the outstanding discrepancies; two new 
models—Models II and II 1-——were derived in this way (Fig. 7) and the resulting 
equivalent widths were calculated by a simplified procedure based on the residual 
intensity at one point of the line profile (AA =o-os A) and on the line-width data 
obtained from the earlier computations. The results derived from Models II, 
III are shown in Fig. 9; the equivalent widths derived from Model II show 
satisfactory agreement with the observational data (within the r.m.s. error), 
and this model has therefore been adopted. Model II (with log P, log P, 
calculated for our standard chemical composition) is given in Table VIII, and the 
corresponding continuum intensities and equivalent widths are given in Table IX 
Taste VIII 
Suggested photospheric model (Model II) calculated assuming A~10*; no helium 
T 18600 T 5000 © log P log P, 
0'000 ‘000 196 
0-005 227 
o’o10 ‘ 180 
o’olrs 077 
054 
028 


002 
O89 


976 
963 
95° 
925 


giz 


574 


57 


762 5 1°93 

747 5 2°04 

1°7 734 5 2°14 

2 727 5 2°18 
6. Conclusion.—'Vhe calculations of Sections 4, 5 have shown that the present 
series of observational data is consistent with the photospheric temperature 
distribution in the deeper layers (7, 20°15; tT, 20°30) derived from the observa- 
tional data on continuous limb-darkening by Michard; the temperature 
distribution of Aller and Pierce is in close agreement with that of Michard in these 

layers if we assume the relation 
=O 77 
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given by Aller and Pierce (1952). Closer to the surface (r,<0°15; 7 <0°30), 

the situation is more complicated. Our Model I begins to depart markedly 

from the temperature distributions of de Jager and Michard (Fig. 7) for 7, <0-07 

(7,<0°16), while Model II continues to coincide with Michard’s distribution 

down to 7,= 0°02 (t,= 0°06). In the range 0°02 <7r,<0°12, Model II and the 
Tasie 1X 


Intensity of continuum and equivalent widths calculated from Model II 


Perro p¢ 08 pw@o6 p05 fPrO4 PK PKO2® f=Oo'l 
(calculated) 4°084 3°942 9 3°755 3636 3°492 3°312 3071 2°699 
Continuum 
(observed) 4°080 3'930 3°737 5 7° 3255 3°O55 
9 31°85 34 35°3 
(1/4)n 383 41°4 45°53 53 


(1/3)n 21°00 23°2 25 
56 
5s" 57° 


/ 


(1/2n 64°7 69°3 74°0 7 2 ; 
n 99°6 104°6 109°7 115 117°! 119°6 
2n 140°2 144°3 148 153 154°3 154°5 
4n 11's 185°2 188-2 1914 189°8 182°0 
models of Michard and de Jager are all closely coincident, so that we may also 
regard the temperature distribution in this region as well established. More 
ambiguity arises in the layers for which 7r,<0°02 (r,<0°06), which make a 
negligible contribution to the intensity of the continuum away from the extreme 
limb, but are strongly weighted when a Fraunhofer line due to a neutral metal 
is observed. ‘The temperature inversion obtained by us cannot be dismissed 
as the result of accidental observational (or truncation) errors, since the sudden 
drop in equivalent width of the four strongest 'ines observed from p =0-3 to 
jp. =0°'2 seems to point rather directly to the presence of such a feature although 
its position at 7, 0°005 may be wrong by a considerable margin owing to errors 
in the adopted values of » near the surface. Such errors are to be expected from 
the difference between the local kinetic temperature near the surface and the 
effective temperature 7, of radiation incident from deeper solar layers (Woolley 
and Stibbs, 1953; Béhm, 1954) which results in an enhancement of excitation 
and ionization in comparison with local thermodynamic equilibrium at the surface 
temperature 7). In a photospheric model with a steep temperature gradient, 
somewhat steeper than in de Jager’s model and in our Model II, Béhm finds an 
enhancement in the ionization of Fe by a factor of 3 for r,<o-o1. ‘The effect 
on the occupation of the highly excited levels of Fe 1, which are relevant in our 
case, has not been calculated in detail, but the undoubted presence of intense 
resonance absorption lines in the U.V. region (5°60 evA 2200) suggests that the 
radiation density at the frequencies relevant for excitation (as opposed to ionization) 
is likely to depart very little from that corresponding to local thermodynamic 
equilibrium. As far as the ionization is concerned, it would be of great interest 
to carry out similar calculations for all the elements, since Béhm’s factor 3 is 
strictly applicable only if there is no significant increase in the electron pressure ; 
if the overall number of ionizing processes in the outer layers of the photosphere 
were increased by a factor 3, the degree of ionization would be increased somewhat 
less than this. We conclude that departures from the Saha equation may require 
an increase by a factor of about 2 in our values of 7, close to the surface 
(tT, < 0°02, T, < 0°01) 
If our conclusion is correct, it would appear that the inversion in the 
temperature gradient that must occur somewhere in the neighbourhood of the 
hase of the chromosphere actually occurs within the photospheric layers 
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(7) = 0°015) and that the average temperature increases monotonically with height 
above this level. An interesting consequence of this is that we obtain the steep tem- 
perature gradient required to account for the increase in Fraunhofer line intensities 
towards the limb without being compelled to assume a particularly low limb 
temperature (7, = 4210 in Model Il). In this way, our temperature distribution 
may make it easier to account for the discrepancy noted by de Jager (1952) between 
his 7, and the value deduced from observations of the continuum close to the 
limb, but additional computations are needed to test this possibility ; furthermore, 
these data are subject to considerable observational errors. 

A further consequence of the new temperature distribution would be the 
concentration of molecules such as CO in the layer corresponding to the minimum 
temperature (t,=0°015; 7T'=4100°); the abundances of molecules would be 
enhanced by the presence of a substantial gas pressure at this level. However, 
molecular absorption lines would be somewhat masked by the substantial amount 
of continuous absorption present (especially in the far infra-red; this prediction 
seems to be consistent with the observed characteristics of the fundamenta! 
vibrational band of CO at A46000 given by Goldberg and Miller, 1953 It 
would be of interest to compare the centre-limb variations in the equivalent 
widths of CO lines with predictions based on Model II 

Finally, if Model II is correct, it must reproduce the observed curves of growth 
and excitation temperatures of metallic lines \ strict test from this point of 
view would require a series of elaborate calculations, but the excitation 
temperatures—at least at the centre of the disk—-cannot differ substantially from 
those of Michard, who expressly corrected his model in the outer photospheric 
layers to obtain agreement with observation in this respect. It is therefore 
reasonable to conclude that Model II will give good agreement with observations 
of line intensities at the centre of the solar disk, although a complete discussion 
will have to include the effects of temperature heterogeneities. 
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FINITE ATMOSPHERES WITH ISOTROPIC SCATTERING 
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Summary 


Some transfer problems in plane parallel atmospheres of finite optical 
thickness 7, are considered, the scattering being isotropic with albedo mw 
(o<m<1). The equations for Chandrasekhar’s X- and Y-functions are 
obtained by a simple mathematical method, an extension of one originally 
due to Ambartsumian. The same method is then used to obtain expressions 
for the emergent intensities at both faces of the atmosphere when there is 
thermal emission and when it is subject to incident radiation 





1. Introduction.—Problems concerning the transfer of radiation through 
plane parallel finite atmospheres are more complicated than the corresponding 
ones for semi-infinite atmospheres because there are two faces on which incident 
radiation may fall and through which radiation will emerge. Hitherto the only 
successful method of tackling such problems has been by means of Chandra 
sekhar’s Principles of Invariance and it is significant of the difficulties involved, 
that the only people who have worked on such problems are Chandrasekha: 
himself and a few astrophysicists associated with him at Yerkes. 

In two recent papers ((3) and (4)) | have developed a method which was 
originally given by Ambartsumian in (1) and of which there is an account in (9), 
Section 28. ‘This method has proved to be unexpectedly powerful and in this 
paper | show that it can be applied to finite atmospheres as well as to semi- 
infinite ones. Only isotropic scattering with albedo w(o<@ <1) is considered, 
but it should not be difficult to adapt the method for other laws of scattering. 
(See, for example, (4) where the scattering is non-coherent. I have also used 
the method for interlocked multiplet lines.) 

Ambartsumian’s method only involves simple mathematics. It starts from 
the form of Milne’s first integral equation (Section 2) appropriate to the problem 
and it depends upon the solution of an auxiliary integral equation (Section 4), 
which leads directly to the X- and Y-functions, Solving the auxiliary equation 
is equivalent, in Chandrasekhar’s work, to finding expressions for the reflected 
and transmitted intensities in the “ diffuse radiation field” when one surface of 
the atmosphere is illuminated by axially symmetric radiation inclined at a fixed 
angle to the inwards drawn normal (see Section 5). 

The solution found for a given problem is unique only if the corresponding 
homogeneous form of Milne’s integral equation has no non-null solution of the 
appropriate order of magnitude. ‘The solutions of the homogeneous equation 
are investigated in detail in the appendix and the results are summarized in 
Section 6. 

When o<w<1, the X- and Y-functions exist and are unambiguous, but 
when w=1 (the conservative case), if X() and Y() are solutions of the X- and 
Y-equations then so are 


X(u) + Op[X(u)+ Y(w)} and = Y(pn)— Op[X(yu) + Y(p)] 


for any value of the constant 0. This ambiguity is related to the question of 
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uniqueness and it is also discussed in the appendix. ‘Throughout the body of 
the paper, X(u) and Y() will represent the standard solutions when w= 1, i.e 
those solutions which are such that,* 


“i “1 
xy | X(u)du = 1, yom | Y¥(u)du =o. (1.1) 
0 0 


In every problem it is necessary to discuss the conservative case separately. 

With the aid of the auxiliary equation, a number of problems can be solved 
in terms of the X- and Y-functions. In Sections 7(0< m<1) and 8 (w=1) the 
emergent radiation from an emitting atmosphere is considered. A polynomial! 
representation is used for the Planck function and from this the results found by 
Horak and Lundquist in (8) for a linear Planck function are deduced. 

In Sections 9 (o<m<1) and 10 (w= 1), one surface of a non-emitting atmo- 
sphere is illuminated by isotropic radiation and the reflected and transmitted 
intensities are found. The results of these sections can also be used to give the 
intensities of the emergent radiation from an atmosphere emitting according to 
the law E,{7) (the second exponential integral), where + is the optical depth below 
the upper surface of the atmosphere. By combining the results of Sections 


7-10, it is possible to deal with a representation such as 

a+br+cE{r) 
for the Planck function, a form which has been found to give a good represen- 
tation in the grey case. 

‘Tables of the X- and Y-functions for values of 7, (the optical thickness of the 
atmosphere) and ow such that o<r,<1, 0-5 <m<1t have been compiled by 
Chandrasekhar, Elbert and Franklin (6). It is probable that tables for larger 
values of 7, will be needed before much useful numerical work can be done on 
finite atmospheres 

In this paper details of rigour are omitted. In (3) | have treated the same 
problems rigorously for a semi-infinite atmosphere and in general the treatment 
for a finite atmosphere is similar. ‘The existence of the X- and Y-functions, at 
least for sufficiently large 7r,, has been established in (2) and I shall therefore treat 
them as known.functions. 

2. Milne’s first integral equation..-Wlien there is isotropic scattering, the 
equation of transfer is of the form 


pa U(r, ») I(r, 1) ) or | Ae de Br), (2.1 


where B,(r) may be (1 —) times the Planck function or it may be some other 


d 


given function, 

Following the notation of Chandrasekhar in (§), I shall use /(7, +) (0 <p <1) 
to denote the intensity of the radiation directed towards the surface r=o0 and 
I(r,~p) (O<p <1) for that directed towards the surface r=7,. Then if 
radiation of intensity /,(j) is falling on the surface r =o at an angle cos~'y« with 
the inwards normal, and if radiation of intensity /,(j) is similarly falling on the 
surface r= 7,, the boundary conditions are 


K(o,—p)=fo(u) (O<p <1), 


(2 
I(r, +p)=f(p) (o<p <1). 


* All solutions satisfy the equation x, + 4 
equations (1.1) therefore determines O 


owt in the conservative case. Either of the 
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If we write (2.1) mm the form 


d 
oP I(r, wp) = I(r, 4) — Wr), 


“1 
Ur)=ia I(r, u')dp' + Br), 
a 


then the solution of (2.3) subject to the boundary conditions (2.2) is given by (cf 
(5), Chap. I, equations (64) and (65)) 


" dt 
I(r, +) =f(u) exp[—(7,—17) ue] +exp(r.) | O(thexp(—t/u)—, (2.5) 
fe 


r 


(2.6) 


, dt 
I(r, —p) =fo(u) exp (—7/4) + exp( 7) | Athexp(ty)—, 
0 i 


where 0<yp <1, 0 <7 <7, in each equation. Hence 


“1 “1 l 
I(r, ) dy | U(r, +p)dp 4 I(r, — x) dd 
l 0 0 


1 1 
ffpwjexp(—r p)du+ fi (ujexp | —(r, 
0 


+ | B(VE Y(t —7]) ae, (2.7) 


tee d 
Ey (r) exp (—rTx) > 
J} 


(2.5) 
x 

On substituting from (2.7) into (2.4), we have the following integral equation for 
the source function &(r): 


r"% i 
Ur)=ho | BAHL (\t r|) dt + ho] fle')exp( —r/u') du 


i 
t ber | f(y jexp[—(t,—7) we’ | du’ + Br) 
0 


This is the form taken by Milne’s first integral equation. 
We shall introduce the operator A defined by 


Lif(op=4 | SE Mt ~7| dt. (2.10) 
0 


In this, and in all subsequent equations, o<7r<r,. When 71, is infinite, A is 
identical with Hopf’s operator A. (See (7), and (g), Chap. II.) Equation (2.9) 
can now be written 
Br) = wA,{B(t)} + Br), (2.11) 
where 
1 
B(r) = Br) 4 ar | folie )exp ( rp’) dy’ 
J0 
1 
+ io f(y )exp| (tT, —T) pw’ | du’. 
Jo 


This is a known function in any particular problem. 
The emergent intensities from the atmosphere are /(o, +) and 
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By (2.5) and (2.6), these are given by 
Ko, +p)=f/lvexp(—ty 4) tide *), 2.13) 
Nr, h)=Sleyexp(—ty 4) +h’), (2.14) 


j,(s)=5 ; B(t) exp ( — st) dt, (2.15) 


/0 


j,(s) = s exp ( 78) | B(t)exp (st) dt | Or, thexp( —st) dt. (2.16) 


0 
In the following sections we shall aim at finding j,(u~") and j,(u~') when B(r) is 
some standard function. 

The functions j,(s) and j,(s) will be called the 7,-Laplace transforms, or 
nore briefly the 7,-transforms, of B(r). From (2.15) and (2.16) we have the 
important relations 

Io(S) exp (—s7,)hy(—5),  $4(5) exp (— ST, hol —5), (2.17) 
of which we shall make frequent use. We shall always use the corresponding 
lower case letter with subscripts o and 1 to denote the 7,-transforms. For 
example, the 7,-transforms of J(r) will be jo(s) and /,(s). 

Note that, if * 


Br) = aB'(r) + bB*(r) 
and if 

JI'(r) oA {J *(t)} Br), 

J*(r) o\_{J*(t)} ' B*(r), 
then a solution of (2.11) is 

U(r) =aJl'*(r) +b) %X(r), 
and jo( 5) = Gjq'(5) + bjg*%(s), 

1,(5) aj ,*(s) + b7,%(s). 
hus the solutions corresponding to standard functions can be added to give 
solutions (or approximations to the solutions) of more complex problems. 


3. Properties of the operator A In deducing (2.21) from (2.19) and (2.20), we 
have used the fact that the operator A is linear, i.e. if a and b are constants, then 


\, {af(t) + be(t)} = aX, { f(t)} + bA{g(t)}. 
in addition to this \ has the following important properties : 


I 1, f())} =A, f{ f(r, — 9}, 


d d ’ . 
II gy ttt) 0,4 Sf) + 4 OEW(7)—BF (rs)E rs Tr), 


the function f(t) being subject to appropriate conditions. (Cf. (7), p.34.) 
To prove I, we have, by the definition (2.10), 


\,,~ A S(t)} = | "NE (r, - 7—t)dt+ sf" I(DEKt—7, +7) dt. 


0 sd eed 





* Because of the subscripts o and 1, it is necessary to write other subscripts as superscripts. 
As no powers of functions are involved, there is no risk of confusion. 
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On putting t=7,—u, we get 
X,,-fO}=3 [" fers —w)E (ur) du 4 {fry —u)By(r —u) du 
0 


- A f(r, ag u)}. 
To prove II, put t=7—w in the integral over (0, 7) and t= 7 + in the integral 
over (7, 7,) in the definition (2.10) of A,{f(t)}. Then 


NO} =b [fr upEy(uddu+ § [flr +u)By(u) de (3-4) 


By the usual formula for differentiation under the integral sign, 


, vd / 
FAAMO} = | Ff —u)B iu) du + M0) (7) 


2. ind 


‘ if. £ f(r + u)E,(u)du— f(r, Er, —7), 


and this (cf. (3-4)) is (3.3). 
4. The auxiliary equationand the X- and Y-functions.—The auxiliary equation is 
J(r, 0) =wA,{J(t, c)} + exp(—or) (o<m <1), (4.1) 
Just as in (3), Theorem III, it can be shown that, when o<@ <1, (4.1) hasa 
solution which is bounded for o >o and o <7 <7,, the bounds being independent 
of 7, We shall use J(7, 0) to denote this solution. When w«1, J(r, @) is 
ambiguous in the same way that X(y) and Y(yu) are ambiguous. This difficulty 
will be dealt with in the appendix. For the present it is enough to say that J(r, a) 
is a solution of (4.1). 
Differentiate (4.1) with respect to r, using (3.3). ‘Then 
fa] é } ‘ 
a I(r, a)=a\, 15H vr + hard (0, 0)E,(7) 
bad (7,,0)E,(7, —7) — a exp(—aor). 
Hence 0 J 0 
5, (7,9) +ad(r,a)=a\, 29) ' ad (t,o) 


+ jer J(0, 0)E,(r)—}arI (14, 0)E\(7, ~7). (4.2) 

We next obtain expressions for £,(r) and E,(r,—7) from (4.1). Multiplying 

(4.1) by ao! and integrating with respect to o over (1, ©), we have (on writing x 
for a)* 


G thew dx) F 
ee —} +E,(r). 
w\, , I(t, x) a| E(t) (4-3) 


On replacing 7 by 7,—7 and using (3.2), we get 


fa) d f d : 
J (1, 7,3) = -oK,4 | J (1, nS} + E,(r,~—7). (4.4) 


From (4.2) subtract }wJ(o, 0) times (4.3) and add |wJ(r,, 7) times (4.4). 
Since the operator A is linear, we see that 


¢ Ir,0)+0I(1,0) }wJ(0, 0) | "Jee 
CT i x 


‘@ d 
+ jer J(7,,0) | I(r, —7,*) - = CJ(r), (4.5) 
1 





* An integration with respect to « commutes with A, which is independent of o. In general 
we shall invert orders of integration without comment, 


37 
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where J(r) satisfies the equation 
I(r) =wK,{J(t)} (4.6) 


and (' is a constant which may depend on a. 


The solutions of (4.6) are discussed in the appendix. As we shall show that 
C= 0, the value of J(r) does not matter. 


Let R(s,a) = (” exp(—st)J(t,a) dt. (4.7) 


70 


Then we have the important relations 
R(s,0) = R(a, 5), (4.8) 
exp (—s7r,) R( —s, 0) = exp(—or,) R( —a, 5). (4-9) 


To prove (4.8), multiply (4.1) by J(7,s) and integrate with respect to r over 
(o,7,). Then 


f J(7,8)J (7,0) dr = ber \° J(r,s) dr | " I(t, o)E(|t—r| dt + R(o, 8) 
0 0 0 


aie {° J(t, 0) dt [" J(r,s)E(\7 —t|)dr + RCo, 5) 
0 0 


= ° J (t, a) (t, s)dt ~ R(s,0) + R(a,s) 


(by interchange of s and o, ¢ and r in the first line) and hence we have (4.8). 
To prove (4.9), write r,—7 for r in (4.1). Then 


JI(1,—7,¢) wi, iJ (t,0)} + exp [—o(r,—7)]. 
Hence 


: J (1,8) (1,—7,0) dr = jw | J(r,s)dr % I(t, o)E(|t +7 —1,]) dt 
0 0 0 
+ exp(—or,) I. exp (or)/(r, 5) dr 


=}o f° J(t,0) dt [ I(r, s)E (|r +t—1,]) dr 
+ exp (—or,)R(—a,5) 
= " J(t,a)J (7, —t, 5) dt — exp ( —#7,)R(—5,¢) 


+ exp (—or,)R(—<o,5). 
But (writing t= 7r,—r) 


| " I(t, 0)J(r, —t,s)dt = I I(r,—1, 0)J (1,8) dr, 


and so (4.9) follows. 
Now multiply (4.5) by exp(—sr) and integrate with respect to r over (0, r,). 


Since 


f" J(r,—1, x) exp (—sr) dr= [ J(t, x)exp [—s(r,—2)] dt 


= exp(—sr,) R(—5,2), 
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we have (after integrating the first term by parts) 
exp ( —s7,)/(7,, 0) —J(0, 0) + sR(s, 0) + oR(s, 0) — bo (0, 2) {- R(s, x) = 
1 
+ mJ (7,,0) exp (—s7,) . R(- sx) =Cs¥,(s), 
1 


where j,(s) is the first r,-transform of J(r). (See equation (2.15). ‘This may be 
written 


(s+o)R(s,0)— Cs 4j,(s) =J(0, a) {1 +) i R(s, x) < } 


° d. 
—J (7,0) exp ( or) + jor | R( - 6) Sf 


Put r=0 and write s for o in (4.1). Then we have 


J(o, 8) = he \° I(t, s)E(t) dt +1 


r 


=I+ jer | ' I(t, s) dt [ exp( - xt) 
J1 ; 


= 1+ ae R(x, s) S (by (4.7)) 


=I+ jor | R(s, 2) = 
, : 


by (4.8). 


Again, putting r=7, in (4.1), we have 


J(1,,5)= iw |" J(t, s)E\(r, —t) dt + exp ( —sr,) 
= exp(—s7,)+ jor” I(t, s) dt # exp [—(r,—#)s] < 


= exp (—sr,)+ jor | exp ( —7,x) R( xo 
1 


= exp ( ~sr,)41 +ho| R( 5,2) SI 


/ 


by (4.9). Hence (4.11), (4.12) and (4.13) give 
(s +o)R(s,0)— Cs“ ¥o(s) = J (0, 0) (0, 8) —J(1,, 0S (7,, 5). 


Since, however, R(s,c) = R(o,s), we must have C = Ao~'j,(0), where A may be zero 
Hence 


(s+0)R(s,0)— Asa )(s)j(a) = J (0, 8) (0, 0) —J(74,0)4 (7, 5). (4.14) 
From (4.12) and (4.13) we have the important relations 
J(t,, —5) = exp (sr,)J (0, 5), 
J(o, ~s)= exp (sr,)J(7,, 5). f 
Write —s for s in (4.14) and use (4.15) and (2.17). 
Then on multiplying by exp(—sr,), we have 
(s—o) exp (—s7,)R( —5, 0) + As-'o~}j,(s)jo(a) = J (0, s)J (71, 0) — J (0, 0) I (4,8). 
(4.16) 


(4.15) 


37° 
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Since exp( —s7,) R( —s,0) =exp(—o7,)R(—a,5), either A=o or 
lsd) = —jloYols)- 


In this ‘we can put o =s and we get j,(s)=0 or j,(s) <0 and either implies that /(r) =o. 
Hence (4.14) and (4.16) become 


(s+0)R(s,0)=J(0, s)J (0,0) —J(1,, 58)4 (7, 0), (4.17) 
(s a) exp ( 5:3 57) R( ~$,@) =J(0,s)J(7,,0) ~ J (0,0) I (7, 5). (4.18) 


From (4.12) and (4.17), (4.13) and (4.18), we have the following integral 
equations for J(o,s) and J(7,, 5): 


J(0,s)= 1 + an {J (0, s)J(0, x) —J (7, 8)J (7, *)} qa" (4.19) 


(@ dx 
I(r, xp (— j {J(o,s)J (7, J(o, x) J (71, 8)} ——. - : 
(74,8) = exp(—s7,)+ jor | (0,5 (ry *)-JO, Wr» Yay (4 
Now put s=1/u, x=1/u’ and write 
X(p)=J(0,u™"), ¥(u)=S (71,4 *). 
Then (4.19) and (4.20) become the X- and Y-equations: 


. a - : on tee 
X(u)= 1+ boom | (XW) ¥(u)¥(u’)} ——, 


perp 

, wee on 

Y (m) = exp( — 74/4) + omy | {¥(w)X(p') —X(m) Y¥(u')} —— . 
These are (5), Chap. VIII, equations (1) and (2) with ‘Y(y)= $a. 

The existence of X(u) and Y(u) when o<aw <1 follows from (4.21) and the 
existence of J(r, a). 

In the appendix it is shown that, when m = 1, the function J(7, 0) can be chosen 
so that the functions defined by (4.21) are the standard solutions, i.e. those 
satisfying (1.1). We shall assume that this has been done. 

From (4.17), (4.18) and (4.21), we have 


R(u*, ug?) = SE® {X(u)X(49) — ¥(4) ¥(H40)} 5 
b> Ho 
exp (—7,/W)R(~ 0's wy!) = re (X(4)¥ (40) — ¥(4)X(u9)}- 
Lo L 
Note that the relation 
Y (pm) exp ( Tp) X( ps) 


comes from (4.15). 


5. Diffuse reflection and transmission.—l\f the surface 7 =o is illuminated by 
axially symmetric radiation at an angle cos’ '4,(0 <j, < 1) with the inwards drawn 
normal and if F is the intrinsic flux of the incident radiation, then (see (9), 
section 31.2) with the notation of Section 2 


fle) = } FO(u — Ho), (5-1) 
where 4 is Dirac’s 5-function. By the properties of the 5-function 
A 
| folu)exp(—7/u) du = 4 F exp(—7/1H9)- 
“9% 


If f,(«)=0 and B,(r) = 0, (2.11) and (2.12) give 
B(r) = aA, {B(t)} + Jo F exp( —r/j19). 
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Comparing this with (4.1) we see that a solution of (5.2) is 
B(r) = Jor FI(r, p49"). 
From (2.13), (2.15) and (5.3), 


I(o, +) = bork | "S(t, Ho ‘)exp(—t/m) : 
0 


7 jo Fu *R( 1 ie 1) 
by (4.7). Hence (4.24) gives 
— } ‘ . , 
I(0, + 4) = fo P —*— [X(u)X (Ho) — ¥(H) ¥(H0)}- 
+ Ho 

From (2.14), (2.16), (5.1) and (5.3), 

I(r), —#) = § FO(u — no) exp (— 71/49) + exp (—7,/4) Fae *R(— 4, wg). 
The first term on the right gives the “ reduced incident radiation” and the second 
term, by (4.25), gives the following expression for the transmitted intensity in 
the ‘‘ diffuse radiation field”’ : 

P V ‘iy F 4 ‘ i - « 
| Fo —+— (X (49) ¥ (x) Y(H9)X(4)}. (5.5) 
=e 
The expressions (5.4) and (5.5) are those given in (5), Chap. VIII, equations (108) 
and (109). 
6. The homogeneous Milne equation.—\t is important to determine whether 
the equation 
I(r) =a, {J(t)} (6.1) 
has any solutions other than J(r) =o, for if 3\(7) and 3*(r) are solutions of the non- 
homogeneous equation (2.11), then 3'(r) — 3%(r) satisfies (6.1) and 
B47) = Br) + C/(r), (6.2) 
where C is an arbitrary constant. 

The investigation of (6.1) is rather long and it depends on some relations 
between the X- and Y-functions which I have proved in (2). I give below the 
main results, leaving all details to the appendix. 

As in the semi-infinite case, we have to consider the characteristic equation 

wm 1+€ 

5c In —a =I, (6.3) 
If o<m <1, this has two roots C= +k (o<k<1). Whenw«= 1 it has the single 
root C=o. The values of k corresponding to different values of w are given in 
(5), Table I. 

When o<m~< 1, the functions X(z) and Y(z) of the complex variable z have 

simple poles at the points z= +1/k. Let 
. X(z) 
L = lim 37 9 (6.4) 
zoi/k y (2) 
so that L is the quotient of the residues of X(z) and Y(z) at z=1/k. Then the 
7,-transforms corresponding to the non-null solution J(r) of (6.1) are given by 
B 
——— {LX(yn)— Y(h)}, 6.5 
14 hey ( ) () yy ( 5) 


B : 
k {X()- LY()}, (6.6) 
fA 


A 
l otpnie, aD r $4 
Jolt ) 1 hy X( ) LY (y)} 


A 
——— J "(yu (u)' 4 
"Tok [LX (un) — ¥(u)} 


where A and B are arbitrary constants 
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For all A and B other than some special values such as A = B=o0, J(r) has a 
multiple of exp (kr) as its term of largest magnitude for large r. ‘Thus (6.5) and 
(6.6) must be added to the emergent intensities found for any given problem 
unless it is physically impossible for the source function 8(r) to increase at expo- 
nential rate. ‘The constants A and B must be determined in some other way, 
e.g. by the flux through the surfaces 7 =0 and +r =7,. 

When w = 1, the 7,-transforms are given by 

ji *) = AX(n) - BY(p), (6.7) 
ju") = BX(u) — AY(y), (6.8) 
where A and Bare againarbitrary. Here X(u)and Y() are the standard solutions. 

In this case J(r) is only O(r) for large + and the addition of (6.7) and (6.8) to 
the emergent intensities found for a given problem is likely to be more important 
than in the non-conservative case. 

7. An emitting atmosphere (o<m—1).—In this section we shall take 
fl) =O, fie) =o and r” ) 

By(r) = a9 “G1! Teor TG, (7.1) 
in (2.12). Because we can deal with recurrence relations, it is a8 easy to work 
with (7.1) as with B,(r)=agr+a,. In the latter case, we have to deal with 
B,(1r) = ay first, and the work is nearly doubled. 

Let B(r)=J"(r) when B,(r) has the value (7.1). ‘Then the Milne equation 
(2.11) becomes 


” nl 
I(r) = oA {J"(t)} + a, — 


P 
+ @,>— 
'(n—1)! 
where n=0, 1, 2,.... For convenience define J~"(r) to be identically zero. 
From (7.2), by the differentiation formula (3.3), 


d d ; , 
ZI=wr,{ 5M} + bed "(O)E, (17) — bord "(7,)E\(7, — 7) 


the polynomial being absent when n=o0. Hence 


d mw ite ; ; 
= I"(r) J" Mn) =o, 4 5,I%0) J" XN) + bed "(0)E,(1) ~ food "(rE (ry ~ 7). 


Because J~'(r)m0, this holds for n=o. If we now subtract }aJ"(o) times (4.3) 
and add 4mJ"(r,) times (4.4), we get the homogeneous Milne equation 


d aad d: aad dx 
5 "(r)-I" Mr) - jonJ™(0) | I(r, x) = rs bos™(r,) I(r, —7,%) ~ 


1 
x {a ° d: + d 
wi Bd") —J* t) — jars (0) | Mt VS + hod") | It) =I, (7.3) 


and a solution of this is 


d iw dx 
z7 (r)—J "r)— ho/ (0) ; J(r,%) - 


+ bed "(7,) Ir-79)2 =0. (7.4) 
| 


This holds for n =o, 1, 2,... 
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Multiply (7.4) by exp(—sr) and integrate with respect to r over (0, 7,), in- 
tegrating the first term by parts. Then, with the usual notation for r,-transforms, 


i? d 
exp ( —s7,)J"(7;) -J™(0) +Jo"(s) -s Yo" (s) jon) ™(0) | R(s, x) = 


* d. 
+ }aJ/"(7,) exp ( 7)) R( sx) =O. (7.5) 
, : 


In this we have used (4.7) and (4.10). By (4.12) and (4.13), we have 

Jo"(s) —5 Yo" (8) =I "(0)J (0, 5) —J "(7,4 (7, 5). (7.6) 
On putting —s for s and using (2.17) and (4.15), we get 

"(sh +9," (s) =S"(7,)/ (0, 8) —J"(0)/ (74,8). (7-7) 
If we now write s=1/« and use (4.21), (7.6) and (7.7) become 

jo"(v ") Hq” ys ') =J"(0)X(p) J") ¥(p), (7.8) 

Ae") + ay” (ee *) =I (7) X(n) —F™(0) ¥ (pn). (7-9) 


Since J~\(r)=0, therefore j, (s)-0, j, (s)>o and these equations hold for 
t=0,1,2,... Solving, we get 


Jo"(u 1) = X(p)XS"(0)u” r V(p) LI", yo” . (7.10) 


r= re0 


Ne *) = X(w)XI (7, — pw)” * — Y(e) LI"(0)( pe)". (7.11) 
r=0 r=0 
We have now to find the 2m +2 constants/’(o), J’(r,)(r=0,1,...,"). Put 
7=0 and r=7, in (7.2). Then 


J") = hor (" J"(1)E,(t) dt +a, 
Jo 


pe dx 
, n 


=hw (” I™(t) dt | exp (— xt) 5 +a 
J0 J 


*@ 1 
=tor| joMx) S +a, (by (2.15)) 


rl 
= lw Jo" (4 ) du T@.. 
“0 


JI%1,) = hor % IMO)E (r,t) dt-+ La, > 


*, «o dx n rT,” ’ 
jor | J” - Pe oa eae 
be | (t) dt [ exp[ —(7, — ¢)x] : +a yl 
‘oO dx " 7," j 
= jer | Ans) +2a,~t (by (2.16) 
cl ” » &-0 
wi “(u")\du+ La,-+.. 
jor) ( )d Soo (n ry! 


‘1 
x= ber | u'X(i) du 
0 


rl 
y= jor | oY (u) dy 
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These definitions agree with those of (5), Chap. VIII, with V(u)=4n. By 
Theorem 2 of the same reference 

(1 —%) — yer =1 —w. (7-15) 
(When mw = 1, this leads to the relation x, + y, = 1 and hence to the conditions (1.1) 
for the standard solutions.) 


On substituting from (7.10) and (7.11) into (7.12) and (7.13), we obtain the 
equations 


J"(o)— vJ"(0)x, 4 vIn), 7 (7.16) 


r~0 r= 


J 74) — LI" 7 — 1)" x, + LI"(OM —1)""y,_-= X 
r=0 r=0 ' " ee 


a: (7-17) 


From these equations with n=o, 1, 2 J%o), J%r,), Jo), J4{74),.-. can be 
evaluated in turn. Then (by (2.13) and (2.14)) the corresponding emergent 
intensities are given by (7.10) and (7.11). 

When n= 0, equations (7.16) and (7.17) become 


IO 1 — x9) +I 1 )¥q = Gp, 
I%o)yg +I%7, (1 — Xp) = ap. 
Hence, by (7.15), 
JS%0) = I%7,) = ag( 1 — #9 — Vo) (1 —m). 
From (7.16) and (7.17), with n= 1, 
JO) 1 ~ X%) +S "7, )yg —F%(0)x, +I %(7, 1, =, 
J0)yg +I (7,1 — %) —F%(O)y, +S 7 ey = Gg7, + 4. 


From these equations, (7.19) and (7.15), we get 


(7.20) 


(1 —a)J (0) = a,(1 — %y — Vo) + Gol %, — Vy — TLV 0) } 


(1 —a)J(7,) = a,(1 ~ Xp — V9) ~ Gol *, — Vy — 7, (1 — %o)]- 


Hence (7.10) and (7.11), with n= 1, give for the emergent intensities correspon- 
ding to linear emission 


X(u) 


I(o, +) = ES (agp + Ay 1 —%y— Vo) + Gol *%, — ¥y —71¥0)} 


¥(u) 


XH) 
I 


= tol + Ay) 1 — X%q — Vo) — Gol, — Vi — T1(1 — %)]}, (7-21) 


I(r,, = mh Gols + aI — Xp — Vo) - a,[ x; ~F7i —7,(1 —%X»)|} 
¥ (4) 


: = i oft + 4,1 ~ Xy— Yq) + Gol % — Vy To}: (7.22) 


These values agree with those found by Horak and Lundquist in (8). ‘They have 
used the moments «, and £, instead of x, and y, and have not made use of (7.15). 

8. An emitting atmosphere (aw = 1).—When w= 1, the coefficients of J"(o) and 
J"(r,) in (7.16) and (7.17) are zero by (1.1). The equations are then obviously 
wrong. For example, the equations (7.18) can only be satisfied if a, =o. 
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In this case we have to use the general solution of the homogeneous Milne 
equation (7.3), *viz. 


a 
, x 


IM) IM) — 1d) | 
fe dx 
+i") | J(7,—1,%)— =S(r), (8.1) 
where, from (6.7) and (6.8), 
jhe )=A,X(u)-B,Y(u), | 
j(u)=B,X(u)—A,¥(u). J 
Then an additional term s~j,(s) has to be inserted on the right of (7.5), and 


Pid"), — J (e*) have to be added, respectively, on the right of (7.8) and (7.9) 
Equations (7.10) and (7.11) become, by (8.2), 


(5.2) 


Jo" (ee) = X(p) LI "(0)" — Y¥(u) I(r, ue” 


r«0 r=-0 


teX(p)d Ap’ —pY (pn) UX Bp", 


r=0 r#0 


Ae) =X(p) LI — py" — Y¥(u) LION — py” 


r=0 r=0 


n n 
eX(pyX BL py" +eY(w)L Alp 
re ra 
Write 
Co Ag, Do By, C, 4,+J° (0), D,=B J’ Yr) et) 
Then these equations become 


joe) =I*(O)X (pn) —IM( 7) V(H) + XW) UC wt V)EDerr, (8.3) 


rT=0 r= 


fie") =I (7 )X(u) -I"(0)V(u) + X(H) 2 Dp)" = VLCC —pyrr?. 
r=-0 ra 
(3.4) 


The equations (7.12) and (7.13) (with w =1) are unchanged. On substituting 
for jo" *) and j,"(~*) from (8.3) and (8.4) and using the conditions x, = 1, vg =0, 
the terms containing J"(o) and /"(r,) disappear and we are left with 

n 


n 
SC Sy-2417 © DV e041 re (8.5) 


r=0 r=0 
ad n vad Tr n 7 
LDA — 1)? +2, 43 —- DC l(—1)" tat La, +. (8.6) 


r=0 r=0 r=0 (n r)! 
These hold for n=0,1,2,... and the constants Cy, Dy, Cy, D,,... can be deter- 
mined successively. ‘Then the 7,-transforms of /"(r) are given by (8.3) and (8.4), 
J™o) and J"(r,) being arbitrary. 
When n =o, (8.5) and (8.6) become 


Coe, — Dey; dp, 


Coy, — Dox, dy. 





* Had this been done in Section 7, after elaborate reductions, we should have obtained the 
same emergent intensities plus those given by (6.5) and (6.6). 
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Hence 
Ca= — Dy= — o/(*; + 9). 
However, by (5), Chapter VIII, Theorem 8, with ‘¥() =}, 
x? —y?=4, 
and so (8.7) can be written 
Uo Do = ~ 34(f%, ~ ¥1)- 
When n= 1, (8.5) and (8.6) become 
C,%, + Corts — Diy, — Doya= — 4) 
Cy — Cove. — Dy, + Dore Gg7, — 4, 
and on using (8.9) and (8.8) we get 
Cy = 34% 2 + V2) ~ 34,(%1 ~ VY) + 3407 
Dy = 34% + V2) + 344(%) ~ Yi) + 34071 *1- 
Hence (8.3) and (8.4) give for the emergent intensities corresponding to linear 
emission 


Io, +p) = Jou *) =I (O)X(u) —I 7) V (4) + 34X(H) idol Xe + Ye + 711) 
(4, + Gy (%y ~ ¥4)} — 3H V (ue) ol%g + Ve + 7y%1) + (Gy + Gop — Ys)}s (8.11) 

U(r, —~ w= hye *) =I 7) X (4) ~ I (0) (4) — 34-X (4) dol %2 +2 + 71%) 
+ (Gy ~ Gop (%, ~ Y1)} + 34 V (4) lol Xe + Vo + 711) — (4, — Go 1 —Yy)}- (8.12) 
The constants J'(o) and J'(r,) are arbitrary as long as the Milne equation 
(7.2) alone is considered. ‘They can be evaluated by the method used by Horak 


and Lundquist in (8), of which the following is a modification. 
Let 


(8.10) 


1 1 
F(r) 2) llr u) de, K(r) | HAG pw) dy. (8.13) 


Then it follows from the equation of transfer (2.1), in the usual way, that 


ad 


d 
; K(r) =} F(r). (8.14) 
aT 


dr 


From these we easily obtain the relations 


F(r) = —4B,(r), 


F(r,) — F(o) 4| B,(t) dt, (8.15) 
0 


2[K(r,)— K(0)} — }7,[ F(7,) + F(0)] 
| 


. \" B,(t) dt |" By(e) dt} a (8.16) 


~ 0 


With the initial conditions given by (2.2), these can be written 


‘1 
\ ptT(o, +e) + Mr, — 4)} dp 
1 wr, 
/ | 1h fol) 4 h()} dp + B,(t) dt, (8.17) 
0 0 
1 l 
bf (4 bry){ U(r), — 2) — Ho, + 4) }dp = 4 (dre *)fi(4) ~SolH)} dy 
( “0 


i" {° B,(t) dt - I. Bt) at) dr. (8.18) 
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On putting B,(r) = dor + @,,fo(u)=f(u) =0, and substituting for J(o, +) 
and /(r,,—) from (8.11) and (8.12), we obtain two equations giving, respectively, 
J*(0)+J"(r,) and J(o)—J(r,) from which we obtain the following values ; 

J") = dol 3(%y — Yi )[%q + Va + b7y(%5 + Va) Ari} A 
+ 34171 V2 + d71%(%, + ¥1) — Xs + Ya) + 344, (3.19) 
Jr) 4y{3(% — Yi L%a + Yat b7i(%s + Ys) — yrs }/A 
+ 3Go[7 1% + h7y°(%, + Y1) + %y— Ya) + 34,4, (3.20) 
where A = X_+ Vo + 47,(%, + y,). (8.21) 
The equations (8.11), (8.12) and (8.19) —(8.21) give the final solution for linear 
emission. 

g. Isotropic illumination of the surface r=0(0<—m—1).—-Suppose that the 
surface r=o is illuminated by isotropic radiation of unit intensity. ‘Then 
1) = 1, f\(4) =0 and, if the atmosphere is not emitting, B,(r) =o. 

Hence (2.12) gives 
rl 
B(r) jer | exp ( T pe’) dy’ 
0 
- dx 
2 jor exp ( 1) 
\wE fr). 
Let O(r)=0%7r) when B(r) has this value. ‘Then the Milne equation (2,11) 1s 
D7) =aA,{O*(1)} + bw fr). (9.1) 
We could solve the general equation 
D"(r) = aA, {O"(t)} + hw, (7), (9.2) 
as has been done in (3) for the semi-infinite case, but little use is likely to be made 
of the solution for n>2. When n=1, the solution is easily obtained from the 
equation (4.3), but as E\(r)-- © as r-+0, this solution has no practical value. 
Since (d/dr)E,{r) E,(r), (9.1) gives, on differentiating with respect to r, 
d (d : 
; (7) wX, 4 7X0) + harD*%(0)E (7) — bor O*(7, )E (7, — 7) — bw Ey (7). 
( a. 


aT 
Hence by (4.3) and (4.4) 
dx 
= 


¢ *(7) — bm[D*(0) — 1] Je. ' jor) [ J (7, —7,%) 
dz ? 1 x = 1 


, ad 1x « dx 
wA, 4 7, Ot) | w[@%(0) | I(t, x) = + her%r,) | Sry t,x) — }, (9.3) 


i 
and.a solution of this is 


5 (7) — bw[O*%X(0) — 1] , Jd (r, ys + barr) : J(1,—7,x) - Oo. (9.4) 
aT 1 71 


x 
Multiply (9.4) by exp (—sr) and integrate with respect to r over (0, 7,), inte- 
grating the first term by parts. On introducing the r,-transforms of ®*(r), we 


have 


“« dx 
4%(r,) exp (—s7,)— 40) + bys) — h[X0)~ 1] { Rs, x)S 
. ie | ‘ 


+ }w@(7,) exp ( srs) | R( sx) oO, 
1 
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and by (4.12) and (4.13) this can be written 
bo%(s) = D*%(0)J (0, s) — D*(7,)J(7,, 5) —[J(0, 5) — 1]. (9.6) 
On writing —s for s and using (2.17) and (4.15), we get 
$,%(8) = D%(7,)J(0, 5) -D%(0)J(7,, 5) +I (74,8) — exp(—745), (9-7) 
and hence, by (4.21), 
dbo *) = D(0)X() — OX 7,) ¥(n) +1 —X(p), (9.8) 
by") = O47, )X(u) - 90) ¥(u) + ¥(4) —exp(—74/r)- (9-9) 


Put r=o and r=r, in (9.1). Since E,(o)=1, we have (proceeding as in 
(7.12) and (7.13)) 


‘1 
D0) = bor | dou) du + bor, (9-10) 
0 


rl 
*(7,) 7 jer | bu ') du + bw E(7;). (9.11) 
0 


On substituting for 4,%"') and 4,%"") from (9.8) and (9.9), we get 
Po) 1 — x9) + D7, )yv_ =o — Xp, (9.12) 
D*(0)Vv9 + D7, (1 ~ Xp) = Vo — FEA 7,) + bWEL7,) 
= Vo. (9.13) 
Hence, by (7.15), 
P*%(0)=xy, D*(7,)— Vo, (9.14) 
and (9.8) and (9.9) become 
by(w*)=1—(1 ~ X%)X(n) — yo ¥ (4), (9.15) 
b(e*) = yoX (4) + (1 — %) ¥(u) — exp ( — 7/4). (9.16) 
From (2.13) and (2.14), the emergent intensities are 
(0, +4) = bo") = 1 ~ (1 ~ %)X (4) — Yo ¥ (4), (9-17) 
U(r, — ) = bye") + exp ( — 71/4) =VoX (4) + (1 — X%) ¥ (4). (9.15) 
If we have an atmosphere emitting according to the law }wE,(r) and with no 


incident radiation, the Milne equation is again (9.1), but the emergent inten- 
sities are simply do!) and ¢,%'). These are given by (9.15) and (9.16). 
10. Isotropic illumination of the surface +r =0 (w =1).—When w=1, equation 
(9.4) has to be replaced by 
d iho d > dx 
q, PX") \(@*(o) 1) | I(r, x)= + j@%(r,) | J(r, 1,8) = =J(r), (10.1) 
where the r,-transforms of J(r) are given by (6.7) and (6.8). ‘Then s 'jo(s) has to 
be inserted on the right of (9.5), and (9.8) and (9.9) have to be replaced by 
bX u-?) = O40)X (ux) ~ ©%(r,) ¥(u) + 1 — X(u) + n[AX(u) — BY(u)], (10.2) 
4,4 u-?) =X r,)X (2) — @Y0) ¥(u) + ¥(u) — exp(—7/4) — n[ BX(u) — A¥(u)- 
(10.3) 


Equations (9.10) and (9.11) (with w=<1) are unchanged. On substituting 
from (10.2) and (10.3) and using the conditions x,= 1, yyg=0, we get 


Ax,—By,=0, Ay,—Bx,=0 
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and hence A= B=o. ‘The emergent intensities are therefore 
I(o, +p) = dof 2) =x ( '‘X(y) —-D Y() + I, ( 10.4) 
I(r, — 4) = b(*) + exp (— 7,4) = DX(n) — CY (yp), (10.5) 
where C = ®*%(0)—1, D=@*(r,). 
The constants C and D can be evaluated by (8.17) and (8.18). In this case 
Bz) = 0, fol) = 1, f(4) = 0 and it is easily found that 
C= —D= — fa, +y2+ 474(%,+y9,)]". (10.6) 
Equations (10.4) — (10.6) give the solution of the problem. 
For an atmosphere emitting according to the law }£,(r) and with no incident 
radiation, the emergent intensities are 4,%(°') and 4, "). ‘These are given 


by (10.4) and (10.5) but with different constants C and D, which have to be 
evaluated similarly from (8.17) and (8.18). 





APPENDIX 
11. The homogeneous Milne equation.—¥or the homogeneous equation 
J(r)=wA fJ(t)} (o<w <1), (11.1) 
we proceed as in the preceding sections. Differentiating with respect to r, we 


have 


\ 


d id ' , 
7g, 2) oX, GI} + bod (O)E(r) — bod (7, )E\(7, — 7). 


( 


On subtracting }w/J(o) times (4.3) and adding }w/(r,) times (4.4), we get another 
homogeneous Milne equation, from which it follows that 


d “ d. @ dx 
< I(r) jorJ(0) | Irsx) = + hod(rs) | Ir rx) = CI(r), (11.2) 


where C is an arbitrary constant. 
Multiply (11.2) by exp(—sr) and integrate with respect to r over (0, 7,), 


integrating the first term by parts. ‘Then by (4.7) and (4.10) 


; . d: 
J(7,) exp (—87,) ~J(0) + jo(s) — bad (0) R(s, x) _ 
Jy 


+ od rsdexp(~sr){ ° R(—5,2)F = Crd) 


By (4.12) and (4.13) this can be written 
jo(s(1 — Cx) =J(0)J(0, 8) —J(r,)}/ (74,8). 
Writing —s for s and using (2.17) and (4.15), we have 
lst + CH) =JI(r,)d(0, 8) —J(OI(14,5). 
Finally, on putting s=1/p and using (4.21), (11.3) and (11.4) become 
jhe "Yt — Cu) = (0)X(n) —J (7) ¥(v), 
j(e* 1 + Cu) =J (7, )X(n) —J (0) ¥(u). 
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If we now put r =o and + =7, in (11.1) and proceed as in (7.12) and (7.13), we 
get 


| 
J(o) jor | joe) dp, (11.7) 
0 


1 
Mry= bor | ile) de. (11.8) 


Hence on substituting for jj ') and j,(u°') from (11.5) and (11.6), we obtain 
the equations 


X(u rl = A 
H(o)41 4m oe dy} Ir) a| = aa non (11.9) 
“0 4 
ht Vu ‘ae | } 
I(o)hen | aor du tr) 4 so} — ve) MI ' (11.10) 
Unless C is a root of the equation, 
F(1/C) =o, (11.11) 
to . 
7 “ieee 
o =~p { 
“1 y 4 4’ 


Jgz— ys 


where 


be, (11.12) 
the only solution of (11.9) and (11.10) is J(o)=J(r,)=0. This gives j,(s) <0 
and hence J(r)0. 

We shall now show that F(z) =@(z), where 


E(2)=1 — hme ln- ——, (11.13) 
so that (11.11) is equation (6.3). Ifo<cm<1, then C= +k, and if w=1, then 
C =o, 

In (2) | have obtained the following relations between the X- and Y- functions 
(see (2), equations (2.4), (2.6) and (2.7) with f(z) = 4m): 
I 'sX(u) , |  exp(—7,/2 
L« a -4 | —— A> — - mt 
= XG)" J, sn Sf ~ XK) 
P21 1 { Se) 

: Y(z) y jo [ 5" 

ae Cin 4! 

X(z)E(2z)=1 jen | — : jw exp(—7, =) MH) ay’, (11.15) 


9 trp 


) 


dy’, (11.14) 


"he rt Xi ‘ ’ l 'sY() 
Y(2)@(2)=exp(—r,/2 iow “y ber | de. (11.16) 
0 


sp 


From (11.15) and (11.16), 
e¥(u), (s¥lu') fv vee tzX(u) , | 
jw? { ee ds ah at Me’ { X(2)@(=) t+)o | —— de Me 


Jo ®~p z+ = 
x { ¥(2)€(2)exp(r, s)-1+)e@ o| Pa X(p' a FM ‘lL 


oth J" 
Hence, by (11.12), 


5 2X 4 
F(s) X(2)(s) {1 — bor [ ( vay} ~ X(z) ¥(2)[€(2)]* exp (r,/2) 


o trp 


le " 
t ¥(2)@(=)exp(r,/2)4 ~ ber | eH) ay |, 


9 =~ 
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On multiplying (11.14) by X(z)¥(s)€(z) exp(r,/z), it is seen that the right-hand 
side of (11.17) reduces to Ez). 

First consider the case o<m<1. Then @(z) has simple zeros at s= +1/k 


and (see (2)) both X(z) and Y(z) have simple poles at z= +1/k. Let r,, r, be 
the residues at z= 1/k and let 
L=— =lim X(s) 
M. 


1 Vz) (11.18) 


The equations (4.22) and (4.23) hold when , is replaced by the complex variable 
z. Writing them in the form 


X(z) 12X(y') , , 1 z¥(u') I 
VG}! jo | a + ho [ s+ dy . Y(z)’ 
X(z),_'  ¥(u) f *sX(u) |  exp(—r,/z) 
vo} du + ¢ lo - ) du L Rn 
Va)!” |, 3 m ed |" . |, sp , j Y(z) 
and letting z+1/k, we have 


X(p') wi eee 
LA “for f 2 + hy’ dy’ 4 dl + hy’ du’ =0, (11.19) 


* Y(u#) X() | 
| 3 . a) 4. 
Lhe} . hy + {1 her [ hy: d f oO. (11.20) 


On comparing these with (11.9) and (11.10), we see that J/(r,)/J(o)=L if Ck, 
and J(o)J(r,)= Lif C=—k. If C=k, (11.5) and (11.6) give 

A ; 

ile ')= ey [X(n)—- LY(p)}, 


A 
= f "(us i)} 
: + hy, XH) Y(z)}, 


where A =J(o) and A is arbitrary. If C= —k, then 


Ayu ”) 


: B . , 
jl *) = —>— {LX(n) — Y(u)}, 
I+ ku 

RB (11.22) 
h(e-") = — hy [X(u)- LY(«)}, 

where B=J(r,) and B is arbitrary. But the sum of two solutions of (11.1) is 

also a solution and therefore the general non-null solution of (11.1) has the r,- 

transforms 


<. . B : , 
Jol Ya —— ~ (X(#) ~LY(pz)} +4 + Ie {LX(u)— Y(u)}, (11.23) 


ite ')= ra (LX(u)— Y(u)} + : ~ 7, (X(H) LY(y)}. (11.24) 
These are (6.5) and (6.6). 

Next consider the function J(r) corresponding to (11.23) and (11.24). Since 
js) is the ordinary Laplace transform of the function equal to J(r) in (0, 7,) and 
to zero in (r,, «), therefore 

J(r) = mil, es * Jol 1 exp (rs) ds (0 <r<7,), (11.25) 


where the line R(s) = lies to the right of all the poles of j,(s)/s. 
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Near z= 1/k, 


X(z)= tS +a(z), Y(z)= —, + (2), (11.26) 


where x(z) and y(z) are regular at z=1/k. Hence, by (11.18), 
X(s")—LY(s"")=x(s-")— Lys"), 
r2—r,* ks 
eI ath Se ee Hi adant 
LX(s*)— Y(s-"")= 7 ee Lx(s'*)— y(s"*). 

Thus X(s°')— LY(s"') is regular at s=k, while LX(s'')— Y(s"") has a simple 
pole there. It now follows from (11.23) with » replaced by 1/s that j,(s)/s has a 
simple pole at s= k except when A and B are such that the coefficient of (s—k)"' 
is zero, On evaluating J/(r) from (11.25) by contour integration, we get a mul- 
tiple of exp (kr) arising from the residue at s = and this is the term of greatest 
magnitude for large r. 

When w=1, the working is simpler. Since Co, equations (11.9) and 
(11.10) become 

J (01 ~ %) +4 (7 )¥o =, 


J (0)vg + J (4,1 — Xp) = 0. f (11.27) 


Let X(u), Y(#) denote any solutions of the X- and Y-equations and let X,(,), 
Y (4) denote the standard solutions. For X(u) and Y(y) we have the relation 
vy + ¥o = 1 and the equations (11.27) give 
JS (r,) J(0), (11.28) 
and so (11.5) and (11.6) give 
Jol ') A{X(p) t Y(“)}, 
: > (11.29) 
ie ')= — A{X(u) + ¥(w)}, J 
where A=J(o) and A is arbitrary. For Xu) and Y,(u) we have x,=1, yy=0 
and the constants J(o) and J(r,) become indeterminate by (11.27). 
Then 
jhe *)= AX (u)— BY (x), 
: ; (11.30) 
ju )=BX(u)-AY (pn), J 
where A=J(o0), B=J(r,) and these are arbitrary. Since X(yu)+ Y(u)=X fu) + 
Y (4) (see Section 1), the solution (11.29) is only a special case of (11.30). In 
the body of the paper we have dropped the suffix s and (11.30) agrees with (6.7) 
and (6.8). 

In (2) I have shown that X,(z) and Y,(z) are O(z) for large z. Hence j,(s) 
is O(s ') for small s and j,(s)/s has a double pole at s=o. ‘The contribution to 
J(r) arising from the residue at this pole, when evaluating (11.25) by contour 
integration, is of the form ar + 8 and so J(r) is O(r) for large r. 

We now turn to the question of the choice of J(7, ) in Section 4, when w= 1, 
and the corresponding ambiguity in X(y) and Y(,). 

The equation (4.1) has the solution 

J(r,0)+CJ(r), (11.31) 
where ( may depend upon a. Is it possible to choose C so that the subsequent 
analysis will hold with J(r,0) replaced by (11.31)? By (4.7), R(s,c) will be 
replaced by 


R,(s, 2) = R(s,a) + Cs Yofs), 
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and in order that we may have 
R,(s,0)= Re, 5), 
we must have C = A’o'j,(a), where A’ is independent of o ands. Hence 
R,(s,0) = R(s, 2) + A's or Yols\jolo). (11.32) 
If jo(*) is given by (11.29), then 


jos) 1(9) = Joos) 
and it is easily verified, with the aid of (2.17), that the condition 
exp (—s7,)R\( — 8, 0) = exp (—o7,)R,( — 2,5) 
is also satisfied by (11.32). Hence J(7, 0) may be replaced by 
J(r,0)+ A’a Yla)J (7), (11.33) 
where j,(o) is given by (11.29) and J(r) is the corresponding source function, 
It now follows from (4.21) that X(y) and Y() may be replaced, respectively, 
by 
X (pe) + A’pjle "J (0) and V(x) + Apfel *)A(7,). 
Hence, by (11.28) and (11.29), X(y) and Y(yz) may be replaced by 
X(u) + Qe[X(H) + V(H)), YC) — Ov [X(H) + ¥(H)], 
where O= AA'J(0). By choice of QO, these functions are X,(u) and Y (4), and 


then (11.33), with A’ = O/J(0)A, becomes the correct J(7, 0) for the analysis of 
section 4- 


St Hugh's College, 
Oxford : 
ross July 26 
Re fere nCceS 


V. A. Ambartsumian, Ast. 7. of the Soviet Union, 19, 30, 1942 

I. W. Busbridge, Ap. 7., 122, 327, 1955 

I. W. Busbridge, Quart. 7. (Oxford) (2), 6, 218, 1955 

I. W. Busbridge, M._N., 118, 661, 1955 

S. Chandrasekhar, Radiative Transfer, Oxford, 1950 

S. Chandrasekhar, D. Elbert and A. Franklin, Ap. 7., 115, 244, 1952 

E. Hopf, Mathematical Problems of Radiative Equilibrium, Cambridge Tracts, No. 31, 
1934 

8) H. G. Horak and C. A. Lundquist, Ap. 7., 116, 477, 1952, and 119, 42, 1954. 

(9) V. Kourganoff, Basic Methods in Transfer Problems, Oxford, 1952 


é 





THE OXFORD 35m SOLAR TELESCOPE 
H. H. Plaskett 


(Received 1955 October 7) 


Summary 


A new telescope for solar work with an equivalent focal length of 45:2 m 
is described. It is an ordinary Cassegrain reflector used at the coudé focus 
The optical requirements for satisfactory performance are listed and som« 
of the mechanical features briefly described. The performance is assessed 
from a Hartmann test 


Introduction.-Yor some time there has been a need at Oxford for a solar 
telescope giving a larger scale than the existing instrument (4). As a result of 
savings made during the war and a generous benefaction from Professor D. A. 
Jackson, enough money was available to build a high-dispersion spectroscope, 
and in 1948 the University Grants Committee made available a sum of {16,300 
for a telescope and a building to house it and the spectroscope. Specifications 
were drawn up for an enlarged version of the 19°8 m telescope with 81 cm diameter 
flats for the coelostat and secondary mirror. Fused silica blanks of this size 
were not, however, available in this country at that time, and a provisional 
estimate indicated that two blanks of this size in low-expansion glass would cost 
{2000-—too high a proportion of the money available. Accordingly it was 
decided to depart from conventional solar-telescope design, and to build an 
{70 Cassegrain reflector with a principal mirror 51 cm in diameter. 

Optical requirements.—For a variety of reasons it is essential in a reflector 
designed for solar work to have large-aperture Cassegrain and coudé mirrors. 
These act as central stops in the optical system and affect the diffraction pattern. 
Theodore Dunham* has recently calculated the distribution of light in the 
diffraction disk for various relative diameters of the central stop. With his 
permission | quote in Table I under each central stop first the value of z,/n 


Taste I 
Diffraction patterns 


Relative diameter of central stop 


‘ ™ rt) 


Z,/77 
ad od | he /o 


Dark ring 1 83:8 


1°0°o 47°35 o'9Os } o’S9 26° 

2°29 83:0 2°16 2°05 56°0 
Dark ring 3 ‘ 93°5 3°50 go°2 3°39 3°21 76:2 
Dark ring 4 4°2 95°3 4°13 go"4 4°52 D°2 438 86'5 


Dark ring 2 : g1°o 


where 2,A/7D is the angular radius of the dark ring, and secondly the percentage 
of the total light of the image which is concentrated inside the dark ring. For 
no central stop Rayleigh arbitrarily defines the observable image as that part of 
the Airy disk which lies within the first dark ring and contains 83-8 per cent of 
the light. If the relative diameter of the central stop is o-6 we get a little more 


* An abstract appears in 7.0.5S.A., 41, 290, 1951. 
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than Rayleigh’s percentage of light within the third dark ring, and consequently 
an observable image with angular radius 


3°39\/D=0"-78 


for a diameter, D, of the principal mirror of 50°8 cm and a wave-length of 5600 A. 
In solar work we cannot afford to have a resolving power smaller than 0°-78 
(comparable with the radius of the tremor disk in superlatively good seeing), 
and consequently the effective diameters of the Cassegrain and coudé mirrors 
must not exceed 30°4.cm, that is 0-6 of the diameter of the main mirror. 

For reasons of economy the dome has to be kept as small as possible, and this in 
turn means that the focal ratio of the main mirror must not exceed f/5. At this 
focal ratio the flux of radiation in the focal plane of the main mirror, that is in 
the vicinity of the Cassegrain mirror and its supports, is sufficiently intense to 
make the paint-work smoke and thus to produce bad seeing in the telescope from 
the convection currents resulting from this local heating. ‘lo eliminate this 
possibility the Cassegrain must have a clear aperture so large that no matter what 
part of the solar image is on the axis of the telescope, that is on the slit of the 
spectroscope, the extra-focal patch from the main mirror will never fall off the 
reflecting surface of the Cassegrain mirror. In this way, apart from the minute 
fraction absorbed by the aluminium coating, all the high-intensity flux from the 
main mirror is reflected by the Cassegrain and none falls on the rim of its cell 
or its supports to produce local heating. ‘lo meet this requirement the Cassegrain 
mirror must have a minimum clear aperture of 21-7cem. ‘The actual clear 
aperture given the mirror is 22:1 cm (8-70in.) with an outside diameter of the 
cell of 24°1cm (9°50 in.). 

In a telescope designed for solar spectroscopy it is necessary only to ensure 
full illumination of the image over the height of the slit. As indicated above, 
the Cassegrain mirror is (though for quite other reasons) more than sufficiently 
large to give full illumination over the whole solar image. ‘The limitation on full 
illumination comes from the first coudé mirror, distant 11-94m (470-265 in.) 
from the focal plane of the //70 telescope. The diameter of the extra-focal 
patch on the first coudé for an axial image will therefore be 171m (6°73 in.) 
To ensure full illumination over the height of the spectroscope slit, approximately 
5 cm, the coudé mirror must be large enough to accept the full extra-focal patch 
for a central ray making an angle of 2'-35 with the optical axis. ‘This requires 
a minimum clear (projected) circular aperture of 17-7em. Further, with this 
clear aperture points on the image radially distant as much as 16'-3 from the 
optical axis will suffer only 12-2 per cent vignetting; in other words the limbs 
of the centred solar image in the focal plane will always be clearly visible. 

The first coudé mirror which directs the light from the Cassegrain down the 
polar axis is used at angles of incidence as high as 56 -75, corresponding to a 
declination of +23°:50. ‘The reflecting surface is therefore given an elliptical 
cross section with a major axis 32°8cm (12-goin.) and a minor axis 17°8cm 
(7-00 in.), that is, just more than sufficient to give full illumination over the whole 
height of the spectroscope slit for this declination and a fortiori for all smaller 
declinations. A point to be watched is that with decreasing declinations the 
ends of the major axis of the coudé mirror move out of the shadow of the Cassegrain 
cell and reflect parallel light direct from the Sun into the focal plane. To avoid 
this, lobes are attached to the north and south sides of the Cassegrain cell, giving 


33° 
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to it an elliptical cross section with a major axis of 29-2 cm (11°50in. ) and a minor 
axis of 24:1cm (g’s5oin.). The shadow of this enlarged Cassegrain cell prevents 
parallel light from the Sun reaching the coudé, no matter what part of the solar 
image is on the slit, provided 86> —21°-8. This declination limitation could be 
removed by slightly increasing the lobes on the Cassegrain cell, but at Oxford it 
is no limitation since at these declinations the Sun has too small an altitude for 
useful observation. 

The elliptical central stop with axes of 29:2cm and 24:1cm given by the 
Cassegrain cell corresponds closely enough, as far as the diffraction image is 
concerned, to a circular stop with a diameter of (29-2 x 24°1)"*=26-5cm, or 
052 of the diameter of the main mirror. The effective central stop of the 
telescope is, however, provided by the coudé mirror and cell. ‘The mirror, one 
fifth of the way towards the focus of the paraboloid, provides a central stop 
25 per cent larger than its geometrical dimensions. ‘The edge of this mirror 
towards the Cassegrain is cut away, so that for = +23°:5 the mirror and cell 
together present an effective central stop of 18-8cm (7:40in.) diameter which is 
lost inside the central stop of the Cassegrain cell and its lobes. The axis about 
which the coudé rotates lies in the plane of the reflecting surface so that, as the 
declination decreases, the mirror remains centred in the beam, but the north and 
south edges project into the beam. ‘The extreme occurs at = — 23°:5 when the 
effective stop presented by the coudé has a major axis of 37°2cm (14°64in.) and 
a minor axis, provided by the Cassegrain cell, of 24°1cm (g°50in.). This 
corresponds closely enough, as far as diffraction is concerned, to a circular stop 
of the same area and a diameter of 29°4cm(11°7gin.). ‘This is 0-59 of the diameter 
of the main mirror, and is just within our maximum allowable central stop which 
gives a diffraction image of 0”-78 angular radius. 

The fourth and final mirror, the second coudé, is not in the telescope proper 
but is mounted in the south prolongation of its polar axis. It accepts the beam 
from the first coudé at an angle of incidence of 45° and turns it to run horizontally 
to the west. ‘This second coudé is some 3 m from the focal plane of the telescope, 
that is the slit of the spectroscope, and has to be large enough to give full illumina- 
tion over the full height of the slit, which is easy, but also large enough to accept 
at least part of the extra-focal patch from the limb of the centred Sun. ‘This is 
necessary both for guiding and for finding the heliographic coordinates of the 
point under observation on the Sun. The second coudé is given a clear aperture 
of 38-1 cm (15 in.), and this is more than sufficient, except in the plane of incidence, 
to accept all the light from the first coudé. The limitation of aperture in the plane 
of incidence does not prevent the limbs of the Sun being seen in these position 
angles, but they are highly vignetted. 

Telescope design.—A telescope to meet the foregoing requirements was 
designed and constructed by Messrs. Grubb, Parsons & Co. under the supervision 
of its managing director, Mr G. M. Sisson. Plate 7 shows a photograph of the 
telescope in position in the dome, Plate 8 a photograph of the tower with the 
spectroscope tunnel seen stretching westwards and Fig. 1 is a meridian section 
through the tower showing the general arrangement. ‘The telescope is designed 
to do just what is required but no more; hence the large polar axis to pass 
unimpeded the beam from the first coudé, the short tines of the fork since the 
telescope is only to be used within declination limits + 23°-5, and the dome shutter 
pivoting at the top. 
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The principal mirror is figured on a fused silica blank of 54-6cm (21-5 in.) 
diameter and 8-gcm (3:5oin.) thick. It is supported on six “wobbly” pads 
mounted in pairs at the vertices of a triangular frame. ‘This frame is adjustable 
for squaring-on, and the three edge supports are made as far as possible friction 
free so that the mirror does not bind when this adjustment is being made. The 
supports and the mirror are mounted in a massive cell which is bolted to the centre 
piece. ‘I'wo trunnions (one of which can be seen in Plate 7) in the side of the 
cell fit into a yoke carried by a block and tackle from the dome. The mirror and 


cell can thus be easily lowered through a trap door down to the ground floor of the 
tower for aluminizing. 
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Fic. 1 1 diagrammatic sketch to scale of the whole telescope. The height from the basement floor 
to the top of the dome is 10°2 m, or nearly 44 feet 


‘The Cassegrain mirror cell and its focusing mount are supported by webs 
from a stubby and very stiff structural steel-work tube. The mirror is figured 
from a fused silica blank of 23°2cm (g°12in.) diameter, 4:1cem (1°62in.) thick. 
The cell which carries it is provided with squaring-on screws. 

The first coudé mirror is carried in a cell on a stout spindle through the 
declination axis. ‘lhe half-speed drive is provided by phosphor-bronze tapes, 
spring loaded to a high tension, and moving over cylindrical rollers. ‘The mirror 
with an optically plane surface of elliptical cross section (axes 32°8 and 17°8 cm) 
is made from a fused silica blank 60cm (2°35 in.) thick with its leading edge cut 
away so as to offer the minimum central obstruction to the light beam. Many 
of the problems of design and construction of the telescope seemed to centre about 
the coudé mirror. Its awkward shape made figuring with full-size tools very 
difficult, and that there was ultimately a succesful outcome must be attributed 
to the skill and experience of Mr G. E. Manville, the chief optician of the firm. 
Mechanically too the large mass of the coudé and cell increases sensitivity to 
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deflectional and torsional oscillations of the spindle. When in mid-1953 the 
mirrors were finally installed in the telescope, low-period oscillations of the coudé 
with an amplitude of a few seconds of arc were in fact found. ‘They appeared 
to arise from heavy traffic on a nearby road. These oscillations were finally 
eliminated by Mr Sisson, partly as a result of increasing the spring loading of the 
tapes and so increasing the natural frequency of the system, and partly by suppor- 
ting the spindle near the cell by two flat steel strips carried from the centre piece. 

The telescope is driven by a synchronous motor through a large worm and 
focm diameter worm wheel. ‘The long focal length of 352m is a severe test 
of the accuracy of the cutting and the mounting of these parts, but in actual use 
no evidence has been found of periodic or other errors in the drive. By the 
usual differential gear a guiding motion (approximately diurnal rate) can be 
applied in right ascension, or a setting motion of ten times this rate can be brought 
into action by an auxiliary motor. ‘These gears are enclosed in an oil bath. 
The guiding and setting motions are controlled from a panel at the focal plane; 
this panel also carries buttons for a guiding motion in declination, for focusing 
the Cassegrain mirror and for rotating the dome. 

The second coudé mirror is figured on a fused silica blank of 41 cm (16-1 in.) 
diameter, 7°6cm (3in.) thick. The mirror and cell have an altazimuth mounting 
with fine motions. Like the telescope proper this mounting is fixed to a concrete 
pier, isolated from the building, and erected on a large concrete mat 2°6m below 
ground level. 

Hartmann test.—'Vhough the associated long-focus spectroscope has not yet 
been completed, the telescope has been in use since 1954 July for low-dispersion 
spectroscopy and interferometry of the solar chromosphere. ‘The large solar 
image of 32°4cem (12°8in.) diameter with a scale of 1mm-=s5"-9 is well suited 
for such work and for spectroscopic studies of solar granulation, as well as many 
other problems, and such investigations are likely to be aided by the brilliance 
of the image resulting from an effective focal ratio of {//88. ‘The rotation of the 
image, though an inconvenience for some purposes, is a gain in others, and can 
in any event be arrested by an image rotator. No changes of focal length are 
observed with continuous exposure to the Sun, showing both the effectiveness 
of the low-expansion fused silica mirrors and the small amount of heat absorbed 
by the aluminium coating of the Cassegrain. A noteworthy feature is the small 
depth of focus, and this coupled with solar definition varying from moderate to 
good suggests that performance is limited by the seeing rather than the optical 
properties of the telescope. 

To find the latter, Hartmann plates were made in 1955 April of Procyon, 
Regulus and Arcturus. The Hartmann diaphragm had holes 
opposite each other in pairs on one or two azimuths per zone. ‘The holes were 
made as small as possible (5-1 cm diameter), only 1/7ooth the focal length instead 
of the minimum 1/400th recommended by Hartmann (2), but even so several 
trials were made before a satisfactory diaphragm could be constructed in which 
the holes were not obstructed by the Cassegrain webs, or the supports of the 
coudé spindle. ‘The positions of the holes are shown in the first two columns 
of ‘Table Il, where r is the radius of the zone in cm and ¢ is the azimuth measured 
counterclockwise from the east declination axis of the telescope. ‘The remaining 
six columns, one for each photographic plate, give the residuals f—/f in mm, 
where / is the focal setting for a given r and ¢, and f is a provisional mean focal 


diametrically 
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setting for the plate as a whole. ‘The column headings give the star, the angle of 
incidence, /,, of the first coudé, and the angle, 4», of the west side of the plane of 
incidence of coudé II measured counterclockwise from the main-mirror side of 
the plane of incidence of coudé I. 


Taste Il 
‘alues of ({ —f) in mm as a function of r andd 


Procyon Regulus Arcturus 


7 j 5 54°7 
7°3 7's 42°3 


5° . 375 
130 55 45°7 
75 at 34°7 2 14°9 
15 19 y , 22°'5 
105 49 45°! 5 44°7 


In an ideal telescope these residuals would be constant within the errors of the 
determination. An inspection of ‘Table II reveals, however, two systematic 
trends. First the focal settings for azimuth 105° are always greater than for 
azimuth 15°, indicating astigmatism, and secondly the focal settings steadily 
increase from the innermost to the outermost zone, indicating zonal aberration. 

The most likely origin of the astigmatism is to be found in the two coudé 
mirrors since these may have a slight curvature. If such a mirror has a radius of 
curvature p, its focal length in the plane of incidence is }p cos1, and perpendicular 
to the plane of incidence is p, 2 cosi, where 7 is the angle of incidence (5). Inany 
other azimuth (3) making an angle « with the plane of incidence the focal length is 


I 2 /cos* x . a 
* sin” ~ Cosi 
T . p cost p 


Knowing the focal lengths in terms of the radii of curvature, p, for coudé I and 
ps for coudé II, the object distance, /, for coudé | and the distance c between 
coudé | and II, we can apply the paraxial formula to find the image position for the 
azimuth «. lor another azimuth § the focal length of a given coudé will be 


I 2 /cos* p a a 
+ sin* 6 cost 
ly p cosi 
Denoting the difference in image positions for these two azimuths by y, and 
neglecting squares and higher powers of the small quantities p, ' and p, ', we 
find an equation of condition 


apy §+tp_q = y 


where o=(b,—a,)l? + y(a,+,)/ r= (a,—b,\1+cP —~ylag+b,\l +c), 

the suffix 1 referring to coudé I and 2 to coudé I]. From zones r= 16-3 cm and 
21*8cm we thus have 12 observation equations from ‘Table I] to solve by least 
squares to find the curvatures of the two coudé mirrors. The coefficients o and r 
may be readily calculated for each equation from the data given in Table II, 
bearing in mind that the angle of incidence, i,, for coudé I] is always 45°. Solving 
we find 


py t= —(1'24+0°28) x 10 * cm"! ; Ds + (0°00 + 041) x 10°* cm. 
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Coudé II is thus flat within the errors of measurement, and coudé I is convex with 
a radius of curvature of some 8km—a good “ flat” on any ordinary standards. 
Unfortunately coudé | is used at high angles of incidence and is at a considerable 
distance from the focal plane. In Table III are given the differences, y,, in focal 
setting for the resulting line images and the path differences in units of A= 5600 A, 
computed from the appropriate formula, for various declinations. It is unlikely 
Taare Ill Taste IV 
Effect of convexity of coudé I Zonal aberration 
6 iy Ye € f,—f ai | T’/p\ 
mm mm mm 
20 55 42°0 1°2A 16°% 45°0 0°22 1°6 
o 45 25°4 o°7A 19°3 9°8 0°05 o"'4 

20 35 14°4 o°4A 21°8 34°6 o’2! 1°6 
that solar seeing is sufficiently good for astigmatic path differences of as much as 
A/2 ever to impair the definition or affect the contrast, and it is at least fortunate 
that the best solar seeing at Oxford occurs in the early spring and the late autumn 
when the astigmatism from Table III approaches this limit. Nevertheless there 
can be little doubt that an improvement in the performance of the telescope, 
particularly in the three summer months, would result if the radius of curvature 
of coudé I could be at least doubled. 

To find the zonal aberration we must eliminate from the focal settings in 
Table II the known effect of astigmatism. The minimum focal setting, /,(r), 
for any zone r occurs at right angles to the plane of incidence of coudé I, that is at 
d@=o. At any other azimuth we shall have 


lny=hlr)4 Yq sin® d. 

Knowing y, for the angle of incidence 7,, as computed for example for Table III, 
we can find for each entry in Table II the value of /,(r). Then f,=f,(r) + by, 
is the desired focal setting for the zone r which is independent of the astigmatism. 
Taking the mean of these values /, for each zone on a given plate and selecting 
a focal setting f for which the disk of least confusion is a minimum, we can find 
mean values of /,—/ for the six plates. ‘The results are given in Table IV, where 
T’ =r(f,—f)/F is the lateral aberration and 7"/p is its ratio to the expected Airy 
disk inside the third dark ring. ‘The path differences are giver by 


T'dr, 


¢= 


I Tr 

PS ses 
and are found never to exceed 0-2 A, that is just inside the Rayleigh limit. Further, 
since 7”/p is comparable with unity, both of the criteria of Danjon and Couder (1) 


are satisfied and the telescope may be regarded as having negligible zonal 
aberration. 


In conclusion | wish to express to Mr Sisson, Mr Manville and their assistants 
at Grubb, Parsons my appreciation for the thought and care which they have 
given to making the telescope a success. ‘lo Mr E. H. Swainand Mr J. M. Miller, 
former University Surveyors, and to Mr T. R. Davies, the Clerk of Works, 
I wish to add my thanks for the design and construction of a tower and tunnel 
perfectly tailored to fit the telescope and its accompanying spectroscope. 

Umiversity Observatory, 


Oxford ; 
1955 October 6, 
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THE SECOND-ORDER THEORY OF THE FIGURE OF JUPITER 
P. J. Message 
(Communicated by Sir Harold Jeffreys 
(Received 1955 April 12) 


Summary 

The second-order hydrostatic theory of a rotating planet is applied to 
each of the four models of the mass distribution in Jupiter proposed by 
W. H. Ramsey and B. Miles. The coefficients in the external field of the 
planet are found in each case. An inequality relating these coefficients 
given by Sarmmpson for the ellipsoidal planet is extended to cover the case of 
a planet with second-order departure from the true ellipsoidal form, it 
being shown that a real planet must be of this type. Using the general 
nature of the dependence of the coefficient of the fourth harmonic on that 
of the second found from the models, and the value of the latter coefficient 
derived from observations of the great satellite motions, an estimate of 
the secular motion of the node of the fifth satellite is made which is close 
to the observed value. 


1. The hydrostatic theory of a rotating planet to the second order has been 
treated by Sir Harold Jeffreys.* The planet is taken as symmetrical about its 
axis of revolution, the general surface of constant density having the equation 

4 pt 


rer'\1—~ Set +(e’ + e*)L, + (fe* + 4x) L, + Ofe*)}, (1) 
where, with ¢ the distance from the centre and 4’ the latitude, 
L,= 4 — sin? ¢’ qP,, (2) 
: 2 
L,=sint ¢' --°sin? d+ 2= SP, 
r’ is the radius of the sphere of equal volume, equal to 4, say, for the outer 
surface of the planet. e’(r’) is related to the ellipticity e of the surface defined 
by (1) by 
e ” in > (3) 
and its square may be interchanged with e*, to the accuracy desired. «(r’) 
measures the departure of the surface (1) from the true ellipsoidal form. 
With w the angular velocity of rotation, M the mass of the planet, and / the 
constant of gravitation put 
yi w*h* f Vv. (4) 
The potential function is then given by 
4 5 


‘. } 3 /Sr,° re 3 Pr! 
pd(y *) + - Pot ( = + PT) Ly + (3 r 


r, 5 


» 47 2} 
-" Sia), 


* Sir Harold Jeffreys, M.N., 113, 97-105, 1953. 
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where s (* ee 
Po= -3] pd(r®)=p for r,=b 
Ty Jr=d 


, 
: pd{r’(e’ t te*)}, 
“g =i) : 
b 


pd(e’ + ie*), 
r=¢ 


Sd fy’7(p2 4 8) 
r? Poss ‘(e + ae) }> 
h 


| 

| 

O=r? | pd(5). | 
T= ’ } 


The condition that the surface (1) is a level surface leads to de Sitter’s equations, 
as corrected by Jeffreys: 
(e’ + =e" )py (S+T)— bpm = Se(pm 37°), (7) 
(3¢? — 8x)p, — 6e’'S + 3P +20 =o. (8) 
2. The mass distribution models.—W. 1H. Ramsey and B. Miles* have given four 
models for the distribution of mass within Jupiter. ‘They were found from the 
theoretical relation between the pressure and density of solid hydrogen, allowance 
being made for admixture of helium, the proportion varying between the models, 
In the models J,, J,, J,, the heavier elements are supposed concentrated at the 
centre, J, having the most massive central particle and J, the least. J, has no 
massive central particle. Sir Harold Jeffreys allowed me to use values of the 
density in the models which had been sent to him by W. H. Ramsey. I found it 
necessary to extrapolate the values to values of r’ smaller than those for which it 


had been given. ‘The form of the compressibility used in the extrapolation was 


k a. 


i =ap + ho, (9) 
ap 


as given by W. H. Ramseyt, where « varies from 5 3 for a perfect Fermi gas 
to 5 for heavy atoms. «=2 was adopted as being reasonable for hydrogen, 
and leading to a case of Laplace’s hypothesis : 


kx p*. (10) 
In the process the mass of the central particle in J, was re-estimated as 
(47/3) x 70°9 x 10%’ gm. ‘Those in J,, J, were unaltered. 

For J,, J, and J, the density is given to three decimal places for the greater 
part of the planet, so the uncertainty of its specification may be taken as + 0°oo1. 
At the surface po-1o, and is thus specified to +1 per cent. For most of the 
planet p=1, and so is specified to +o-1 per cent. In the case of J, only two 
decimal places are given, so the surface value is specified to only + 10 per cent. 

3. The ellipticity.—lf we put 

de’ 
(11) 


two differentiations of (7) lead to 


P dy . r dp, 
r ; + bY , = 10A, (12 
dr , f* dy ) 


*W. H. Ramsey, M.N., 111, 427-447, 1951; and B. Miles and W. H,. Ramsey, M.N., 112 
234°243, 1952 
+ W. H. Ramsey, M.N., Geophys. Suppl., 6, 42, 1950 
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where 


2e' , dy 3 2 1 dpo 1 — 
= — (17 + an) gy +2’ + 130 +t — qi (O+ 70-20 Vf (13) 


which in turn leads to de Sitter’s equation : 


—— 
Vit +n) = se | por’ F(n)dr’, (™4) 
Po/o 


r 
where F(n) = (1 + 4n — fn? + 4)/V (1 +9). (15) 
An examination of the influence of A on 7» shows it to affect it by o-2 per cent or 


less, so it was neglected. 
For the models with a central particle we have, if m* is its mass, 


3m* 


Py™ 4nr® 


as >O, 


which, substituted in (12), gives 


, 


yr3 as r->o. 


For the models J,, J,, Jy, (14) was solved by iteration, using this initial value of ». 
The solution is given in Tables I-IV. ‘The values for r' =6 thus obtained were 
(the suffix “‘s”’ indicating surface values): 


iF Js Js 
"hs 1658 1455 1*306 
Where there is no central particle (12) shows 


, 


yO as r—>0. (13) 


‘The integration for J, with this initial value has been carried out by Ramsey and 
Miles. They obtain 7, = 1-099. 
Differentiating (6) and eliminating S, with the original equation leads to 


5m 5(2%5" — He + 5)m* 


bs aa T 
2(2 + %) 21(2 +7)" (19) 


De Sitter*® gives the ratio of the mass of the system of Jupiter to that of the Sun 
as 1047°40+0°03. Sampsont, from the semi-durations of satellite eclipses, 
gives the angular equatorial semi-diameter of Jupiter at mean distance as 18"-g2. 
Micrometer measures give about 19"-2, but these are probably increased by 
irradiation. I take a,=18"-92+0":20. The period of rotation is observed from 
spots to be g" 50" 30° at the equator and g"55™"5* to g"55™ 42° in temperate 
latitudes. By analogy with the Earth there is a strong atmospheric current 
at high altitudes, in the same direction as the axial rotation, and as the spots 
observed are almost certainly atmospheric phenomena the period of rotation 
of the solid planet is likely to be nearer the larger of these two. Hence I take 
g"ss"+2™s5. With Struve’s value 0-065 for the ellipticity, (5) gives 


w*a,*( I @s) 


{M 


m 


= 0°08339(1 + O’o10). 


* W. de Sitter, Proc. Roy. Akad. Sci, Amst., 10, 710-729, 1908 
1 R. A. Sampson, Harvard Annals, §2, Pt. 2, 334, 1909 
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(19) then gives : 
J J, Js iP 
€s. = 0°05729(I + O°011), 0°06013(1 + 0°01 1), 0°06339(1 + 0011), 0°06763(1 + 0-020), 
and hence by means of the values of e, from (3), m is corrected to: 
m = 0°08404(1 + 0°010), 0°08379(1 + 0°010), 0°083.49(1 + 0°010), 0°08318(1 + O-o10), 
These values give from (19), correct to the second order, 
€s, = 0°05774(1 + O'O11), 0°'06042(1 + 0°01 1), 0°06347(1 + 0'011), 0°06739(1 + O-Or!). 


From (11) we have a, 7 

(11) Inés’—Ine’= | dr’, (20) 
rv! 
whence e'(r’) is determined throughout the planet, the values being given in 
Tables I-IV. 

4. The departure from the ellipsoidal form.—Put 
1” 
P,=- 7 


rT; ’ 


pd(r’’e*), 


Fw. " 
P,= > | pd(r''«); 


ry r ) 


then equation (8) becomes 
I 
«(r’ = K (r’) ' ~(P,+4 VY), 
, 3pol? ) . 
where , I T 
K,(r') = 3) 624 (P,—2e'S)>. 
ra | 
Since P,, QO are linear operators on the function «(r’), we may write (22) as 
(rx!) =,(1') + L(x) (22.4) 
where # is a linear operator. Put 
84x =L(x,), 
b,K=L(8, x) form 1. 


Then if A is the greatest eigenvalue of 7 we have 


(24) 


5/0, eA as n> D, 


If |A|<1, we have 5,«-+0 as n-> », and hence, putting 
Kn : (25) 
K,—k,~L(x,)= —38,,.« which 0 a8 n-> ow, 


so that since «x, converges, it converges to «, the solution of (22). It was found 
that 5,«/5,« was constant to the desired accuracy (5,«/5,« in the case of J,), so 


(26) 
Therefore KOK, + By + Bye + Bye /(1 —A). (27) 
An examination of (22) as r’-»o shows «(o0)=0. Hence the functions «(r’) were 
found, and are given in Tables I-IV, the surface values being 
iF J Js I, 
r/A= 615 57 5°2 5°0 
10%x,= 5°76+0°13 5°33 + 0°12 4°65 + 0°10 P59 O15 
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Tasie I Tassie If 
Model J ,, central particle 15°5 per cent of Model J,, central particle 10°5 per cent of 


total mass total mass 


10'K r’ p n 
10°cm gem“ 
658 0°05774 5°76 ‘99 095 1°485 006042 
642 05719 ‘60 ‘9s 132 1°466 o-osgg! 
604 0°05653 544 °1BS5 1°426 005929 
S12 0°05400 “78 261 '353 0°05692 
449 0°05164 ‘20 , 393 243 0°°5473 
‘354 0°04942 3°73 335 1°1§8 0°05271 
261 04735 35 361 046 0°05090 
04550 i! "007 Sa 
04252 73 024 C¢ 0°04682 
03940 R35 | 245 0°04376 
03600 2 4° 451 04072 
03255 76 722 « 03764 
g69 03442 


02535 24 2°222 


©°05007 


02905 


23095 
02133 2490 | O2714 
01727 73 a°774 1 02291 
O1313 47 o82 01827 

1§0 0°00907 23 2°C 432 1°92 01326 
b*4ll ©'OOS45 o°10 863 2 00533 
2633 00264 03 2 4°45 
2835 OC 0005S oo 


1956 o'00011 


00417 
52 2 00135 
Oooo It) 


3,000 00000 00000 


Tasre Ill TAaBLe I\ 
\lodel J, central particle 5:2 per cent of Model J , 


total mass yr’ p 1) ‘ 
1o*cm gem 


no central particle 


6°99 o'll 099 0°060739 

064347 6°96 °o 065 0°06705 
06400 6°9 2: 034 0°06647 
06243 6°35 7 977 O'06550 
06026 6 
oss25 6 
0°05050 6 


05495 


926 0°06459 
Bor o'06292 
691 ©0°060147 


06024 


a 


05406 05970 
0573! 
05534 
05 369 
05230 
OSIIS 
OSOls 
04912 
045341 
04755 
04742 
04088 


05195 
04954 
04712 
04409 
04215 
©3939 
03633 
03240 
02753 


coooeoeoo YY OOK KK = 


~ 
-~ = C 


02213 
01544 
OO551 

"00304 

"00045 
00000 


Gouduoundouownd 
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5. The external field.—This is given from (5) with r, = 4: 


’ 5 JF Db* 
i fis + — L,+ my L, + Ofes*) >, 


we ee fy'5(e’ +2 
ee aan =p) _ edly te’ + Fe"), 


D= 3 P,= pd{r'(e* + “«)}; 
2p ) 


a. £ 
2pb? | y 


whence by direct calculation we have the following values : 


J Js Js J 

J =0°01569(1 + 0014), 0°01853(1 + O'014), 0°'02193(1 + O'014), 0°02627(1 + 0°06) 

D =0°00177(1 + 0°023), 0°00217(1 + 0°023), 0°00267(1 + 0°023), 0°00334(1 + 0°04) 

Now from (7), (8), with surface values, using (29), 
J = e',— m+ te", — Ame, 

Pe ae (30) 

D= ie*, yme's + 4Ks 

whence we have 


J, Je J, Ja 
J =0°01575(1 + O'O12), 0'01860(1 + O'012), 0'02187(1 + O'012), 0'02603(1 + 0°05) 
D = 0°00184(1 + 0°022), 0°00225(1 + 0°022), 0°00271(1 + 0022), 0°00332(1 + 0084) 
Weighted, these give 


Jy Je J J, 
J = 0°01§72 + 0°00020, 001857 + 0°00024, 0°021g0 + 0°'00028, 0°02613 + 000131 
D = 000181 + 0°00004, 0'00221 + 0°00005,0°0026g + 0°00006, 0°00334 + O’00015 
Dj J* 7°324+0°35 6-41 +031 5°61 + 0°27 4°89 40°71 
It is noticed that D/J* varies approximately linearly with J~'. Equation (1), 
with r' =b, ¢' =0, gives 
b= ag(1 — he’, — e*, — (31) 


giving, in astronomical units, 


J, Je J 
10h = 4°677 + 0°01! 4°073 + 0°01! 4°068 + ovor! 4°662 + o011 
10°/ b* = 3°439 + 0088 4°055 + O°104 4°772 + O'122 5°078 + 0-310 
10'* Dh* = 0°866 + 0°043 1°054 + 0°052 1°277 + 0°063 1577 + O’070 
6. Comparison with observation.—De Sitter collected and reduced several 
series of observations of the great satellites from 1891 to 1924, and certain older 
eclipse observations. In his 1931 George Darwin lecture* he gave, from the 
secular motions, the coefficients of the great inequalities and the period of the 
libration, equations for the masses of the satellites and the coefficient Jb*. A 
least-squares solution by Professor Brouwer, communicated by letter to Sir Harold 
Jeffreys, gives 
J b* = 5°028 x 10°”. 
This corresponds to a degree of central condensation intermediate between 
J, and J,, being much closer to J,._ Linear interpolation of D/J* gives 5:41 and 
hence Db‘ = 1-37 x 10° ™. 
The most significant dynamical effect of the figure of the planet is the secular 
motion of the orbit of the fifth satellite. We have, where mw, @ are the longitudes 


* W. de Sitter, M_N., 91, 706-738, tot 
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of the apse and node respectively, from the figure of the planet 


= = n{.Jb*/a® — 4J*%*/a4 + 2Db4/a") 
© = —n{db*/a? — 37%/a* + 2b /a4), 
dt . ' . ; 


H. Struve*, from Barnard’s observations between 1892 and 1905, gives the semi- 
major axis as 47°°965, and the mean daily motion in longitude as 722°-63175. 


Hence, with Jb* = 5-028 x 10°*, Db* = 1-37 x 10°", the contributions to the secular 
motions are, per Julian year, 


dw dé 

dt dt 

from Jb? go6’-61 

from J*b* 1-61 4°72 
from Db* 14°°5 - 14°5 

from satellite | 0°39 0°39 

from the other satellites 0°04 0°04 


gob”: 


Total 919 °°9 ~916°8 


A, J. J. van Woerkomt, discussing observations from 1892 to 1921 and in 1949, 
considers the eccentricity too small for any estimate of the apse motion to be made, 
and finds for the node 


do é é 
= O14 *02+0 °30. 
dt “ 3” 


The difference between the observed and calculated values would correspond to 
o’o12 in Jb* x10 or to o-21 in Db*x 10°". The former is a tenth of the 
uncertainty of the value calculated for the model J,, and the latter less than four 
times the corresponding uncertainty. The agreement offers support for the 
dependence of D)J* deduced from Miles and Ramsey’s models, and strongly 
favours the model J,. 


od 


7. The nature of « and D/J*.—-Sampson,{ in the introduction to his theory 
of the four great satellites, shows that the equation of hydrostatic equilibrium 
cannot be satisfied for a perfectly ellipsoidal body with uniform rotation unless 
it is homogeneous in density. ‘Thus for a real body, which has positive com- 
pressibility and whose density therefore increases towards the centre, there 
will be departure from the ellipsoidal form. That this departure must appear 
in the second order can be deduced from the equation for «. From (23), 


alr)= goa] pdt rer) 


' (r’ Pirye'(r,) —1'e'(r’))*dp 


r’ 


ae 
Spor’ / 
Te) for o<r’ <b, 


H. Struve, Sitz. Ber. Preuss, Akad. Wiss., 2, 790, 1906. 


. J. J. van Woerkom, Astron. Papers of the Amer. Ephem., 12, Pt. 1, 1950. 
. A. Sampson, Mem. R.A.S., 63, 14-16, 1921. 
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since p is a strictly decreasing function of r’ and the form e’ « 1/r’ for the ellipticity 
violates both the possible initial values for 7. Also, from (22), 

inn I I r w? P 4 
LUT = - 17 | pd((r’)'f) +r, 
v0 


SPo i r 


pd( fi(r')’) 
17 f(b)ps r2|  filr'dp——<|  (r’)'fdp 


3Po r=, ry 
o for o<r'<b if f>o for o<r'<h, 


Thus in turn 6,« >0, 5,« >o,..., and hence 


r 


K(r') >o for o<r' <b. (33) 


Sampson derives an inequality (n > m* in his notation) for which the following 
proof is extended to cover the case of non-zero «. 
arbitrary variable A 


F(A) =p, +2AS +2P 
r S\2 
. } pd G <I-+ 2A (? ) (e ? 
ia r ) ry; 
: , r\3 
af ore d(nsa(eY (e 
To” Jotent r 


Consider the function of the 


olr.)| rales -¢) ‘ 


= / 


r’\8 r’\2 2 2 r’\4 
( )| +a y (es -@*) »( )« dp. 
r, ry; ; y ry 

The second term is strictly positive, unless the body is homogeneous, when it 
vanishes for all r,. ‘The first term is positive or zero. For the homogeneous 
body, having « = 0, it vanishes, for given r,, for just one value of A. Hence F =o 
has no real root for A except in the case of the homogeneous body, when the root 
are equal, ‘Thus 

Pp, > S*, (34) 
with equality only if the body is homogeneous. ‘Taking surface values and using 
(29) we have 


D)J* 2 25/6, (35) 


with equality only if the body is homogeneous 

I should like to record my gratitude to Sir Harold Jeffreys for having initiated 
and guided my work, and to the Council of Gonville and Caius College for a 
Research Studentship which made it possible. 

Gonville and Caius College, 


Cambridec 
1955 April 5 





GREENWICH OBSERVATIONS OF THE HORIZONTAL 
AND VERTICAL DIAMETERS OF THE SUN 


P. J. D, Gething 
(Communicated by the Astronomer Royal) 
(Received 1955 June 17) 


Summary 


It has been claimed that the Greenwich observations of the horizontal 


diameter of the Sun are in agreement with those made at Campidoglio and 
Monte Mario in indicating a real variation of diameter with a period of about 
22 years. It is shown in the present paper that this conclusion cannot be 


justified from the data considered. A more detailed analysis of the Greenwich 


measurements of both the horizontal and vertical diameters in the period 
1915 to 1936 is given, in which corrections are applied for irradiation and 
personality differences between observers. ‘The resulting diameters show 


comparatively large variations from year to year. ‘These variations do not 
necessarily represent real changes in the size of the Sun: 


they may be due 
to variations in the mean personality correction 


No definite relationship can 


be established between the apparent variations in diameters and other 


phe nomena, such as the sunspot cyc le 


1. Introduction.— Numerous authors claim to have detected variations in the 


diameter of the Sun, whereas other investigators, such as Auwers (1), failed to 
find any variation, or to detect any significant difference between the equatorial 
and polar diameters. ‘The results of various series of observations have not 
been in good agreement, some authors finding an eight-year period, others an 
eleven-year period and others a twenty-two-year period, while estimates of the 
amplitude of the fundamental oscillation have ranged from a few hundredths 
to a few tenths of a second of arc, with extreme ranges of up to two seconds. 
Some authors (1-4) have examined heliometer measurements of the polar and 
equatorial diameters and others (§~7) the horizontal diameters attransit. Poor (8) 
has discussed measurements of photographic plates. 

The series of measurements made with the transit circle at Campidoglio 
(and later at Monte Mario) deserves special mention (6,7,9-16). An image 
of the Sun about one metre in diameter is projected onto a screen, and the passage 


of each R.A. limb is timed across seven wires by four observers using the 


eye-and-ear method. ‘The technique of obtaining four observations at each 


transit, introduced by Respighi in 1876, has two advantages: it decreases the 
probable error of a yearly mean because of the increased number of observations, 
and it enables differences of personal equation between observers to be calculated 
easily. On the other hand, the eye-and-ear observations across fixed wires 
would be expected to be inferior in accuracy to observations obtained with a 


self-recording moving-wire micrometer, and to involve greater personality 
differences 
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‘This series of observations seems to provide strong evidence for an oscillation 
with a period of 22 years and an amplitude of a few tenths of a second of arc in the 
horizontal diameter. It has been claimed (17) that the Greenwich observations 
support this conclusion, providing evidence of an oscillation of about the same 


period and phase. A comparison of the two sets of observations, given below, 
offers no support to this claim. 


The Greenwich results for the period 1851 to 1939 obtained with the Airy 
Transit Circle are given in Table |. The table exhibits the mean horizontal and 
vertical semi-diameters obtained from all the O-C diameter residuals published 
in the annual volumes of Greenwich Observations. 


It must be emphasized that 
corrections for personality have not been applied. 


Some of the data have been 
published previously (18,19), some are taken from computations undertaken 
by R. ‘T. Cullen for an unpublished investigation. The results have been 
tabulated in the form of semi-diameters to facilitate comparison with the Italian 
results. 


Taare I 
Yearly mean horizontal and vertical radii of the Sun observed at Greenwich 


Year Horizontal Number of Vertical Number of 


semi-diameter observations semi-diameter observations 


g61°46 57 g61°26 
1°74 103 
1°32 78 


106 
112 

SH 
rit 

93 
104y 
123 


I 109 
33 84 
12 104 
33 113 
33 126 

109 ; 125 
19 72 
12 108 
SI 
103 
105 
104 


132 
72 
113 
87 


OS 


73 


79 
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TasLe |-—cont 


hHiorizontal Number of Vertical Number of 


semi-diameter observations semi-diameter observations 


18KO's 961 °46 11s 961°43 121 
1°30 132 

33 34 1°26 96 

15459 1°23 1°03 115 

14g0 24 

rS9!1 

159g2 

1593 


1557 1°34 125 


1444 I 


"23 141 
15 130 
11 139 
44 172 
2 155 
24 185 
20 147 


1594 
1595 
1596 
1397 
1495 
1599 


170 
150 
179 
154 
173 
145 


[goo 
1gol 
1go02 
1903 
1904'* 
1905" 
1906 
1907 
190% 
1go9 
1gio’s 
igi! 
1gi2 
gt 
1914 
19I§ 
1916's 
19gt7 
1915 
1gt9g 
1920 
Ig2i 
1922 
1924 


1926 
1927 
1925 
1929 
1930 
Ig yl 
Ig 42 
1933 
1934'5 
1935 
1930 
1937'S 
1935 
1939 
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The values given in Table I do not agree exactly with those given by 


Miss Giannuzzi (17). ‘The means in Table I were obtained by giving each 


observation unit weight, whereas Miss Giannuzzi first formed monthly means 
and then gave these equal weight in forming the annual means. ‘The differences 
between the two sets of values obtained by these two different procedures are 


small in general. ‘The Italian results for the horizontal semi-diameter are given 


in Table XI of the paper by Gialanella (7) 

The two series of observations are compared in Fig. 1. ‘The upper graph 
shows the two series of measurements of horizontal semi-diameters and the 
lower graph shows the Greenwich measurements in the vertical co-ordinate ; 
measurements of the vertical diameter were not made at Campidoglio. ‘The 


diagram also shows the dates of sunspot maxima and minima for the purposes 


of comparison. 

An examination of the upper graph fails to show any obvious agreement 
between the variations found at the two observatories Viiss Giannuzzi, however, 
claimed that the Greenwich observations of horizontal diameter contained an 
oscillation with a period of between 22 and 24 years, agreeing in period and phase 
with an oscillation found from the Italian results I'he conclusions she reached 


from her study of the Greenwich observations were as follows: 


(1) ‘There appears to be a secular drift in the Greenwich semi-diameters 
from about 13g0 onwards, the semi-diameter dec reasing at the rate of about 
o”-o1 per year; the drift is removed before the results are analysed 

(11) A periodic analysis by a special method suggested the existence of two 
oscillations, one with a period between 22 and 24 years (mentioned above), 
and the other with a period between 6 and 8 year Periods of about 22 years 
and 8 years have already been detected in the Campidoglio results ; it is suggested 
by Miss Giannuzzi that the 22-year period represents a real oscillation of the 


horizontal diameter, since there is agreement in the phase as well as in period, 


whereas the 8-year period may be of atmospheric origin, since there 1s no agreement 
in phase. 


(iii) Although the 22-year oscillations agree approximately in period and 


phase, there is a discrepancy in amplitude, the Greenwich oscillation having an 
amplitude between o”:1 and o”-2 and 


the oscillation found from the Italian 
results having an amplitude between 0-2 and o’ 
In view of Miss Giannuzzi’s results 


, the following tests were applied to the 
results plotted in Fig. 1: 


(i) Correlation.—The correlation coefficient between the horizontal semi 


diameters for the years 1876 to 1936 and the Italian results for the same period 


is —0o°26, or —oo7 after removal of the secular drift from the 


(sreenwich 
observations 


(ii) Analysis for period.—I\Inconclusive results were obtained from a normal 


periodogram analysis. ‘The horizontal semi-diameters showed a suggestion of 
a seven-year period, and the vertical semi-diameters showed a suggestion of a 
ten-year period, but the evidence was not convincing in either case 


(iii) Amplitude of variations.—\t is clear from a cursory inspection of Fig. 1 


that the variations in the Greenwich semi-diameters are considerably smaller 


than those in the Campidoglio results ‘J he standard deviation of the ¢ ampidoglio 
measures from their mean is nearly three times the standard deviation of the 


Greenwich measures from the Greenwich mean. ‘This result is in reasonable 
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agreement with Miss Giannuzzi’s conclusion given as (iii) above, but it is important 
evidence against the hypothesis that the variations represent real changes in the 
Sun’s size, which would be expected to affect both series of measurements equally. 

In the above discussion of the Greenwich observations, no account has been 
taken of the differing personal equations of observers and the differing distributions 
of the dates of observations from year to year; apparent changes in the annual 
means might arise from either of these two causes. ‘The personal equations are 
systematic differences between observations of diameter made by different 
observers. Clearly, if an observer with a personal equation which is large 
compared with the mean of other observers’ personalities obtains a large number 
of observations in one year and a small number in the next year, the mean semi- 
diameter will show a consequent change which does not represent a real change 
in the size of the Sun. Similarly, the annual term in the measurements made 
by each observer, i.e. the “ irradiation’’ term found by Cullen (20), can give rise 
to apparent changes in mean diameter if the distribution of observations changes 
appreciably from year to year. 

2. The Greenwich observations from 1915 to 1937.—-\n view of the importance 
of personality corrections and of the irradiation term, it was decided to attempt 
an analysis in which allowance could be made for their effects. Apart from 
these difficulties, the transit circle appears to be particularly suitable for the 
measurement of the horizontal diameter since the quantity measured, effectively 
the time interval between the passage of the two vertical limbs past a fixed point, 
is not affected by instrumental changes. ‘The measured vertical diameters, 
however, could be affected to a small extent by changes in the scale values of the 
ZD micrometer and of the microscope micrometers used in obtaining the circle 
readings. 

The observations obtained by five regular observers with exceptionally long 
periods of service between the years 1915 and 1936 provided the most suitable 
data. ‘The introduction of the impersonal R.A. micrometer in 1915 forms 
a suitable starting point, since the personality of each observer in measuring 
the R.A. limbs probably changed with the new method of observing; it is 
possible that the personalities for measurements in the vertical direction 
also changed because of the consequent changes in observing technique. ‘The 
total number of observations obtained by the five observers showed a sharp drop 
after 1936, owing to retirements from observing duties or to transfers to other 
instruments. ‘There is therefore a poor “ connection” with the new observers 
who replaced them. The names of the observers, the initials by which their 
observations are identified, and their periods of service as regular transit-circle 
observers are as follows: 


Mr Witchell (W), 1902-1934 
Mr Cullen (RC), 1903-1935 
Mr Acton (HA), IG12-1937 
Mr Barton (HB), 1921~-1937 
Mr Symms(LS), 1922-1937 


it was decided to collect the observations of both the horizontal and vertical 
diameters in order to compare the variations in the two co-ordinates. 

The published observations of solar diameter are in the form of residuals, 
either (C—O) or (O-C), between the observed diameter O, and the computed 
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diameter C published in the Nautical Almanac for the time of transit. ‘The 
first step was to correct each residual for the irradiation term found by Cullen (20) 
‘This is an annual term in the measured diameters, both horizontal and vertical, 
varying approximately with the secant of the Sun’s zenith distance; the Sun is 
observed large in summer and small in winter. 

After corrections for irradiation had been applied, yeariy means were formed 
for each observer. ‘These means will be referred to as “individual means”’ 
‘The yearly means were then examined for personality differences by the method 
described below. ‘The appropriate personality corrections required to refer 
all observations to a standard observer were then applied, and a mean formed 
for each year, by combining with appropriate weights the corrected individual 
means; this mean will be referred to as the “‘ combined mean”’ for that year. 
‘The results are tabulated and discussed in a later section. 

3. The irradiation correction.—Throughout this section results given refer 
to diameters, not semi-diameters; this is to facilitate comparison with Cullen’s 
results. ‘The correction applied for irradiation was —0°:55 + 0°*308€c 2, where 2 
is the Sun’s zenith distance The correction was not computed for each 
observation; instead, a mean value of 2 for each month was used, and the same 


correction applied to each observation in that month, according to Table II 


Tasie Il 


Irradiation corrections 


Month | Correction to be added Month | Correction to be added | 
to (OC) | | to (O-C) 
| 


July 
Aug 


Oct | 
Nov 31 
Dex | ‘$7 


Apr. 
May 
| June 





| 
| 

Mar | Sept 
| 


ae 
These values are taken from Cullen’s paper. ‘They were obtained from an 
analysis of the observations made by seven regular observers in the period 
1gi§-1925. Some of the observations used are, of course, contained in the 
present investigation. ‘The constant term, +0":55, was chosen so that the mean 
of all the observations considered was not altered by the application of the 
irradiation corrections. 

The present work provided a good opportunity of confirming the values 
obtained by Cullen. ‘The corrections given in Table II were first applied and 
corrections for personality obtained, ina manner to be described. ‘The personality 
corrections were then applied, and all the observations were classified according 
to month. A least-squares solution, to fit the new diameter residuals for each 
month (horizontal and vertical combined) to a curve of the form a+dhsecz, 
showed that an additional correction of +0°21 +0”-05secs was required, 
‘The constant term is of no importance; it arises mainly because the personality 
corrections in the present investigation have been referred to observer W as 
standard, whereas Cullen used the mean (O—C) of his seven observers as the 
standard. ‘The coefficient of the second term indicates that the coefficient of 
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Sec in the complete correction shou 76 Ihe probabl error ot this 
coetncient however, is about +o ob and it was therefore not consider d worth- 
while to apply additional corrections for irradiation, the effect of which would 


in any case be small in the present investigation 


Miss Giannuzzi, in a recent paper (21), has also investigated the irradiation 


term in the Greenwich observation She ci ill the observations mad 
between 1851 and 1937 and arranged them according to month he monthly 


means have not been analyse d for a sec term in her ] I yut this can be done 


very easily and yields a coefficient of + 0°:23, a‘ nediate between the two 


quoted above. 


No significant differences in the irradiation terms of different observers could 


be detected nor was there any noticeable ditte between th horizontal 


and vertical diameters in thi respect It must be emphasized, however, that 


the data are not really sufficient for detailed considerat i sc points 


Afte rr moval of the irradiatio i term the probal le noel obser ation 


is about 87 for a horizontal diameter and a ; ertical diamete 
As is to be expected thes probable errors are smatlet 


value (about 1 } 
for both diameters before correcting for irradiation 


j. Personality corrections Che difference between the individual mean 
semi-diameters obtained by two observers in the same year was assumed to be 
due to personality differences and a weighted mean of the differences for each 
pair of observers was formed. Comparison of the weighted mean with the 
individual differences showed that the first on two differences obtained with 
each ne observer te nded to show ather al cattet pre umably due to 
inexperience It was therefore decided to revise the personality olution by 
omitting the first tw years’ observations of each new observer and the first two 


years’ observations of every observer after the introduction of the impersonal 


micrometer, te. no observations mad prior to 1gi7 were kept 


Taste Ill 


Personality rrectt to semi-diameter 


Horizontal pr H Vertu 


Ihe per onality correction sgiven by the we rohte dmeansare not self consistent 
becaus f the effects of weighting For example, the (horizontal) value 
W R¢ 24 and W HB 0-02 would suggest RC-HB 0°27 
whereas the actual value is 0 28 \ compromise olution wa therefore 
obtained by the method of | ast-squares The results are hown in ‘Table ITI, 
which gives the correction required to refer an observed horizontal or vertical 
semi-diameter to W as standard 

Che personality corrections obtained by this procedur independent of 


any changes in the Sun’s diameter from year to yeat vever, questionable 
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Taste IV 
Horizontal semi-diameters corrected for personality 


Combined 
Year RC HA HB LS mean 


1917 g60°54 
°° 
° 


IQis’s 
1gt9 
1g20 
Ig2i 
1922's 
1924'S 
1924 
1925 
1 g2f 
1927°4 
1926 o’60 
0°54 
0°45 
1°37 
1°22 
1°04 
o°96 


1929 
i194o 
193! 
1932 
1933" 
1934 
1935 
1930 


116 
0°35 


TaBLe V 


Vertical semi-diameters corrected for personality 


Combine 


Year RC HB 


mean 


1917's 960°sS9 


g6o 32 
1915 


0°97 
1gI9's 1°O8 1°02 
1920 1°0o 0'95 
1g2t 93 
1922 
1923 
1924 
1925 
1926 


04 

gf ft 

gO1°o4 04 

1°37 g61°26 20 
1°14 o's2 
o'S2 1°34 


1927's 
1g25°s 
1929 


O92 1°35 
1°32 
0°55 


o's! 17 


1930 
Ig yt 
1932 4 


1°45 
56 1°36 36 
1933 25 1°40 *30 1°45 ‘30 
1934 : 45 1°44 28 1°56 40 
1935° - 1°33 1°35 26 1°60 1°33 
1936 _ 1°89 1°35 1°66 1°55 1°59 
whether personality differences remain strictly constant from year to year, and 
whether the mean personality of all observers is constant. If the personalities 


of all observers are subject to changes of the same sign the combined means will 


obviously be affected. ‘The apparent variations in the solar diameters therefore 


include the effects of possible real variations and changes in mean personality. 
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The results obtained by each observer, corrected for personality, for the 
horizontal and vertical semi-diameters are tabulated in Tables IV and V. ‘The 
combined mean for each year has also been tabulated for the two co-ordinates. 
It must be emphasized that only constant corrections have been applied for 
personality ; no attempt has been made to remove slopes or “ waves’’ in differences 
between two observers. Changes in personality peculiar to one observer show 
up as a systematic drift of his results from the mean. For example, the measure- 
ments of horizontal semi-diameter obtained by RC show a tendency to increase 
over the entire period covered by his observations. 

An indication of the magnitude of changes of personality to be expected can 
be obtained by dividing the data in these tables into groups. For example, we 
can compare the mean personality differences of the observers W, RC and HA for 
the years 1919°5 to 1926°5 with those for the years 1927°5 to 1934°5. ‘The results 
are given in Table VI. _ It will be seen that there are long-term changes of two 
or three tenths of a second in personality differences, and the possibility of 
comparable changes in combined mean semi-diameters, due to personality changes, 
cannot be excluded. 

Taste Vi 


Personality diffe rences tn two perto ds of emht vears 


| Horizontal semi-diameter Vertical semi-diameters 
| 
W-R¢ W-HA W-R¢ | W-HA 
1919°5 to 1920's | o'19g 
1927°5 to 1934°5 0°35 


5. Discussion of results.—The mean values for the horizontal and vertical 
semi-diameters have been meaned in pairs to give points corresponding to the 
beginning of each year, ‘These points are plotted on a larger scale in Figs. 2 and 3 ; 
big. 4 isa plot of mean sunspot areas for comparison. 

The variations in these values are rather greater than can be accounted for 
as random scatter of the observations; however, it will be shown that they could 
be explained as variations in personality and do not necessarily represent real 
variations in the size of the Sun. 

Working from the probable errors of +0°°87 for a horizontal diameter and 
+ 0-63 for a vertical diameter obtained above, and ignoring possible variation 
in personality, we should expect the probable error of one of the plotted points 
to be +0”-05 at most. ‘The total ranges in semi-diameter are rather larget 
The greatest difference in the horizontal semi-diameters is between the years 
1923°0 and 1933°0, the values and probable errors being as follows : 
1923°0: 60°68 +0"-05 
1933°0: 960"-98 + 0-03. 


Thus the range is nearly four times the sum of the probable errors. For the 


vertical semi-diameters the extreme values are 
19180: g60"-89+0"-02 
1936°0: g61"-48 ; °” 04 


and the range is nearly ten times the sum of the probable errors 
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It must be emphasized that the probable errors quoted have been deduced 


from the internal accordance of each observer's measures in any one year; they 


do not, therefore, include the effects of possible random and systematic errors 


in the personality corrections applied. 


The effect of uncertainties in the person- 
ality corrections can be studied by considering the internal agreement of the 
diameters obtained by different observers in the same year. 


For example, 
working from the scatter of the individual means for each year, we can estimate 
the standard deviation from a constant value (nearly enough the mean of all 


observations) to be expected in the yearly means for all observers combined, on 


corrected vertical semi-diameters 


Relative sunspol area 
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the assumption that there are no real variations in the diameter of the Sun Ihe 
problem of weighting is awkward, since the probable error of an observer's mean 
tor a year 1s only partly dependent on the number of observations obtained by 
him; any systematic error due to a departure of his personality correction for that 
year from its mean value over several years will not depend on the number of 
observations. Since the systematic errors seemed likely to be of more importance 
than the random errors, it was decided to give equal weight to each annual mean, 
independent of the number of observations. 

‘The standard deviations calculated in this way from the discordances between 
annual means for individual observers were: 


o”:26 for the horizontal semi-diameters 
o0”:20 for the vertical semi-diameters 


The actual standard deviations from a constant value of the combined annual 


means are. o”-22 for the horizontal semi-diameters 


0”: 38 for the vertical semi-diameters 
‘These values show that, if anything, the combined yearly means for the horizontal 
semi-diameters are more nearly constant than would be 


expected his 
conclusion is the opposite of Miss Giannuzzi's he used 


the observations 
obtained by a larger number of observers without applying personality corrections, 
on the assumption that changes in the mean personality corrections would then 
be sufficiently small to be ignored. ‘The present 
assumption may not be justified. It can be 


analysis suggests that this 
seen from ‘lable IV that the scatter 
of the individual means, corrected for personality, about the combined mean i 


too large to inspire any confidence in the reality of the variations from year to year 


With the vertical semi-diameters the position seems to be quit ditferent 


the means seem to show larger variations than would be expected Che main 


contribution to the comparatively large calculated standard deviation comes from 


the progressive change in the means An inspection of either Table V or Fig. 3 
shows fairly clear evidence of a positive slope in the yearly means 


‘J his slope 


appears to be confirmed by most, if not all, observer rhe slope is equivalent 


to an increase of about 0°03 per year in the semi-diameter, i.e. about three times 


that found by Miss Giannuzzi in the horizontal semi-diameters, and of the 
opposite sense. ‘There is no evidence to indicate whether the drift should be 


interpreted as a real change in the diameter of the Sun, as an instrumental or 


atmospheric effect, or as a drift in the mean personality 
There does not appear to be any marked drift in the uncorrected vertical 
emi-diameters meaned for all observers Indeed, the 


uncorrected result 
tabulated in ‘lable | for both co-ordinate 


behave rather differently in the period 
1917 to 1937 from those in ‘Tables 1V and V for the same period 


f he uncorrected 
horizontal semi-diameters for the dates corre ponding to the extremes of the 


corrected horizontal semi-diameters are 
1923 '#) g6o’ 59 
1933°0: 960°°77. 


Che change 1s of Opposite sign to that in the corrected result lor the uncorrected 


i 


vertical semi-diameters the results are 
Igis’o G60’ 97 


19026 O 


gb! ob 
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‘The change is of the same sign as that in the corrected results, but is much smaller. 
The differing shapes of the uncorrected and corrected curves demonstrate 


the importance of personality corrections. It seems probable that at least some 


of the difference can be explained by the fact that two new observers with large 
personalities of the same sign in both co-ordinates started observing in 1923 and 
1925 respectively. Their observations have an increasing effect on the uncorrected 
means from these dates onward. 

An inspection of Figs. 2, 3 and 4 shows that there is no obvious connection 
between sunspot activity and the apparent variations in diameter. 
in fact, no significant correlation between the two phenomena. 

It must be 


There is, 


concluded that the Greenwich observations considered here do 


not provide any clear evidence of the existence of real variation in the size of the 
Sun. 


Royal Greenwich Observatory, 
Herstmonceux Castle, Sussex 
1955 June 16 
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ATOMIC OSCILLATOR STRENGTHS AND EXCITATION 
POTENTIALS 


i I/len and A. S. Asaad 
(Received 1955 October 


Summary 
The ratio between calculated and observed oscillator strengths is 
discussed for atoms of the astrophysically important iron group (Ti to Ni) 
The ratio is found to vary systematically and very greatly with excitation 
potential, being high for small E.P. Reasons are advanced for believing 
that observed values may be preferred to calculated values. ‘The data lead 


to the surprising result that there is a general increase of true oscillator 


strength (factor of 1000 or more) for observed lines as E.P 
or 4 volts 


increases by } 


\ knowledge of the oscillator strengths of spectral lines is essential for deriving 
the abundance and excitation conditions of atoms in astrophysical sources 
Estimates can be made theoretically or by observation and in the simpler atoms 
these usually agree. ‘This paper discusses the case of the more complex atoms 
of the iron group (Ti to Ni) which are of very great astrophysical significance 

The theoretical approach finds a very general expression in a paper by Bat 
and Damgaard (1), which provides rules and tables from which oscillator strengths 
can be calculated. ‘The method was developed primarily for simple atoms and 
was shown to give reasonable results for systems containing one or two electrons 
outside closed shells. As far as we are aware no attempt has been made to 
examine whether the method can be applied to more complex atoms. ‘There is 
no theoretical justification for doing this, but we had expected that the results of 
a rule-of-thumb application of the procedure might at least bear some useful 
relation to the true value 

For comparison with theoretical calculations one needs reliable measure- 


ments of oscillator strengths which should preferably be on an absolute scale 


‘The measurements used for the present work are by King (2, 3) and Kopfermann 
and Wessel (4) for Fe; Hill and King (g§) and Davis, Routly and King (6) for Cr; 
King (7) and Estabrook (8) for Ni; Huldt and Lagerqvist (g) for Mn; King (2) 
for ‘Ti, and King, Parnes, Davis and Olsen (10) for Co. Dzyfferences of the 


order of a factor of three exist between sets of absolute measurements, so that 
satisfactory absolute accuracy cannot be considered established. However, the 
accuracy is quite sufficient to demonstrate the part played by the excitation 
potential as illustrated in Figs. 1 and 2 

In order to compare theory with observation we use the difference log wf 
(theor) — log wf (obs) where / is the oscillator strength and m the statistical weight 
The difference has been plotted against excitation potential (upper) in Fig. 1 
The theoretical values have been computed from the Bates-Damgaard tables 
using the mean excitation potentials of the multiplet to derive the « values 


* Received in original form 1955 August 19 





C. W. Allen and A. S. Asaad Vol. 115 


cn"? 

»/+ 

by Moore (11, 12), Goldberg (13) and Allen (14) are used in the calcu- 
Observations are available on an absolute scale for Fe 1, Crt, Mn1t and 


‘Table 


lations 
Nit which are plotted in Fig. 1(a) (the King (3) absolute evaluation of Fe 1 is 


lor Cot and ‘Tit values are relative only and the data taken from the 


used) 
tabulation (2, 10) are plotted in Fig. 1(b). 
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ic. 2 Comparison of observed oscillator stre nth wt EP 
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Intercominnation lines 


The most surprising result evident from Fig. 1, is the approximate linear 
log wf (obs) and E.P. The systematic change 
It is seen that, relative to the observed 


and down to 


relation between log a/ (theor) 


with E.P. is far greater than the scatte: 
values, the theoretical values are up to so times larger for small E.1 
There is thus a total error of over 1000 in the 


30 times smaller for large E.P 
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theoretical values as compared with the observed. ‘The slope of the curve is so 
constant as to suggest that some temperature error has crept into the observations, 
but the correction is impossibly large. It would imply that the true temperature 
of the King furnace was nearly twice as great as estimated. 

Since both wf (theor) and wf (obs) in Fig. 1(a) are intended to be absolute, 
the graph would indicate agreement by a horizontal curve at zero ordinate. We 
find that neither the high E.P. nor the low E.P. lines are in agreement but that 
agreement comes arbitrarily at a medium E.P. of 4°8 volts. There is no ten- 
dency for the curve to flatten out at this level. 

Fig. 1(a) shows that the absolute evaluations of the four elements are in 
reasonable agreement with one another. ‘Two of the elements, Fe 1 and Cr 1, 
show the E.P. effect with approximately the same slope, while Mn 1 and Nit do 
not have sufficient range to show this effect. 

The E.P. effect can be demonstrated just as readily from relative measure- 
ments of wf. Fig. 1(b) shows that for Cot and ‘Ti 1 the slope of the curves 1s much 
the same as in Fig. 1 (a). Within likely error the slope for all the curves is 1-0 
(using 10-based logarithms and E.P. in volts). By comparing Figs. 1(a@) and 
(b) and assuming the curves to coincide we obtain the following conversion factors 
by which the relative values should be multiplied to make them absolute 

for Cot (2): 17x10 
for Tit (10): 85x10, 

That the E.P. effect is not caused solely by the application of the Bates 
Damgaard procedure may be seen by considering the application of the /-sum 
rule. This rule is equivalent to the statement that the strongest lines asso- 
ciated with each level have oscillator strengths of the order unity. In Fig. 2 
we plot simply log wf (obs) against E.P. In this case the scatter is much larger 
because of intrinsic differences in line intensities. Absolute values are in 
Fig. 2(a) and relative in Fig. 2(b). We see that there is a systematic tendency for 
mf (obs) to increase with E.P. with about the same slope as in Fig. 1. ‘This means 
that log wf (f-sum theor)—log wf (obs) is much the same as log wf (Bates 
Damgaard theor)—loga@/ (obs). No doubt a similar result would be obtained 
by using any existing method of calculating oscillator strengths for more complex 
atoms. It is worth noting from Fig. 2 that the intercombination lines show the 
effect as well as the lines allowed in LS coupling. 

With such large discrepancies between theoretical and observed values one 
has to face the quandary of which to regard as right and which wrong. It is get 
likely that any compromise would be correct. Fortunately there are three con 
vincing arguments in favour of the correctness of the observed values: (a) a 
furnace temperature error large enough to explain the result is unthinkable, (4) 
using theoretical values with some copper arc measurements now in progress 
would yield an impossible excitation temperature of zoooodeg.K, and (c) 
independent measurements of lines by emission (15) as well as absorption all 
show the effect. 

If we conclude that the observations are correct, we are led to the following 
rather surprising results: 

(a) Oscillator strengths for lines of high E.P. may have values well above 1o 

(6) Oscillator strengths of lines of low E.P. (the strongest permitted lines with 
LS coupling) have values of the order o-o1, which contrasts strangely with the 
weil-established D lines of Na, Ca 1 4227, and the H and K lines of Ca 1, 


/ 


40 
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(c) The existing theories give no a priori indication of absolute oscillator 
strength since the E.P. of 4-8 volts where the values appear to be correct is entirely 
fortuitous. We are therefore entirely dependent on observations for estimates of 
absolute oscillator strengths. 

A partial explanation of the E.P. effect may be that, as the E.P. increases, 
more and more electrons contribute to s, the effective number of electrons. Also 
for the absorption oscillator strengths of high-level lines the negative emission 
components in the f-sum rule Uf), =s 


~ fem May add progressively to give a 
high value. 


This however would not explain the surprising regularity and 
generality of the effect, nor the reason why the low-level lines should have such 
small oscillator strengths. 


One of us (A. S. A.) expresses his gratitude to the Egyptian Government for 
sending him to London University to carry out research in spectroscopy. 


University of London Observatory, 
Mill Hill Park, 
London, N.W.7 : 
1955 August 17 
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NOTE CONCERNING “THE CONSTANTS OF THE VELOCITY 
ELLIPSOID FROM THE RADIAL VELOCITIES OF 820 STARS” 


D. G. Ewart 
(Communicated by the Director, University Observatory, Glasgow) 
(Received 1955 June 3) 


Summary 


Erroneous values of the solar motion and of the constants of the velocity 
ellipsoid published in a previous paper are corrected. The revised results 
indicate no definite variation in the constants of the velocity ellipsoid for 


the faint stars, but the variation in the solar motion is confirmed 


1. Ina previous paper (1), the derivation of the constants of the solar motion 
and of the velocity ellipsoid from the radial velocities of a sample group of faint 
stars—HD magnitudes 8-5-8-6, spectral types F to M—obtained at Lick 
Observatory by Moore and Paddock was described. It was recently brought 
to my notice by Dr T. R. Tannahill that the solution for ‘“‘ALL” stars was 
incorrect, the figures given not being the means of the values for the three zones. 
On examination, it was found that it had been assumed that the equations could 
be solved for the whole sky by taking the mean galactic latitude as 0°, which 
procedure is fallacious 

A new solution was consequently required, and it was decided to re-chech 
the previous solutions simultaneously. In so doing, two errors of sign were 
discovered in the solution for the solar motion. A fresh analysis was therefore 
performed, the details of which are given below. 

2. No changes were made in the methods of analysis. For the solar motion, 
the normal equations were formed, and for the calculation of the areal corrections, 
their approximate solutions found to be: 


U = 21-05 km sec, G= 41°, g= + 22° with a K-term of + 1-0okm/sec. 


The complete solutions of the normal equations—including the corrections 


were found to be: 
l 22°31 + 1:20 km/sec, A = 280-0 + 5°'8, D +47°-9+2°°8, 
G= 44°33 '5,8 +20 °4+3°1, anda K-term of + 1°17 + 0°79 km/sec. 


] 
The weighted mean results of Moore and Paddock are [ 22'44km sec, A = 280", 
D= +44°°5. 

The constants of the velocity ellipsoid were again determined separately for 
the latitude zones 0°, 20° and 45°. For the whole sky solution, the iterative 
method of solution described by Smart and Chandrasekhar (2) was used, the 
initial values being the means of the zone results. ‘The results, together with 
those of the previous solution and those of Smart's analysis (3) of the radial 
velocities of Schlesinger’s Catalogue of Bright Stars, 1930 are given in ‘Table | 
It must be noted that the values for the whole sky from Smart's analysis are mean 
of the three zones—-weighted in the case of the longitude of the major axis. In 
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the table, the first line gives the revised values of the constants, the second line 
the “old”’ values and the third line Smart’s values, except for the case of “‘ ALL”’ 
stars, where the doubly erroneous “‘ old’’ values are not listed. 


Tase | 
No. of 


Zone : /K WH K/H 


stars 


° 167 37°78 32°24 0°853 298°: 
38°89 32°50 0°337 292 °*' 
(619) (33°70) (20°90) (0619) (340 ° 
277 39°63 4 31°73t2°5 o'Bo1 + 0°159 27 
41°45 26°11 0°630 324 
(931) (30°10) (26°40) (0-879) (332 
296 §3°03+5 19°17 + 7° 0°3577t 0'165 344 
§3°52 17°94 0°336 343 °°7 
(10638) (37°20) (22°10) (0°594) (338 °8) 
820° 40°95 + 4°61 27°742 “ 0°677 + O'117 331 4704 
(2668) (33°70) (23°10) (0°697) (339 ‘5 1 *4) 


* This solution includes the zones centred at b= + 69° for which no separate solution was 
made. 


3. As was noted in the previous paper, the only solution to which any weight 
can be attached is that for ‘‘ ALL” stars. Again the value of the axis-ratio of the 
velocity ellipsoid agrees well with Smart’s determination, but the velocity 
dispersions along the axes are greater for the faint stars. ‘This is not unexpected, 
both because of the lesser amount of the material available for analysis and because 
of the greater probable errors of the solar motion. For the longitude of the 


major axis of the ellipsoid, the result is inconclusive because of the large probable 
error, as regards the vertex. 


Summarizing the revised results, it appears that although there are no definite 
indications of variation in the ellipsoidal constants between the brighter and the 
fainter stars-—-excepting the increased dispersions along the axes—there 


a variation in the solar motion, in velocity, and in the longitude of the apex. 


i8 


L/meversity Observatory, 
Glasgow, W.2 


19gss May 24 
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MEAN DAILY AREAS AND HELIOGRAPHIC LATITUDES 
OF SUNSPOTS IN THE YEARS 1953 AND 1954 


Royal Greenwich Observatory 


(Communicated by the Astronomer Royal) 


(Received 1955 September 27) 


The following results are in continuation of those for the years 1951 and 1952 


in M.N., 114, 606, 1954. They are derived from the measurement at 
Herstmonceux of photographs taken at the Royal Greenwich Observatory, the 
Royal Observatory, Cape of Good Hope, and at the Kodaikanal Observatory, 
India. 

Table I gives the mean daily areas of umbrae, whole spots and faculae for 
each synodic rotation of the Sun included in the years 1953 and 1954. ‘The 


‘TABLE | 
Mean daily areas 
Projected" (Corrected 
Notation Days 
commenced photo Whole 
thy y graphed Umbre spots 
1952-3°-4 
Dec 
Jan 
eb 
Viar 
\pr 
Apr 
Viay 
Jume 
Jul 
Aug 
Sept 
(et 
Nov 
Dex 
Jan 
Jan 
Feb 
Mar 
Apr 
May 
June 
July 
Aug 
Sept 
Sept 
Oct 
Nov 
* Expressed in millionths « Sun’s dish 


t Expressed in millionths Sun's hermaphere 
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numbers of the rotations are in continuation of Carrington’s series. ‘The annual 
means for the two years are included in Table I1, which summarizes similar 
figures for the years since the last sunspot minimum of 1944. As 1954 was also 
a minimum year Table II covers the whole of the most recent 11-year cycle. 


Taare Il 
Mean daily areas 


sate _ 


No. of days Projected * 


Aer, = A 


Corrected for foreshorteningt 


r 


» Photo- Without , Whole . ; Whole : 
Year Umbrae Faculae Umbrae Faculae 
graphed spots spots spots 


1944 366 7 30 160 284 23 126 
1945 365 é 102 560 


344 
774 75 429 94° 
1940 395 ? 359 2351 1794 2g! 1517 2155 
1947 365 u 555 3559 2326 405 2637 2594 
1945 366 419 2618 1849 314 1977 2331 
1949 305 ) 435 2873 2140 350 2129 2597 
1950 305 274 1633 142! 203 1222 1754 
1951 305 2 259 1552 1145 155 1136 1379 
1952 366 24 g6 545 594 404 7il 
1953 305 145 15 197 275 2 146 331 
1954 365 237 9 49 117 35 135 
* Expressed in millionths of the Sun's disk 


t Expressed in millionths of the Sun’s hemispher« 


lable III] gives for each rotation in the years 1953 and 1954 the mean daily 
area of the whole spots (corrected for foreshortening) and mean heliographi 
latitude of the spotted areas separately for the northern and southern hemispheres 
of the Sun. ‘The mean heliographic latitude of the entire spotted area and the 
mean distance of all spots from the equator are also tabulated. Corresponding 
annual mean values of all these quantities for the years 1944 to 1954 are given 
in ‘Table IV. ‘The years of minimum, 1944 and 1954, are each given twice in 
‘Table IV: once for all spots occurring during each of the years and once with 
distinction between spots of the old cycle, appearing at low latitudes, and those, 
at high latitudes, of the new cycle. 

Tables If and IV are in continuation of similar tables in Monthly Notices 
for the years 1874 to 1888, 49, 381, 1889: 1889 to 1902, 63, 465, 1903: 1901 
to 1914, 76, 402, 1916: 1913 to 1924, B§, 1007, 1925: 1923 to 1933, 94, 870, 
1934 and 1933 to 1945, T1O0, SOI, 1950 

lhe chief features of the record for the years 1953 and 1954 are as follows: 

(1) Throughout 1953 sunspot activity continued to decrease, reaching a 
minimum during the first half of 1954. ‘The mean daily area of sunspots during 
1953 was less than half that for 1952 and the mean daily area for 1954 was less than 
half that for 1953. ‘The most prolonged periods of low sunspot activity occurred 
in January-February and April-June of 1954. ‘The exact epoch of sunspot 
minimum may be put as 1954°3. High-latitude spots of the new cycle were 
first recorded in August 1953 but were not frequent until a year later. 

(2) The largest spot group of 1953 crossed the Sun’s central meridian on 
April 28:2. The mean latitude for the disk-passage was 10-9 N. and the mean 
area 839 muillionths of the hemisphere. 





1944 
1944 


1945 
1949 
1947 
1945 
1949 
I1gs§o 
19s! 
1982 
1953 
1954 


1954 


of sunspots in the years 1953 and 1954 


Tase Ill 


Spots north Spots soutl 
of the equator of the equator 


Me an 


Mean Mean latitude 

Rotation Mean hel Mean Hel 
wilO- eno 
commenced daily : daily , 
oT vraphic vrapnic spotted 
UT area i area | ' 
latitude latitude 


ft entire 


irca 


1982 
Dex 
Jan 
Feb 
Mar 
Apr 
Apr 
May 
June 
July 
Aug 
Hept 
Oct 
Nov 
Dec. 
Jan 
Jan. 
Feb 
Mar 
Apr 
May 
June 
July 
Aug 
Sept 
Sept 
Oct 
Nov 


Spots north Spots south 
of the equator of the equator 
Mean 


Mean Mes Me an Mean latitude 
Days 


heli helio o 
daily s daily f entire 


photographed graphic graphic potted 
area area 


latitude latituce area 


366 


bs Old Cyele 
- New Cyck 


395 


365 
305 
366 
365 
395 
365 
366 
365 
365 
(id Cyck 


305 New ¢ y< le 


Mean 
distance 
from 
equator 
of all 
spots 


Mean 
diatance 
from 
equator 
of all 


spots 
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The C.M.P. of the largest spot of 1954 was March 17-1. Its latitude was 
8°-2 5. and mean area 403 millionths. ‘The largest new-cycle spot during 1954 
oceurred in latitude 33°°7 N. and had its C.M.P. on December 15-5; its mean 
area was 317 millionths. 

(3) On 148 days in 1953 and 237 days in 1954 no sunspots could be seen on 
the Sun, while on 69 days in 1953 and 180 days in 1954 faculae were also absent. 
January 1954 and May 1954 each had only two days with spots present. 

(4) The ratio of mean corrected areas of faculae/sunspots was 2°27 in 1953 and 
39 in 1954. ‘The ratio of mean corrected areas of umbra/whole spot was 0-171 
in 1953 and o-171 in 1954. 

(5) In 1953 the northern hemisphere was more spotted than the southern, 
while in 1954, when new-cycle spots first predominated, the reverse was true. 
In 1954 the mean distance from the equator of all new-cycle spots (28°-96) was 
the highest similar figure in the Greenwich records, the previous highest being 
27°-6 for the new-cycle spots in 1934. ‘The old-cycle spots of 1953 and 1954 
were, on the average, at approximately the same latitudes (8°~g°) as the spots of 
1952. 

The number and distribution, northern and southern hemispheres, of spot 
groups of 

(a) two days’ duration or longer, 

(6) one day’s duration 
were as follows: 


1954 
(b) (a) (b) 
Northern spots 35 iI 20 10 
Southern spots 13 13 8 


Total 61 24 33 18 


(6) The following table gives the mean daily areas of sunspots and of faculae, 
corrected for foreshortening and expressed in millionths of the Sun’s hemisphere, 
for each calendar month of the years 1953 and 1954. 


1953 1954 
Month Spots Faculae Spots Faculae 
Jan 259 577 
Feb 10 381 
Mar 7° 420 
\pr 406 324 
Viay 196 314 


54 
55 
130 
70 
33 
19 
50° 
188 


wa 


ce - NK BS = 


June 215 370 
July BI 280 
Aug 269 351 
Sept 141 416 144 
Oct 47 299 379 
Nov 9 182 55 295 
Dec 6 64 106 223 


Vv 


= 
cS 


Finally, short descriptive surveys of the sunspots of 1953 and 1954 have been 
given in Monthly Notices, 114, 358, 1954 and 115, 185, 1955, respectively. 
Royal Greenwich Observatory, 


Herstmonceux Castle, Sussex : 
1955 September 16 
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